Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9110GA46.2% 94.8 / 60.7 91intergenic (‑122/+199)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base G (46/3);  new base A (0/42);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 8.97e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001764‑2002059
                                                                                                                                                        |                                                                                                                                                     
ccTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaa                                                                                                                                                    <  1:985183/151‑1 (MQ=255)
  ttCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaa                                                                                                                                                   >  1:1025848/1‑150 (MQ=255)
       cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                             <  1:281105/150‑1 (MQ=255)
          tataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCAGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                         >  2:840460/1‑151 (MQ=255)
            taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgccg                                                                                                                                       >  1:260950/1‑151 (MQ=255)
            taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgccg                                                                                                                                       >  2:1024409/1‑151 (MQ=255)
            taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgccg                                                                                                                                       >  1:924979/1‑151 (MQ=255)
              gagTGAATTCACGATAAACAGCCCTGCGTTATAGGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCccgccg                                                                                                                                       <  2:433974/149‑1 (MQ=255)
               agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATACGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTg                                                                                                                                    >  1:1257283/1‑151 (MQ=255)
                     ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                              >  2:477262/1‑151 (MQ=255)
                       cacaaaaaccAGCCTGGCGTTATATGAGTTATGGGCATGATTTTCCGTTTCTGCGGGGTTTTTAATCGGACGTTTGGTGGGTGAAATGGGGGGTTTTTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGaa                                                                                                                            <  2:568281/142‑1 (MQ=255)
                        aCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                             >  1:513374/1‑149 (MQ=255)
                           aaaaaCGGCCCTGCGTTATATGATTTATCGGCATGATTATCCTTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTGAAATGGGGGGTTATTGGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                        <  1:602848/149‑1 (MQ=255)
                           aTAACCACCCCTGGTTTATATGAGTTATCGGTATAATTATCCGTTTCCGATGGGTTTTAAATCGGATATTTGGTGGGTGAAAAGAGGGGTTTTTGGGGGTTTCAAGGTAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAg                                                                                                                           <  1:291236/148‑1 (MQ=255)
                             aaaCAGCCCGGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGAGGATTAGGGGGTAAAAGGGGGGGTTATTGGGGGGTTAAAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACg                                                                                                                       <  2:378106/150‑1 (MQ=255)
                             aaaCACCCCTGCGTTATATGATTTATCGCATTTATTATCTGTTTCTGCAGGGTTTTAATTCGGACGATAAGTGGGTAAAAGGGGGGGTTATTGGGGGTTTAAAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCa                                                                                                                         <  2:850674/148‑1 (MQ=255)
                                 aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCTGGTGGCGGGGAAGCACGTTGCt                                                                                                                  >  2:1364336/1‑151 (MQ=255)
                                       gctttATATGATTTATCGAAATAATTATCCTTTTCGGCAGGTTTTTTAATCGGACGTTTAGTGGGTAAAAGGGGGGTTTATTGGGGGGTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                              <  1:1142034/146‑1 (MQ=255)
                                         gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                          >  1:596634/1‑151 (MQ=255)
                                          ttATATGAGTTATCGGCATGTTTTTCCTTTTCTGCAGGTTTTTTAATCGGAGGATTAGTGGGGAAAATGAGGGGTTATTTGGGGGTTACAGGAAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                            <  1:1113609/148‑1 (MQ=255)
                                           tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                        >  2:341337/1‑151 (MQ=255)
                                             tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                       >  1:1024081/1‑150 (MQ=255)
                                                gAGTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTAAAATGGGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                      <  1:884256/148‑1 (MQ=255)
                                                gAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                   >  2:134138/1‑151 (MQ=255)
                                                  gTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTGAAAGGGGGGTTTATTGGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                   <  2:375418/149‑1 (MQ=255)
                                                   ttttCGGAATGATTATCCTTTTCTGCGGGGTTTTTAAGCGAACAATTAGGGGGGAAAATGGGGGTTTATTGGGGGTTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                   <  1:1354048/145‑1 (MQ=255)
                                                   ttATCGGAATGATTATCCTTTTCGGCGGGGTTTTAAATCGGACGATTGGTGGGTAAAATGGGGGTTTATTGGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                  <  1:792538/149‑1 (MQ=255)
                                                    tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                               >  1:181863/1‑151 (MQ=255)
                                                      tCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                             >  1:1190337/1‑151 (MQ=255)
                                                         gCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                          >  1:855354/1‑151 (MQ=255)
                                                           aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                         >  1:1150227/1‑150 (MQ=255)
                                                              aTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                      >  2:991335/1‑150 (MQ=255)
                                                               ttttCCGTTTCTGCGGGTTTTTTAATCGGACGTTTGGTGGGTAAAAGGGGGGTTTTTTTGGGGGTTACAGTTAATTTCCAGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                     <  2:931114/147‑1 (MQ=255)
                                                               ttATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                    >  2:582703/1‑151 (MQ=255)
                                                                 aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                   >  2:1282828/1‑150 (MQ=255)
                                                                 aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                   >  2:814957/1‑150 (MQ=255)
                                                                   ccGTTTCTGCGGGTTTTTTAATCGGACGATTGGTGGGTGAAATGAGGGTTTATTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                 <  2:569233/150‑1 (MQ=255)
                                                                   ccGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                >  1:1278609/1‑151 (MQ=255)
                                                                      tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                              >  1:420484/1‑150 (MQ=255)
                                                                          tGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                         >  1:688273/1‑151 (MQ=255)
                                                                             aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                       >  2:295058/1‑150 (MQ=255)
                                                                               ggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                    >  2:700017/1‑151 (MQ=255)
                                                                                 tttttAAGGGGAAGATAAGGGGTTAAAAGGAGGTGTTATTGGGGGTTAACTGATAATTTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                       <  2:1063872/146‑1 (MQ=255)
                                                                                  ttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                 >  2:780373/1‑151 (MQ=255)
                                                                                   tttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                >  2:1379426/1‑151 (MQ=255)
                                                                                          ggaggataggTGGGTAAATTGAGGGTTTATGGGGGGGTTAAGGGTAAATCCAGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                           <  1:633409/140‑1 (MQ=25)
                                                                                          ggaggaTTGGGGGGTAAAATGGGGGGTTTTTGGGGGTTAAAAGTAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                          <  2:1403056/150‑1 (MQ=25)
                                                                                          ggACGATTGGGGGGTGAAATGGGGGTTTTTTGGGGGGTTACGGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                          <  2:1142738/150‑1 (MQ=37)
                                                                                          ggACGATTAGGGGGTGAAATGAGGGGTTTTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                          <  2:395826/150‑1 (MQ=255)
                                                                                           ggacgatggggggTGAAAATGGGGGTTTTTTGGGGGTTTAAAGTAAAATCCAGGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                           <  2:1180403/139‑1 (MQ=25)
                                                                                            aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                        >  2:1039451/1‑150 (MQ=255)
                                                                                            aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                        >  2:815371/1‑150 (MQ=255)
                                                                                              gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                     >  2:166958/1‑151 (MQ=255)
                                                                                               aTTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                     >  1:1207380/1‑150 (MQ=255)
                                                                                                ttAGTGGGGGAAAGGGGGGGTTTTTGGGGGTTTACAGGAAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                    <  1:638426/150‑1 (MQ=25)
                                                                                                ttAGGGGGTAAAAGGGGGGTTTTTTGGGGGGTTACAGATAATTCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                      <  2:39351/148‑1 (MQ=25)
                                                                                                  aGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                  >  1:783240/1‑150 (MQ=255)
                                                                                                   gggggTAAAATGGGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                   <  1:75998/146‑1 (MQ=37)
                                                                                                    gggggaaaaaggggggTTTTTTGGGGGTTAGCAGGAAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                     <  1:1282856/134‑1 (MQ=25)
                                                                                                    ggggaaaaaTGGGGGGTTATTGGGGGTTTACAGGAAAATTCGAGCGAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACAAAAAtgtg                                                   <  1:850950/141‑1 (MQ=37)
                                                                                                     gggaaaaaTGAGGGTTTATTGGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                               <  1:436352/145‑1 (MQ=38)
                                                                                                      ggTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                             >  2:1207073/1‑151 (MQ=255)
                                                                                                      ggTGAAAAGGGGGGTTTTTTGGGGGTTTAAAGTAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                  <  2:662362/147‑1 (MQ=25)
                                                                                                        tGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                            >  1:483701/1‑150 (MQ=255)
                                                                                                           aaaggggggGTATTTTGGGGTTAACAGTAAAATCCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                         <  2:138356/146‑1 (MQ=25)
                                                                                                              gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                     >  2:362132/1‑151 (MQ=255)
                                                                                                                gggtTTTTTTGGGGTTAACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                   <  2:924958/151‑1 (MQ=25)
                                                                                                                ggggTTATTTGGGGGTTACAGGTAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                    <  2:762434/150‑1 (MQ=255)
                                                                                                                 gggTTTTTTGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                  <  1:352210/151‑1 (MQ=37)
                                                                                                                  ggTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                  >  2:924465/1‑150 (MQ=255)
                                                                                                                    ttATTTGGGGTTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                                 <  1:796370/149‑1 (MQ=255)
                                                                                                                     tttttGGGGGTTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                                <  1:1238714/148‑1 (MQ=37)
                                                                                                                     tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                              >  2:954341/1‑151 (MQ=255)
                                                                                                                       tttGGGGTTTACATGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             <  1:1367581/150‑1 (MQ=255)
                                                                                                                       ttggggggTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             <  1:367280/147‑1 (MQ=255)
                                                                                                                       ttggggggTTAAAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             <  2:244217/147‑1 (MQ=37)
                                                                                                                              gTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                      >  1:370268/1‑150 (MQ=255)
                                                                                                                                 tACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                   >  1:1117485/1‑151 (MQ=255)
                                                                                                                                 tACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                   >  2:143227/1‑151 (MQ=255)
                                                                                                                                 tACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg                    <  1:134144/150‑1 (MQ=255)
                                                                                                                                            aattcccAGCCAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt             <  2:1025827/145‑1 (MQ=25)
                                                                                                                                            aaaTTCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt             <  1:1176619/150‑1 (MQ=255)
                                                                                                                                            aaaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt             <  1:1132009/150‑1 (MQ=255)
                                                                                                                                             aaatCCAGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg                    <  2:260943/140‑1 (MQ=255)
                                                                                                                                             aaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa            >  2:1395926/1‑150 (MQ=255)
                                                                                                                                             aaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa           <  1:991936/151‑1 (MQ=255)
                                                                                                                                               ttCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagc         >  1:978365/1‑151 (MQ=255)
                                                                                                                                                 caaGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg        <  1:1395011/148‑1 (MQ=255)
                                                                                                                                                  cAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc       >  1:517350/1‑150 (MQ=255)
                                                                                                                                                      cAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc  >  2:393621/1‑151 (MQ=255)
                                                                                                                                                      cAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc  >  1:66690/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001764‑2002059

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