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

ATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536051‑4536346
                                                                                                                                                   |                                                                                                                                                        
aTACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGaa                                                                                                                                                      <  3:822133/151‑1 (MQ=255)
  aCAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAAt                                                                                                                                                     <  2:427642/150‑1 (MQ=255)
  aCAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATg                                                                                                                                                    <  3:844639/151‑1 (MQ=255)
   cAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTggt                                                                                                                                                       >  3:863894/1‑146 (MQ=255)
   cAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATga                                                                                                                                                   <  2:82917/151‑1 (MQ=255)
    aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATgaa                                                                                                                                                  <  2:283074/151‑1 (MQ=255)
      gAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGGTTTTGGGATgaaa                                                                                                                                                 >  3:1017207/1‑150 (MQ=255)
         tagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGCTTTTTTGTATgaaa                                                                                                                                                 >  3:875528/1‑147 (MQ=255)
          agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATgaaagaaa                                                                                                                                             <  3:86018/150‑1 (MQ=255)
           gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTTGTATGAAAGaaaa                                                                                                                                            >  3:170345/1‑150 (MQ=255)
           gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAAt                                                                                                                                           <  2:171336/151‑1 (MQ=255)
            taAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGCTTTTTTGGATGaaaaaaaa                                                                                                                                            >  2:734500/1‑149 (MQ=255)
             aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATgc                                                                                                                                         <  2:1226092/151‑1 (MQ=255)
              aaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATgca                                                                                                                                        <  3:687185/151‑1 (MQ=255)
              aaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATgca                                                                                                                                        <  2:558771/151‑1 (MQ=255)
               aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATgcaa                                                                                                                                       <  3:674035/151‑1 (MQ=255)
                 ggCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATgcaagc                                                                                                                                     <  2:1528282/151‑1 (MQ=255)
                  gCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGTCATTTAGGGAAGaaaaaaaaaataa                                                                                                                                       >  3:1329260/1‑143 (MQ=255)
                   cATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTGGAATGAAAGaaaaaaaaaaca                                                                                                                                    >  2:90594/1‑142 (MQ=255)
                    atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTTGGAAGAAAAAAAATgaaaaca                                                                                                                                    >  3:1618503/1‑149 (MQ=255)
                      atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGCTTTTTTGTATGaaaaaaaaagaaaaaaaaa                                                                                                                                 >  3:1522125/1‑139 (MQ=255)
                       ttttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTTGAAAGAAAAAAAAAAAAAACaataa                                                                                                                                >  2:258668/1‑150 (MQ=255)
                         ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataata                                                                                                                              <  2:735210/150‑1 (MQ=255)
                         ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTAGTATGGAaaaaaaaaaaaacaaa                                                                                                                                  >  3:1633390/1‑136 (MQ=255)
                           aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGCGTTTTTGGATGAAAAAAAAAAAAAGCaaaaaaa                                                                                                                              >  3:4545/1‑146 (MQ=255)
                            ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatataat                                                                                                                          <  3:874991/151‑1 (MQ=255)
                               aTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATaat‑aat                                                                                                                      <  2:1292438/151‑1 (MQ=255)
                               aTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATaat‑aat                                                                                                                      <  2:973669/151‑1 (MQ=255)
                                tGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTTGAATGAAAAAAAAAGAAAAAAAAAAAAtaat‑aa                                                                                                                       >  3:1398371/1‑149 (MQ=255)
                                  ttAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGCAGCTTTTTTGAATTAAAAAAAAAAGAAACaataaaataaaa                                                                                                                         >  2:820439/1‑144 (MQ=255)
                                   tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGc                                                                                                                   <  2:833676/150‑1 (MQ=255)
                                   tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGc                                                                                                                   <  2:1123259/150‑1 (MQ=255)
                                         gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTTCCTAATCATTTTCCTTTGCATAATGCGCCTAGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGTCCT‑AAtt                                                                                                            <  3:494004/150‑1 (MQ=255)
                                         gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTTGGATGAAAAAAAAAAAAAGCAATAATATAAT‑AAAAGACT‑AAttt                                                                                                           >  2:257273/1‑151 (MQ=255)
                                         gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AAttt                                                                                                           <  2:1243808/151‑1 (MQ=255)
                                             ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCt                                                                                                         <  2:183631/149‑1 (MQ=255)
                                               ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGa                                                                                                     <  2:956864/151‑1 (MQ=255)
                                               ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGa                                                                                                     <  2:930884/151‑1 (MQ=255)
                                                 gTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAAc                                                                                                   <  3:1552099/151‑1 (MQ=255)
                                                    ttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTTTTATGAAAAAAAATGCAAGCAATAATATAAA‑AAAGGGCT‑AATttttttttaaa                                                                                                  >  3:650028/1‑144 (MQ=255)
                                                    ttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAAcat                                                                                                 <  2:168299/150‑1 (MQ=255)
                                                      aTTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAAcatca                                                                                               <  3:113237/150‑1 (MQ=255)
                                                      aTTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAAcatca                                                                                               <  2:1640202/150‑1 (MQ=255)
                                                         gCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTaa                                                                                           <  3:745378/151‑1 (MQ=255)
                                                          cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAg                                                                                          <  3:914758/151‑1 (MQ=255)
                                                           gAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTTGGATGAAAAAAAAAGCAAGAAAAAATATAAA‑AAAAAGCT‑AATTTTTTTTTACATCAtt                                                                                             >  2:1335767/1‑147 (MQ=255)
                                                             gTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTTGGATTAAAAAAAAAGCAAGCCAAAAAATAAT‑AAAGGGCT‑AAATTTTTTTTACATAATTaaaaa                                                                                        >  3:827084/1‑147 (MQ=255)
                                                               cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTGGGATGAAAAAAAAAGCAAGAAATAATATAAA‑AAATGCCT‑AATTTTTTTTAAAATCATTaaaaa                                                                                        >  2:1459101/1‑145 (MQ=255)
                                                               cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCatat                                                                                     <  3:1627474/151‑1 (MQ=255)
                                                                gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTTGGATGAAAAAAAAAAAAAACAAAAAAATAAT‑AATGGGCT‑AATTTTTTTTAAAATCATTAAACata                                                                                      >  3:739499/1‑149 (MQ=255)
                                                                gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTTGGATGAAAGAAAAAACAAGCAAAAATATAAA‑AAAAGGCT‑AATTTTTTTTTAAATCAtttaaca                                                                                        >  2:41421/1‑142 (MQ=255)
                                                                 ccTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATatca                                                                                   <  3:86100/151‑1 (MQ=255)
                                                                    ttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTTGGATGAAAAAAAAAAGAAACAAAAAAAAAAA‑AAAAGACT‑TATTTTTTTTTACATCATTAAACata                                                                                      >  3:524932/1‑145 (MQ=255)
                                                                      ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATatcatcatc                                                                              <  2:1070394/151‑1 (MQ=255)
                                                                      ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATatcatcatc                                                                              <  3:1381683/151‑1 (MQ=255)
                                                                      ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATatcatcat                                                                               <  3:1334436/150‑1 (MQ=255)
                                                                         gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCg                                                                             <  2:730437/149‑1 (MQ=255)
                                                                         gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCg                                                                             <  2:18332/149‑1 (MQ=255)
                                                                         gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTc                                                                           <  2:993860/151‑1 (MQ=255)
                                                                          aTTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCt                                                                          <  2:415714/151‑1 (MQ=255)
                                                                          aTTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAAAGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCt                                                                          >  2:9441/1‑151 (MQ=255)
                                                                           ttATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTTGAATGAAAAAAAAAGCAAACAAAAATATAAT‑AAAATCCT‑CATTTTTTTTTAAATCATTTAACATATTATCATCGTCt                                                                          >  2:1326132/1‑150 (MQ=255)
                                                                            tataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTTGTATGAAAAAAAAAGCAAGCAAAAATATAAT‑AAATGCCT‑AATTTTTTTTAAAATCATTAAACATatcatcatcct                                                                            >  3:552692/1‑145 (MQ=255)
                                                                             ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTTTTATGAAAGAAAATGAAAGCAATAATATAAT‑AATGGTCT‑CATTTTTTTTAACATCATTTAACATATCATCATTGTCT‑At                                                                       >  3:1634161/1‑150 (MQ=255)
                                                                               aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTTGAAAGAAAGAAAAAACAAGCAATAATATAAT‑AATAGCCT‑TATTTTTTTTTACATCATTTAACAAATAATTATTTTCT‑Atttt                                                                    >  2:835965/1‑150 (MQ=255)
                                                                               aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTGGAATGAAAAAAAATGCAAGCAAAAATATAAT‑AATGGCCT‑AATTTTTTTTAACATCATTAAGCATATCATCATCGTCT‑Attt                                                                     >  3:808392/1‑150 (MQ=255)
                                                                               aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTg                                                                    <  3:579200/151‑1 (MQ=255)
                                                                                aagcggaggggggaCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGc                                                                   <  3:1047028/138‑1 (MQ=255)
                                                                                  tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTTGAATGAAAAAAAAAGCAAACAATAATATAAT‑AAAAGGCT‑AATTTTTTTGAAAAAAATTAAAAATATTATCATCGTCT‑ttttttc                                                                  >  2:1211401/1‑147 (MQ=255)
                                                                                    ttATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTTCTTTTTTGAATGAAAGAAAAAGCAAACAAAAATATTAT‑AAAAGCCT‑CATTTTTTTTTAAATCATTAAACATATTATCATCTTTT‑Attt                                                                     >  3:73068/1‑145 (MQ=255)
                                                                                           ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGt                                                         <  3:25002/150‑1 (MQ=255)
                                                                                                tAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTTGGATGAAAAAAAAAACAAGCCATAATATAAA‑AAAAGGCT‑AATTTTTTTTTAAAACATTTAACAAATCATCATCCTCT‑ATTTTCCCGAAATCGTCGACtt                                                   >  3:519124/1‑151 (MQ=255)
                                                                                                  aTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTtat                                                 <  2:982268/151‑1 (MQ=255)
                                                                                                  aTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTtat                                                 <  3:1374187/151‑1 (MQ=255)
                                                                                                     aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTTGGATGAAAAAAAAAACAAACAAAAAAATAAA‑AATAGGCT‑AATTTTTTTTAACATCAATAAACATATCATCATCGTCT‑ATTTTCCCCAAACCCTCGACTTATaa                                               >  2:744986/1‑150 (MQ=255)
                                                                                                     aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTGGGATGAAAAAAAAAGCAAGCAATAAAATAAA‑AAAAGCCT‑AATTTTTTTTTAAAACATTTAGCAAATAATCATCGTCT‑ATTTTTCAGAATCCCTCGACTTTTaa                                               >  2:1085621/1‑150 (MQ=255)
                                                                                                      ttttCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGCTTTTTTGAATTAAAAAAAAAACAAACAATAAAAAAAA‑AATAGCCT‑TAATTTTTTTAAAAACATTAAACATAAAAATATCCTTT‑ATTTACCTAAATACGTCGTCTTTTaa                                               >  2:1065992/1‑149 (MQ=255)
                                                                                                         tCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTTGGATGAAAAAAAAAACAAGCAAAAAAATAATAAATGGCC‑‑AATTTTTTTTAAAAACATTAAGAATATCATCATCGTCT‑ATTTTCCAGAAACCGTCGACCTTTaaatccc                                          >  2:1499253/1‑146 (MQ=255)
                                                                                                         tCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTTGGATTAAAAAAAAAGCAAGCAATAATATAAA‑AAATGCCT‑AATTTTTTTTAAAAACATTAAGCAAATCCTCCTCCTCTCATTTG‑CAGAATCCCTCCGCTTATAAttacc                                          >  2:1278909/1‑147 (MQ=255)
                                                                                                              tGCATAATGCGCCTGGGATTAATGGC‑GGGGGCTTTTTGGGATGAAAGAAAAAGCAAGCAAAAATATAAA‑AAAAGGCT‑AATTTTTTTTAAAATCATTTAACATATACTCATCGTCT‑ATTTTTCCGAAACCGTCGACTTATAATCGCCAGaa                                      >  2:12483/1‑150 (MQ=255)
                                                                                                               gCATAATGCGCCTGGGATTAATGGC‑GGGGGCGTTTTGGAATGAAAAAAAAAGCAAGGAAAAAAATAAA‑AAATGTCT‑AAATTCTTTTAAAAACAATAAACATATAAATAAAGTCT‑TTTTTTCAGAATCCGTCGACTTATAATAGCCAGaa                                      >  3:563617/1‑149 (MQ=255)
                                                                                                                cATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACgg                                   <  3:1543561/151‑1 (MQ=255)
                                                                                                                cATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACgg                                   <  2:1087370/151‑1 (MQ=255)
                                                                                                                 aTAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTGGGAAGAAAGAAAAAACAAACAAAAATATAAT‑AATAGCCTCAATTTTTTTTAAAACAATTAA‑AAAATCATCTTTTTCT‑TTTTTTCCGAATACGTCGACTTTTaaatgccccaaaag                                   >  2:1371326/1‑138 (MQ=255)
                                                                                                                 aTAATGCGCCTGGGATTAATGGC‑GGGTGGCTTTTGGAATGAAAAAAAAAGCAAGCAATAAAATAAT‑AATAGCCT‑AATTTTTTTTAAAAACATTAAACATATCATCATCCTCT‑ATTTTTCAGAATCCGTCGACTTATAATCGCCAGAAAgg                                   >  2:769181/1‑150 (MQ=255)
                                                                                                                 aTAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGc                                  <  2:633044/151‑1 (MQ=255)
                                                                                                                 aTAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGc                                  <  3:108331/151‑1 (MQ=255)
                                                                                                                  tAATGCGCCTGGGATTAATGGCGGGGGGCTTTTTGG‑ATGAAAAAAAAAACAAACAAAAAAATAAA‑AAAAGCCT‑AATTTTTTTTAAAATCATTAAACATATTATCATCGTCT‑ATTTTCCAGAAAACCTCGACTTATAATTGCCAGAACgg                                   >  2:1302472/1‑149 (MQ=255)
                                                                                                                  tAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTTGGATGAAAAAAAAAACAAAAAAAAAAATAAT‑AATAGGCT‑AATTTTTTTTTAAAACAATAAACAAATTATCATTGTCT‑TTTTTCCCGAATACCTCGCCTTATAATATCCAGAACgg                                   >  3:1551212/1‑149 (MQ=255)
                                                                                                                   aaTGCGCCTGGGATTAATGGC‑GGGGGCTTTTAGGGATGAAAGAAAAAGCAAGCAATAAAATAAA‑AATAGCCT‑AATTTTTTTTAAAATCATTAAACATATCATCATCGTCT‑ATTCGCCAGAAACCGTCGCCTTTTAATCTCCACaa                                      >  2:1006843/1‑145 (MQ=255)
                                                                                                                           tGGGATTAATGGC‑GGGGGCTTTTTTGGATGAAAAAAAAAACAAGCAAAAATATAAT‑AATAGGCT‑AATTTTTTTTTACAACATTAAACATATCATCATCGTCT‑ATTTTTCAGAATCCGTCGACTTATAATAGCCAGAACGGCataaaaag                          >  3:1212755/1‑144 (MQ=255)
                                                                                                                              gATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGt                     <  2:1569347/151‑1 (MQ=255)
                                                                                                                              gATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGt                     <  2:9511/151‑1 (MQ=255)
                                                                                                                                      gC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGaa             <  3:90591/151‑1 (MQ=255)
                                                                                                                                       c‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc            <  3:630182/151‑1 (MQ=255)
                                                                                                                                         gggTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACg           >  2:673413/1‑151 (MQ=255)
                                                                                                                                              cTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg      >  2:1070257/1‑151 (MQ=255)
                                                                                                                                                   aGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc  <  3:455462/150‑1 (MQ=255)
                                                                                                                                                   aGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc  <  2:897592/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
ATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAAT‑AATGGCCT‑AATTTCTTTGAACATCATTAAGCATATCATCATCGTCT‑ATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  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: ‑