Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,536,1970AT42.3% 117.6 / 81.3 97intergenic (+166/+417)yjhR/nanSpseudogene, helicase family; putative frameshift suppressor/9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base A (5/51);  new base T (41/0);  total (46/51)
Fisher's exact test for biased strand distribution p-value = 3.39e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536051‑4536346
                                                                                                                                                   |                                                                                                                                                        
aTACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGaa                                                                                                                                                      <  3:822133/151‑1 (MQ=255)
  aCAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAAt                                                                                                                                                     <  2:427642/150‑1 (MQ=255)
  aCAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATg                                                                                                                                                    <  3:844639/151‑1 (MQ=255)
   cAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTggt                                                                                                                                                       >  3:863894/1‑146 (MQ=255)
   cAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATga                                                                                                                                                   <  2:82917/151‑1 (MQ=255)
    aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaa                                                                                                                                                  <  2:283074/151‑1 (MQ=255)
      gAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGGTTTTGGGATgaaa                                                                                                                                                 >  3:1017207/1‑150 (MQ=255)
         tagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGTATgaaa                                                                                                                                                 >  3:875528/1‑147 (MQ=255)
          agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaagaaa                                                                                                                                             <  3:86018/150‑1 (MQ=255)
           gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGTATGAAAGaaaa                                                                                                                                            >  3:170345/1‑150 (MQ=255)
           gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAAt                                                                                                                                           <  2:171336/151‑1 (MQ=255)
            taAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGaaaaaaaa                                                                                                                                            >  2:734500/1‑149 (MQ=255)
             aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgc                                                                                                                                         <  2:1226092/151‑1 (MQ=255)
              aaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgca                                                                                                                                        <  3:687185/151‑1 (MQ=255)
              aaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgca                                                                                                                                        <  2:558771/151‑1 (MQ=255)
               aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaa                                                                                                                                       <  3:674035/151‑1 (MQ=255)
                 ggCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaagc                                                                                                                                     <  2:1528282/151‑1 (MQ=255)
                  gCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGTCATTTAGGGAAGaaaaaaaaaataa                                                                                                                                       >  3:1329260/1‑143 (MQ=255)
                   cATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGAAAGaaaaaaaaaaca                                                                                                                                    >  2:90594/1‑142 (MQ=255)
                    atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGAAGAAAAAAAATgaaaaca                                                                                                                                    >  3:1618503/1‑149 (MQ=255)
                      atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGTATGaaaaaaaaagaaaaaaaaa                                                                                                                                 >  3:1522125/1‑139 (MQ=255)
                       ttttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAAAGAAAAAAAAAAAAAACaataa                                                                                                                                >  2:258668/1‑150 (MQ=255)
                         ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataata                                                                                                                              <  2:735210/150‑1 (MQ=255)
                         ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTAGTATGGAaaaaaaaaaaaacaaa                                                                                                                                  >  3:1633390/1‑136 (MQ=255)
                           aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCGTTTTTGGATGAAAAAAAAAAAAAGCaaaaaaa                                                                                                                              >  3:4545/1‑146 (MQ=255)
                            ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatataat                                                                                                                          <  3:874991/151‑1 (MQ=255)
                               aTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATaataat                                                                                                                      <  2:1292438/151‑1 (MQ=255)
                               aTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATaataat                                                                                                                      <  2:973669/151‑1 (MQ=255)
                                tGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAGAAAAAAAAAAAAtaataa                                                                                                                       >  3:1398371/1‑149 (MQ=255)
                                  ttAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGCAGCTTTTTTGAATTAAAAAAAAAAGAAACaataaaataaaa                                                                                                                         >  2:820439/1‑144 (MQ=255)
                                   tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGc                                                                                                                   <  2:833676/150‑1 (MQ=255)
                                   tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGc                                                                                                                   <  2:1123259/150‑1 (MQ=255)
                                         gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTTCCTAATCATTTTCCTTTGCATAATGCGCCTAGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGTCCTAAtt                                                                                                            <  3:494004/150‑1 (MQ=255)
                                         gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATGAAAAAAAAAAAAAGCAATAATATAATAAAAGACTAAttt                                                                                                           >  2:257273/1‑151 (MQ=255)
                                         gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAAttt                                                                                                           <  2:1243808/151‑1 (MQ=255)
                                             ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCt                                                                                                         <  2:183631/149‑1 (MQ=255)
                                               ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGa                                                                                                     <  2:956864/151‑1 (MQ=255)
                                               ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGa                                                                                                     <  2:930884/151‑1 (MQ=255)
                                                 gTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAc                                                                                                   <  3:1552099/151‑1 (MQ=255)
                                                    ttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTTATGAAAAAAAATGCAAGCAATAATATAAAAAAGGGCTAATttttttttaaa                                                                                                  >  3:650028/1‑144 (MQ=255)
                                                    ttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcat                                                                                                 <  2:168299/150‑1 (MQ=255)
                                                      aTTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcatca                                                                                               <  3:113237/150‑1 (MQ=255)
                                                      aTTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcatca                                                                                               <  2:1640202/150‑1 (MQ=255)
                                                         gCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTaa                                                                                           <  3:745378/151‑1 (MQ=255)
                                                          cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAg                                                                                          <  3:914758/151‑1 (MQ=255)
                                                           gAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAAGCAAGAAAAAATATAAAAAAAAGCTAATTTTTTTTTACATCAtt                                                                                             >  2:1335767/1‑147 (MQ=255)
                                                             gTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATTAAAAAAAAAGCAAGCCAAAAAATAATAAAGGGCTAAATTTTTTTTACATAATTaaaaa                                                                                        >  3:827084/1‑147 (MQ=255)
                                                               cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGGATGAAAAAAAAAGCAAGAAATAATATAAAAAATGCCTAATTTTTTTTAAAATCATTaaaaa                                                                                        >  2:1459101/1‑145 (MQ=255)
                                                               cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCatat                                                                                     <  3:1627474/151‑1 (MQ=255)
                                                                gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAAAAAAACAAAAAAATAATAATGGGCTAATTTTTTTTAAAATCATTAAACata                                                                                      >  3:739499/1‑149 (MQ=255)
                                                                gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGAAAGAAAAAACAAGCAAAAATATAAAAAAAGGCTAATTTTTTTTTAAATCAtttaaca                                                                                        >  2:41421/1‑142 (MQ=255)
                                                                 ccTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatca                                                                                   <  3:86100/151‑1 (MQ=255)
                                                                    ttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGAAAAAAAAAAGAAACAAAAAAAAAAAAAAAGACTTATTTTTTTTTACATCATTAAACata                                                                                      >  3:524932/1‑145 (MQ=255)
                                                                      ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcatc                                                                              <  2:1070394/151‑1 (MQ=255)
                                                                      ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcatc                                                                              <  3:1381683/151‑1 (MQ=255)
                                                                      ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcat                                                                               <  3:1334436/150‑1 (MQ=255)
                                                                         gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCg                                                                             <  2:730437/149‑1 (MQ=255)
                                                                         gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCg                                                                             <  2:18332/149‑1 (MQ=255)
                                                                         gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTc                                                                           <  2:993860/151‑1 (MQ=255)
                                                                          aTTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCt                                                                          <  2:415714/151‑1 (MQ=255)
                                                                          aTTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAAAGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCt                                                                          >  2:9441/1‑151 (MQ=255)
                                                                           ttATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAAAAAAAGCAAACAAAAATATAATAAAATCCTCATTTTTTTTTAAATCATTTAACATATTATCATCGTCt                                                                          >  2:1326132/1‑150 (MQ=255)
                                                                            tataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGTATGAAAAAAAAAGCAAGCAAAAATATAATAAATGCCTAATTTTTTTTAAAATCATTAAACATatcatcatcct                                                                            >  3:552692/1‑145 (MQ=255)
                                                                             ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTTATGAAAGAAAATGAAAGCAATAATATAATAATGGTCTCATTTTTTTTAACATCATTTAACATATCATCATTGTCTAt                                                                       >  3:1634161/1‑150 (MQ=255)
                                                                               aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAAAGAAAGAAAAAACAAGCAATAATATAATAATAGCCTTATTTTTTTTTACATCATTTAACAAATAATTATTTTCTAtttt                                                                    >  2:835965/1‑150 (MQ=255)
                                                                               aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATGAAAAAAAATGCAAGCAAAAATATAATAATGGCCTAATTTTTTTTAACATCATTAAGCATATCATCATCGTCTAttt                                                                     >  3:808392/1‑150 (MQ=255)
                                                                               aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTg                                                                    <  3:579200/151‑1 (MQ=255)
                                                                                aagcggaggggggaCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGc                                                                   <  3:1047028/138‑1 (MQ=255)
                                                                                  tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAGCAAACAATAATATAATAAAAGGCTAATTTTTTTGAAAAAAATTAAAAATATTATCATCGTCTttttttc                                                                  >  2:1211401/1‑147 (MQ=255)
                                                                                    ttATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTTCTTTTTTGAATGAAAGAAAAAGCAAACAAAAATATTATAAAAGCCTCATTTTTTTTTAAATCATTAAACATATTATCATCTTTTAttt                                                                     >  3:73068/1‑145 (MQ=255)
                                                                                           ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGt                                                         <  3:25002/150‑1 (MQ=255)
                                                                                                tAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGAAAAAAAAAACAAGCCATAATATAAAAAAAGGCTAATTTTTTTTTAAAACATTTAACAAATCATCATCCTCTATTTTCCCGAAATCGTCGACtt                                                   >  3:519124/1‑151 (MQ=255)
                                                                                                  aTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTtat                                                 <  2:982268/151‑1 (MQ=255)
                                                                                                  aTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTtat                                                 <  3:1374187/151‑1 (MQ=255)
                                                                                                     aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATGAAAAAAAAAACAAACAAAAAAATAAAAATAGGCTAATTTTTTTTAACATCAATAAACATATCATCATCGTCTATTTTCCCCAAACCCTCGACTTATaa                                               >  2:744986/1‑150 (MQ=255)
                                                                                                     aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGGATGAAAAAAAAAGCAAGCAATAAAATAAAAAAAGCCTAATTTTTTTTTAAAACATTTAGCAAATAATCATCGTCTATTTTTCAGAATCCCTCGACTTTTaa                                               >  2:1085621/1‑150 (MQ=255)
                                                                                                      ttttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATTAAAAAAAAAACAAACAATAAAAAAAAAATAGCCTTAATTTTTTTAAAAACATTAAACATAAAAATATCCTTTATTTACCTAAATACGTCGTCTTTTaa                                               >  2:1065992/1‑149 (MQ=255)
                                                                                                         tCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAAACAAGCAAAAAAATAATAAATGGCCAATTTTTTTTAAAAACATTAAGAATATCATCATCGTCTATTTTCCAGAAACCGTCGACCTTTaaatccc                                          >  2:1499253/1‑146 (MQ=255)
                                                                                                         tCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATTAAAAAAAAAGCAAGCAATAATATAAAAAATGCCTAATTTTTTTTAAAAACATTAAGCAAATCCTCCTCCTCTCATTTGCAGAATCCCTCCGCTTATAAttacc                                          >  2:1278909/1‑147 (MQ=255)
                                                                                                              tGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAGAAAAAGCAAGCAAAAATATAAAAAAAGGCTAATTTTTTTTAAAATCATTTAACATATACTCATCGTCTATTTTTCCGAAACCGTCGACTTATAATCGCCAGaa                                      >  2:12483/1‑150 (MQ=255)
                                                                                                               gCATAATGCGCCTGGGATTAATGGCGGGGGCGTTTTGGAATGAAAAAAAAAGCAAGGAAAAAAATAAAAAATGTCTAAATTCTTTTAAAAACAATAAACATATAAATAAAGTCTTTTTTTCAGAATCCGTCGACTTATAATAGCCAGaa                                      >  3:563617/1‑149 (MQ=255)
                                                                                                                cATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACgg                                   <  3:1543561/151‑1 (MQ=255)
                                                                                                                cATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACgg                                   <  2:1087370/151‑1 (MQ=255)
                                                                                                                 aTAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGGAAGAAAGAAAAAACAAACAAAAATATAATAATAGCCTCAATTTTTTTTAAAACAATTAAAAAATCATCTTTTTCTTTTTTTCCGAATACGTCGACTTTTaaatgccccaaaag                                   >  2:1371326/1‑138 (MQ=255)
                                                                                                                 aTAATGCGCCTGGGATTAATGGCGGGTGGCTTTTGGAATGAAAAAAAAAGCAAGCAATAAAATAATAATAGCCTAATTTTTTTTAAAAACATTAAACATATCATCATCCTCTATTTTTCAGAATCCGTCGACTTATAATCGCCAGAAAgg                                   >  2:769181/1‑150 (MQ=255)
                                                                                                                 aTAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGc                                  <  2:633044/151‑1 (MQ=255)
                                                                                                                 aTAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGc                                  <  3:108331/151‑1 (MQ=255)
                                                                                                                  tAATGCGCCTGGGATTAATGGCGGGGGGCTTTTTGGATGAAAAAAAAAACAAACAAAAAAATAAAAAAAGCCTAATTTTTTTTAAAATCATTAAACATATTATCATCGTCTATTTTCCAGAAAACCTCGACTTATAATTGCCAGAACgg                                   >  2:1302472/1‑149 (MQ=255)
                                                                                                                  tAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGAAAAAAAAAACAAAAAAAAAAATAATAATAGGCTAATTTTTTTTTAAAACAATAAACAAATTATCATTGTCTTTTTTCCCGAATACCTCGCCTTATAATATCCAGAACgg                                   >  3:1551212/1‑149 (MQ=255)
                                                                                                                   aaTGCGCCTGGGATTAATGGCGGGGGCTTTTAGGGATGAAAGAAAAAGCAAGCAATAAAATAAAAATAGCCTAATTTTTTTTAAAATCATTAAACATATCATCATCGTCTATTCGCCAGAAACCGTCGCCTTTTAATCTCCACaa                                      >  2:1006843/1‑145 (MQ=255)
                                                                                                                           tGGGATTAATGGCGGGGGCTTTTTTGGATGAAAAAAAAAACAAGCAAAAATATAATAATAGGCTAATTTTTTTTTACAACATTAAACATATCATCATCGTCTATTTTTCAGAATCCGTCGACTTATAATAGCCAGAACGGCataaaaag                          >  3:1212755/1‑144 (MQ=255)
                                                                                                                              gATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGt                     <  2:1569347/151‑1 (MQ=255)
                                                                                                                              gATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGt                     <  2:9511/151‑1 (MQ=255)
                                                                                                                                      gCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGaa             <  3:90591/151‑1 (MQ=255)
                                                                                                                                       cGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc            <  3:630182/151‑1 (MQ=255)
                                                                                                                                         gggTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACg           >  2:673413/1‑151 (MQ=255)
                                                                                                                                              cTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg      >  2:1070257/1‑151 (MQ=255)
                                                                                                                                                   aGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc  <  3:455462/150‑1 (MQ=255)
                                                                                                                                                   aGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc  <  2:897592/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
ATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536051‑4536346

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