Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,3490CA41.5% 120.7 / 36.7 65V69F (GTT→TTT) yafYlipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (6/32);  new base A (27/0);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 9.04e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAT  >  NC_000913/266208‑266494
                                                                                                                                                |                                                                                                                                                   
tCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAgg                                                                                                                                             <  2:1314079/150‑1 (MQ=255)
      cctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGAc                                                                                                                                      <  1:799178/150‑1 (MQ=255)
         ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGa                                                                                                                                  <  1:566832/151‑1 (MQ=255)
         ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACAGCAAGGGGGGGcc                                                                                                                                     >  2:821058/1‑148 (MQ=255)
          gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGAt                                                                                                                                 <  2:747824/151‑1 (MQ=255)
           tATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAggggggggcgg                                                                                                                                    >  2:305706/1‑143 (MQ=255)
            aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATa                                                                                                                                <  1:1312052/150‑1 (MQ=255)
            aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGGAAAAACCGCAAAGGGGGGCCCGaga                                                                                                                                >  1:561292/1‑148 (MQ=255)
             ttCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAGGGAAAAACCTCAAAGGGGGGCCCgggga                                                                                                                                >  1:378043/1‑147 (MQ=255)
               ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGATGTAAAAACTCCAAAGGGGGGGCGGAGaa                                                                                                                               >  1:1297816/1‑148 (MQ=255)
                cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACTGCAAGGGGGGACCGGATa                                                                                                                                >  1:763100/1‑146 (MQ=255)
                     tGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCaa                                                                                                                       <  2:76439/150‑1 (MQ=255)
                      gATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAAcc                                                                                                                     <  1:453990/151‑1 (MQ=255)
                           tATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCtgtg                                                                                                                 <  1:1169239/150‑1 (MQ=255)
                            aTCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGATAAGAAGGCAACCTGTGAc                                                                                                               >  2:220240/1‑151 (MQ=255)
                             tCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACAGCCAAGGGGGGCCGGGTAAAAAGTCAAACCGTGGcc                                                                                                              >  1:1057395/1‑151 (MQ=255)
                              cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAAGGGGGGCCGGGGAAGTAGGCAACCTGTGAccc                                                                                                             >  2:817747/1‑150 (MQ=255)
                              cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGAAAAAACCACAAGGGGGGGCCGAGTAAGTGGGGAAACTGTGagc                                                                                                              >  2:480287/1‑148 (MQ=255)
                                   aTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCAttc                                                                                                          <  2:1224372/149‑1 (MQ=255)
                                              aTTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCGCAAGGGGGGGCCGGATAAGTAGTCAACCTGTGACCATTTTTCGCTGAcatc                                                                                             >  2:596051/1‑150 (MQ=255)
                                              aTTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGGGAAGAAGTGAAACTGTGAACAttcttttctggca                                                                                                >  2:45942/1‑140 (MQ=255)
                                                tcatcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAACGGAAAACTTCCTCGGGGGGGCCGGGAAAGAGTCAAACCCTGTGCACACTTCTCTGAcatcctt                                                                                          >  1:1072962/1‑148 (MQ=255)
                                                 catcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTTAAAACTGCAAAGGGGGGCCGGGGAAATAGTCAAACCGTGACCaattttttct                                                                                                    >  1:267918/1‑132 (MQ=255)
                                                  atcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACAAAAAAAACTCAAGGGGGGGCCGGGTAAATAGTCAAACTGTGACCCTTCTTCGCTGAcataaaccc                                                                                        >  2:743448/1‑145 (MQ=255)
                                                       tGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCAGCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCtggt                                                                                     <  1:637968/149‑1 (MQ=255)
                                                               cGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACcac                                                                           <  2:71097/151‑1 (MQ=255)
                                                                 gTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAGGGGGGGCCGGGGAAAAAGTCAACCCGTGACCCTTCTTCTCTGACATCATCCTGGGGGGTGACcac                                                                           >  2:931008/1‑149 (MQ=255)
                                                                   tttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACCGAAAAACCGCCAGGGGGGGCCGGGTAAAAAGGTAAACCGTGACCATTCTTTCCCGACATCATCCTTGGGGCTGGcccc                                                                           >  2:204996/1‑145 (MQ=255)
                                                                     tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAAATGCAAGGGGGGACCGGATAAATAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGGTGAcc                                                                             >  2:410753/1‑143 (MQ=21)
                                                                      ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGTACTAAAACCCGCAAGGGGGGAACCGGAAAAAAATCAAACTGTTAACCTTTTTTTCTTTCATCATCCtgctggtctgccccaaccctt                                                                    >  2:74457/1‑136 (MQ=255)
                                                                      ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGTCTGACCACATCACt                                                                     <  1:1289039/150‑1 (MQ=255)
                                                                      ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACt                                                                     <  2:258286/150‑1 (MQ=255)
                                                                       tCCTACATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTg                                                                    <  1:76264/150‑1 (MQ=255)
                                                                        ccTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaa                                                                  <  2:1258681/151‑1 (MQ=255)
                                                                           gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCAGTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaatc                                                                <  1:177423/150‑3 (MQ=255)
                                                                            cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcgc                                                              <  1:495806/151‑1 (MQ=255)
                                                                             atatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAAGGAAAAACTGCAAAGAGGGGCCGGATAAATAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGGTGACCCCATCCCTCAAAcac                                                              >  1:1201521/1‑148 (MQ=18)
                                                                              tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGGCCGGGAAAGTAGTCAACCTGTGGCCCTTCTTTGCCTACCTCATTCTTGTGGGTGACCACATCCCTGaaaaaggg                                                            >  2:131218/1‑146 (MQ=255)
                                                                                 ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTgg                                                         <  2:1096822/151‑1 (MQ=255)
                                                                                 ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAAGGGGGGGCGGGGAAAAAGTCAACCCGTGGACATTCTTTGCTGACATCATCCCTGGGGGTGACCAAAACACTGAAAAGCGgggg                                                         >  2:1297641/1‑151 (MQ=255)
                                                                                   aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACAGCAAAGGGGGAACGGGTAAGAAAGCAACCTTTGAACATTTTTCTGCTGCAATCATCTTGTGGCTGACCACCTTACTAAAAAACGgggggg                                                       >  1:898842/1‑151 (MQ=255)
                                                                                     gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCCCACGGGGGGGGCGGATAAGAAGTCAACCCGTGAACCTTCTTCGCTGTCATCATCCTTGTGGCTGGCCACACCACTGAAAAGCGgggggggg                                                     >  1:971876/1‑149 (MQ=255)
                                                                                     gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGCGGTAAAAACTGAAAGTAGGGACCCGGAAAGCGATAAACCCCTCACCATTTTTCTCTTCGACAATCCCGGGGGCTGGCCCCATCCCTGAAAACCGGTggggg                                                      >  2:1110177/1‑149 (MQ=255)
                                                                                      aGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGa                                                    <  2:125851/151‑1 (MQ=255)
                                                                                       ggTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAg                                                   <  2:852972/151‑1 (MQ=255)
                                                                                       ggTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGACCGGAAAAAAAGTCAAACTGTGACCATTCTTCCCTGACAACATCCCGGTGGCTGACCACATCACTGAAACGCGgggggggg                                                     >  2:907149/1‑147 (MQ=255)
                                                                                          agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTc                                                <  2:139417/151‑1 (MQ=255)
                                                                                            tagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAAACTGCAAGGGGGGGCCGGGGAAGTAGGCAAACCGTGACCCTTCTTCGCCGACATCATCCCGGGGGGTGGGCCCATCACTGAAACGCGggggggggggggcc                                               >  1:1121858/1‑142 (MQ=255)
                                                                                              gtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAAGGGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGCTGACCACATCACTTAAAAGCGGGGGGGTGGGGGCCtt                                             >  2:514947/1‑150 (MQ=18)
                                                                                                aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                           <  2:635894/150‑1 (MQ=255)
                                                                                                aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                           <  1:200497/150‑1 (MQ=255)
                                                                                                     ccATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCTGCAAGGGGGGGCCGGGTAAGAAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGCGGGGGGAGGGCCTTTAAttt                                       >  2:26040/1‑149 (MQ=18)
                                                                                                       aTCCTCGAAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTg                                   <  2:1087161/151‑1 (MQ=255)
                                                                                                          cTCGCAGCCATCGGGCGACTGGCGGGTACGTCAAACCTGCAAGGTGGGGCCGGATAAGTAGTCAACCCGTGACCATTCTTTGCTGACATCATCCTGGTGGGTGACCACAACACTGAAACGCGGGGGGGGGGGGGCCTTTAATTTACt                                    >  2:504208/1‑147 (MQ=255)
                                                                                                              cAGCCATCGGGCGACTGGCGGGGACGTAAAACCCGCAAGGGGGGAGCGGATAAGAAGTCAACCCGTGACCCTTCTTCGCTTACCTCCTCCCGGTGGGTGGCCACATCACTGAAACGCGGGGGGGGGGGGGGCTTTAATTTACTGATgggc                             >  1:849572/1‑148 (MQ=255)
                                                                                                                gCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTc                          <  1:1197946/151‑1 (MQ=255)
                                                                                                                      gggCGACTGGCGGGTACGTAAAAACCGCAAAGGGGGGCCGGGTAAAGAGGCAAACCGTGGCCATTCTTTGCTGGCATCATCCCGGGGGGTGGCCCCAACACTTAAAAGGGGGGGGGGGGGGGGCTTTAATTTTCTGATgggccttagg                       >  1:1265569/1‑140 (MQ=255)
                                                                                                                       ggCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcat                   <  1:408560/151‑1 (MQ=255)
                                                                                                                           aCTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACt               <  2:1130843/151‑1 (MQ=255)
                                                                                                                             tGGCGGGTACGTAAAAACCGCAAGGGGGGACCCGGTAAGTAGGTAACCTGCGGCCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGGGGTGGGGTGTGGTCCTTAAATTTACTGATCGCCTTCCGGtcatc                  >  1:797539/1‑146 (MQ=18)
                                                                                                                              ggCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCt            <  2:428753/151‑1 (MQ=255)
                                                                                                                               gCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg           <  1:784059/151‑1 (MQ=255)
                                                                                                                               gCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg           <  1:1254423/151‑1 (MQ=255)
                                                                                                                                   ggTACGTAAAACCTGCAAGGTGTGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa        <  1:707042/151‑1 (MQ=255)
                                                                                                                                            aaaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGAGCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAt  <  1:1063727/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAT  >  NC_000913/266208‑266494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: