Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,8620GA62.8% 32.6 / 88.7 87intergenic (+31/‑72)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (6/26);  new base A (54/0);  total (61/26)
Fisher's exact test for biased strand distribution p-value = 1.26e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-12

ACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAA  >  NC_000913/2386712‑2387005
                                                                                                                                                        |                                                                                                                                                           
acCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAg                                                                                                                                                             <  2:768725/151‑1 (MQ=255)
   gCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCCCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGa                                                                                                                                                          >  1:1222769/1‑151 (MQ=255)
             gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAAGCACatt                                                                                                                                               >  1:31015/1‑151 (MQ=255)
             gTAATAAACAATACCACTTGCCCGGTGAAACCGCAGACAGATTGAATTATCAGAAAATGGCCCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAGACAtaaaaatt                                                                                                                                               >  1:1567448/1‑144 (MQ=255)
              taataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACattt                                                                                                                                              >  1:1444311/1‑151 (MQ=255)
                   aaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttattt                                                                                                                                        <  2:772487/151‑1 (MQ=255)
                     cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCCGAAAAAAAGAAATTACCattttattttt                                                                                                                                      >  2:1130533/1‑150 (MQ=255)
                     cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACaattttttttt                                                                                                                                      >  1:853465/1‑150 (MQ=255)
                          ccACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAGACAAGCACcttttatttttaccca                                                                                                                                 >  1:1304401/1‑145 (MQ=255)
                            aCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                                                               <  1:497265/151‑1 (MQ=255)
                                 cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGAAAAGCACATTTTTTTTTTCCTCATCTGttt                                                                                                                         >  2:1146980/1‑151 (MQ=255)
                                     gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAACACAATTTTTTTTTCCTCCTCTGtttttaa                                                                                                                     >  1:1182107/1‑148 (MQ=255)
                                     gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCCAAAAAAAAACAATCCCAATTTTTTTTTCCCTCCCTGtttttta                                                                                                                     >  1:829819/1‑148 (MQ=255)
                                     gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAATGAAATGCAAATTTTTTTTTTCCTCATCTGtttttaa                                                                                                                     >  1:1147023/1‑148 (MQ=255)
                                     gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATAAAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTGTGAAAAAAAGACAAGCAAATTTTTTTTTTCCTCTTCTGTTTTaaa                                                                                                                     >  1:680929/1‑150 (MQ=255)
                                      tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                                    <  1:1209308/151‑1 (MQ=255)
                                      tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                                    <  2:478918/151‑1 (MQ=255)
                                          aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTCCTCCCCTTTTTTTAAATGCg                                                                                                              >  1:1190825/1‑151 (MQ=255)
                                                gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAATGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATGCGAAAAAt                                                                                                       >  1:771105/1‑151 (MQ=255)
                                                     aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAAt                                                                                               <  2:1016554/151‑1 (MQ=255)
                                                      ttGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAAtt                                                                                              <  1:1567617/151‑1 (MQ=255)
                                                       tGAATTATCAGCAAATTGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGCGGCGGGTCCAGAAAAAGCGACATGCCCATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTc                                                                                             <  1:168474/151‑1 (MQ=255)
                                                                 gAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTCTGTTTTTAAATGGGAAAAAATAAAATTAACgggggga                                                                                >  1:676609/1‑149 (MQ=255)
                                                                    aaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaaaa                                                                             <  1:1089789/151‑1 (MQ=255)
                                                                       ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCCCATTTTATTTTTCCTCATCTGTTTTTAAATGCGGAAAAAGAAAATTTAAAGAGGaaaaaaaa                                                                          >  1:336265/1‑151 (MQ=255)
                                                                        gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGACCAGGAAAAAAGAAAAGGACAATTTTTTTTTTCTTCTCTTTTTTTAGATGAGGAACATGAAATTTCAACGGGGGaaaaaaaa                                                                          >  2:137293/1‑151 (MQ=255)
                                                                          tCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATg                                                                       <  2:1181556/151‑1 (MQ=255)
                                                                            aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAATGCACAATTTTTTTTTTCTCATCTGTTTTTAGAAACGAAAAATGAAAATTAAAGGAGGAAAAATAAATggg                                                                     >  2:499571/1‑151 (MQ=255)
                                                                            aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAAACAATTAAAATTTTTTTTTCCTTATCTGTTTTTAAAAGGGAAAAATAAAAATTAAAAGAGGGAAAAAAAAATgg                                                                      >  1:491919/1‑151 (MQ=255)
                                                                             ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGGAAAAAAGAAAAACAAAATTTATTTTTCCTCCTCTGTTTTTAAAAGCGAAAAATGAAAATTACAGGAGGGAAAAAAAAAgggg                                                                    >  2:173829/1‑150 (MQ=255)
                                                                             ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGGACATTTTTTTTTTCCTTCTCTTTTTTTAAATGCCGAACCTGAAATTTAACGGAGGAAAAATAAATtggg                                                                    >  2:460907/1‑150 (MQ=255)
                                                                             ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCACGAAAAAGGAAAAGCGCATTTTTTTTTTCCCTATCTCTTTTTATAAGCGCAAAAATTAATATCAAAGGGGAAAAAAAAAATggg                                                                     >  1:1101042/1‑151 (MQ=255)
                                                                             ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                                                               <  2:1495515/102‑1 (MQ=255)
                                                                             ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTTTTTTTAAATTCGAAAAAAAAAAATTAACGGGGGGAAAAAAAAATggg                                                                     >  1:867184/1‑151 (MQ=255)
                                                                             ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCACAAAAAAAGAACAGAACATTTTTTTTTTCCTCATCTTTTTTTTGATGCGACAAATAAAAATAAACAGGGGGAAAAAAAAATggg                                                                     >  1:1209569/1‑151 (MQ=255)
                                                                             ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTACAGAAAAAAGGCCATTCACAATTTTTTTTTCCTTTACTGTTTTTAGAAGTGGAAAAAGAAATTTAAAGGGGGGAAAAAAAATTGGa                                                                    >  1:781684/1‑151 (MQ=255)
                                                                             ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGGCAAGCAAATTTTTTTTTTCCCCCTCTTTTTTTAAAAGGGAAAAAAGAAAATTCACGGGGGGAAAAAAAAATggg                                                                     >  1:69843/1‑151 (MQ=255)
                                                                              gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGGAAAGGACATTTTTTTTTTCCTCATCTGTTTTTAGAAGCGAAAAAAGAAATTTAAAAGGGGAAAAAAAAATAGGaa                                                                   >  1:1469403/1‑151 (MQ=255)
                                                                                 gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGCGCAGAAAAAAAGACACGTAAATTTTATTTTTTCTTATTTGTTTTATGATTCGGAAAACGAAAATCAACAGGGGAAAAATAAATGGGGATGa                                                                >  2:847086/1‑151 (MQ=255)
                                                                                   ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAACGAAATGCAAAATTTTTTTTTCCTCATCTGTTTTTAAAATAGAAAAATGAAAATTAACGGAGGGAAAAAAAAAGGGGAAGAtt                                                              >  1:1460165/1‑151 (MQ=255)
                                                                                    gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAtt                                                              <  1:1355369/150‑1 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGTTTCCGAAAAAAGGGCATGCCCATTTTTTTTTTCCTCATTTTTTTTTAGATTTCAAAAAATGAAATTAACCGGGGGAAAAAAAAAGGGAATGATTgggtt                                                         >  1:218333/1‑148 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGGAAAAAAAACAAACAAATTTTATTTTCTCCCATCTGTTTTTATATGAGGAACATAAAAATCTAACGGGGGAAAAATAAATggaatgttttggta                                                          >  2:1023518/1‑139 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAAACATGAAAAATTTATTTTTCCTCATCTGCTTTAATAAGCGAAAAAAGAAAATCAAAAGAGGAAAAAAAAATAGGGAGGATTgggtt                                                         >  2:1478027/1‑148 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCCCATTTTTTTTTTCCCCATCTTTTTTTAAAAGCCAAACAAAAAATTTAACCGAGGAAAAAAAAAATGGGATGATTGGCTa                                                         >  1:336507/1‑151 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTTTTTTTAGATGGGGAAAAAGAAAATTAAACGGGGGAAAAAAAATTGGGaagatttgggt                                                          >  1:895989/1‑142 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGAAAAGCACATTTTATTTTTCCTCATCTGTTTTTAGAAAGGAAAAAAGAAAATTAAAGGAGGAAAAAAAAATGGGGATGGTTTGCTa                                                         >  2:1253946/1‑151 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAAATGACATTTTTTTTTTCCCCCTCTTTTTTTAAATGGGGAAAAATAAAATCAAAAGGGGGAAAAAAAAAGGGGAAGGTTTggtt                                                          >  1:245181/1‑149 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCAAAAAAAAGGCAAGCACATTTTTTTTTTCCCCCTCTGTTTTTAAAAGGGGAAAAATAAAATTAAAGGGGGGAAAAAAAAAAGGGAAGATTGGGTa                                                         >  1:942132/1‑151 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt                                                          <  2:37048/150‑1 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGACAATTTTTTTTTTCCCCATCTTTTTTTAAAAGCGAAAAAAGAAAATTAACAGAGGGAAAAAAAAAGGGGaagttttggtt                                                          >  1:1163924/1‑142 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCCCATCTTTTTTTAAAAGGGGAAAAAGAAAATTAAAGGAGGGAAAAAAAATGGGGATGAtttggta                                                         >  1:1494375/1‑146 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCAAAATTTATTTTTCCCCCCCCGGTTTTATAATCGGAACATAAAATTCAAAAGAGGGAAAAAAAATAGGGATGATTGGCtt                                                         >  1:1050563/1‑150 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAGGGAAATGCACATTTTTTTTTCCCCCATCTGTTTTTAGATCGGCAAAATGAAATTTAAAAGAGGAAAAAAAAAAAGGAATGATTgggtt                                                         >  1:1231080/1‑148 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCTGAAAAAAGAGAATCACATTTTAATTTTTCTTCTCTTGTTTAATATGAGAAACATGAAAATTAACGGCGGGAAAATAAAATGGAATGGTTTGCt                                                          >  2:1264079/1‑150 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGGAAAAAAGGCATTCACACTTTAATTTTGCTTCTCCTTTTTTAGAGACGGAAAAAGAAAATCATCGAAGGAGAAATTAATgggcaatgttggccg                                                         >  2:749772/1‑139 (MQ=255)
                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAGACAAGCAAATTTTTTTTTTTCTCATCTGTTTTAAGATGCGAAAAAAGAAAATCAACGGAGGTAAAATAAATGGGAATGATTggtgt                                                         >  2:521522/1‑148 (MQ=255)
                                                                                         tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTa                                                         <  2:1469811/150‑1 (MQ=255)
                                                                                          gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCAGAATTTATTTTTTCTCCTCTGTTTTTAGATACGAAAAAAGAAAATTAACCGAGGGAAAAAAAAATGGGATGATTTGCTAtt                                                       >  2:1091188/1‑151 (MQ=255)
                                                                                                aCTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAg                                                 <  2:212002/151‑1 (MQ=255)
                                                                                                aCTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGCACCATTTATTTTTTCTCACCCGGTTTTTGATGCGAAAAATGAAAATTAACAGGGGAAAAAAAAAAGGGAAAGATTTGGTAttttttag                                                 >  2:1614456/1‑147 (MQ=255)
                                                                                                   ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAGATGGGAAAAATGAAAATTAACGGAGGGAAAAAAAATGGGAATGATTTGCTATTTTTTAGaaa                                              >  1:1591229/1‑151 (MQ=255)
                                                                                                   ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGCACAATTTTTTTTTTCTCCTCTTTTTTTAGATGGGGAAAATGAAAATTAAAGGAGGAAAAATAAAAGGGGATGATTGGGTATTTTTCAGaaa                                              >  2:890734/1‑151 (MQ=255)
                                                                                                   ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAACGACACGCACATTTTTTTTTTTCCCCTCCTTTTTTAAAAGAGGAAAACAAAATTCAAAGGGGGAGAAAAAAAAGGGGAGGAATGGGGAttttttcgaaa                                              >  1:343549/1‑144 (MQ=255)
                                                                                                   ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACATTTTATTTTTTCTCCTCTGTTTTTAAATGGGGAAAAAGAAAATTAAAGGAGGAAAAAAAAAAGGGAATGATTGGCTATTTTTTAGaaa                                              >  2:246389/1‑151 (MQ=255)
                                                                                                     gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATGGGAAAAAATAAAATTAAAGGGGGGAAAAAAAAAGGGGATGATTTGGTTTTTTTTAGaaaat                                            >  1:33888/1‑149 (MQ=255)
                                                                                                     gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAATGTACATTTTTTTTTTCCTCCTCTTTTTTAAAATGCGAAAAATGAAATTCAACGGGGGGAAAAAAAATGGGGAAGATTGGCTATTTTGTAGaaaat                                            >  1:933768/1‑149 (MQ=255)
                                                                                                     gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGGACATTTTTTTTTTCCTCATCTGTTTTTAAATGGGGAACAAGAAAATCAAAGGGGGGAAAATAAATGGGGATGATTGGGTATTTTTCAGaaaat                                            >  1:388207/1‑149 (MQ=255)
                                                                                                        ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGACAAGAACATTTTATTTTTTCTCCTCTGTTTTTAGATGCGAAAAATGAAAATTAACGGAGGGAAAATAAATGGGAAAGATTGGGTATTTTTTAGAAAATGag                                         >  2:38684/1‑150 (MQ=255)
                                                                                                          acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc                                       <  1:993846/151‑1 (MQ=255)
                                                                                                           caacaGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGGAAAAACGAACAAAACATTTTTTTTTTCCTCATCTGTTTTTAAAATCGAAAAAAGAAAATAAAAGGGGGGAAAAAAAAAGGGAAAGATTGGGTGTTTTTCAaaaaattaggca                                      >  1:437021/1‑143 (MQ=255)
                                                                                                             acaGTAAATAATAAAAATTTCTCGAGGCGGGGGCCCAAAAAAAGACAAACACACTTTTTTTTTCCCCTTCTTTTTTAAAAAGCGGAACAAAAAATTCAAACGAGGGAAAAAAAATAGGGAAGATTGGCTATTTTTTAGAAAATGACTCAGa                                    >  1:980053/1‑151 (MQ=255)
                                                                                                             acaGCAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGGCAAGGACATTTTTTTTTTCCTCCTCTGTTTTTAGAAGGGAAAAATGAAAATTAAAGGGGGGAAAAAAAATGGGGAAGATTTGCTATTTTTCAGAAAATGACTCaaa                                    >  1:1593256/1‑149 (MQ=255)
                                                                                                                   aataatAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGAt                              <  2:370024/151‑1 (MQ=255)
                                                                                                                   aataatAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGAt                              <  1:88205/151‑1 (MQ=255)
                                                                                                                        taaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc                         <  2:951104/151‑1 (MQ=255)
                                                                                                                         aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa                        <  1:23191/151‑1 (MQ=255)
                                                                                                                             aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTg                    <  1:922608/151‑1 (MQ=255)
                                                                                                                                       gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc            <  2:467245/149‑1 (MQ=255)
                                                                                                                                       gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc            <  2:254761/149‑1 (MQ=255)
                                                                                                                                       gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc            <  2:1583174/149‑1 (MQ=255)
                                                                                                                                           ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc                                       >  1:641913/1‑119 (MQ=255)
                                                                                                                                           ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc                                       <  2:641622/119‑1 (MQ=255)
                                                                                                                                            gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGa      >  2:1203580/1‑151 (MQ=255)
                                                                                                                                            gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGa      >  2:1043894/1‑151 (MQ=255)
                                                                                                                                             tCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGa      <  2:1056194/150‑1 (MQ=255)
                                                                                                                                                  gAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaaa  >  1:907293/1‑150 (MQ=255)
                                                                                                                                                   aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaaa  <  2:291972/149‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                           
ACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAA  >  NC_000913/2386712‑2387005

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