Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,536,1970AT50.0% 46.8 / 88.2 74intergenic (+166/+417)yjhR/nanSpseudogene, helicase family; putative frameshift suppressor/9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base A (6/31);  new base T (37/0);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 6.98e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536054‑4536346
                                                                                                                                                 |                                                                                                                                                       
cAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGG‑CTTTTTggga                                                                                                                                                     >  2:244481/1‑146 (MQ=255)
  ggAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATgaaa                                                                                                                                                <  2:413892/151‑1 (MQ=255)
          aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑GTTTTTGGAATGaaaaaaaa                                                                                                                                           >  2:1017128/1‑148 (MQ=255)
          aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATgc                                                                                                                                        <  2:708892/151‑1 (MQ=255)
            aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATgcaa                                                                                                                                      <  1:1023498/151‑1 (MQ=255)
               gCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATgcaagca                                                                                                                                   <  1:524617/151‑1 (MQ=255)
                 atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCaat                                                                                                                                 <  2:684146/151‑1 (MQ=255)
                  tatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑GTTTTTTGAATTAAAAAAAAAaaaagcaa                                                                                                                                  >  1:345692/1‑149 (MQ=255)
                      ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGG‑GTTTTTGGGATGAAAAAAAAAGAAAGCaataaaa                                                                                                                             >  1:898202/1‑148 (MQ=255)
                       tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTGGGATGAAAAAAAAAGAAAACaaaaaaa                                                                                                                             >  2:722310/1‑147 (MQ=255)
                       tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAataatat                                                                                                                            <  2:686305/150‑1 (MQ=255)
                          cAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑GTTTTTTGGATGAAAAAAAAAACAAACAaaaaaat‑aa                                                                                                                         >  2:701737/1‑149 (MQ=255)
                             tGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑GCTTTTTGAATTAAAAAAAAAGCAAGCAATAAAAa‑aataa                                                                                                                      >  2:1126495/1‑149 (MQ=255)
                                tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGc                                                                                                                  <  2:430570/150‑1 (MQ=255)
                                tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGG‑GTTTTTGGGATGAAAAAAAAAAAAAGAAAAAAAAt‑aataa                                                                                                                      >  1:952508/1‑146 (MQ=255)
                                 aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑GTTTTTTGGATGAAAAAAAAAAGAAAGAATAATAT‑AAAAAAgggc                                                                                                                 >  1:125485/1‑148 (MQ=255)
                                 aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCt                                                                                                                <  1:520334/151‑1 (MQ=255)
                                     cGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGG‑CTTTTTTGAATGAAAGAAAAAGCAAGCAATAATAT‑AAAAAAAGGCTAAtt                                                                                                            >  2:578447/1‑151 (MQ=255)
                                      gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAAttt                                                                                                           <  2:107618/151‑1 (MQ=255)
                                         tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTTTAATTTAAAAAAAATCAAGCAAAAATAT‑AAAAAATGCCTAAtttttt                                                                                                        >  1:324711/1‑151 (MQ=255)
                                          ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCt                                                                                                         <  1:315845/149‑1 (MQ=255)
                                           tCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCttt                                                                                                       <  1:900610/150‑1 (MQ=255)
                                             cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGaa                                                                                                    <  2:449993/151‑1 (MQ=255)
                                             cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGaa                                                                                                    <  1:467958/151‑1 (MQ=255)
                                             cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGG‑CTTTTTGGAATGAAAAAAAAAAAAAACAATAATAA‑AATAAATGCCTAAAtttttt                                                                                                       >  2:604525/1‑148 (MQ=255)
                                               ttttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGG‑CTTTTTTGGAAGAAAGAAAAAACAAGCAAAAAAAT‑AAAAAAAGACTAATttttttt                                                                                                      >  2:372298/1‑146 (MQ=255)
                                               ttttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGATGTTTTTTTTGAAT‑AAAAAAAAAGGAAACAAAAAAAT‑AATAATATTCTAATtttttttt                                                                                                     >  1:795492/1‑146 (MQ=255)
                                                       cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAg                                                                                          <  1:695017/151‑1 (MQ=255)
                                                             gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATGATAT‑AATAATGTCCTAATTTCTTTGAACATCATTAAGCAT‑atc                                                                                   <  2:1177986/151‑1 (MQ=255)
                                                               cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑GTTTTTGGGATGAAAGAAAAAAGAAACAAAAAAAT‑AAAAAATGGCTAATTTTTTTTTAAATCATTaaaaa                                                                                        >  2:597044/1‑142 (MQ=255)
                                                               cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTTGGATTAAAAAAAAAACAAGCAAAAAAAT‑AAAAAAAGCCTAATTTTTTTTTAAAACATTAAACa                                                                                        >  1:170899/1‑145 (MQ=255)
                                                                  tCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑GTTTTTGGAATGAAAAAAAAAGCAAGCAAAAAAAT‑AAAAAAAGGCTAATTTTTTTTTAAAACATTAAAAAA‑atcatca                                                                               >  2:553887/1‑150 (MQ=255)
                                                                   ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑atcatcatc                                                                             <  1:485116/151‑1 (MQ=255)
                                                                   ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑atcatcat                                                                              <  1:501113/150‑1 (MQ=255)
                                                                      gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTc                                                                          <  2:814892/151‑1 (MQ=255)
                                                                        ttATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTGGGATGAAAAAAAAAGCAAACAATAATAT‑AAAAAAAGCCTAATTTTTTTTAACAACATTAAACAA‑ATCATCATCGTCt                                                                         >  2:838860/1‑150 (MQ=255)
                                                                          ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGG‑CTTTTTGGAATTAAAAAAAAAGCAAACAATAATAT‑AAAAATGGGCTAATTTTTTTTTAAATCATTAAAAAA‑ATTATCAACGTCTAtt                                                                      >  2:1022032/1‑151 (MQ=255)
                                                                            aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTTGAATGAAAAAAAAAACAAAGAATAAAAT‑AATAAAAGGCTAATTTTTTTTTACAACATTTAGAAA‑ATCATCCTCGTCTAtttt                                                                    >  1:959321/1‑150 (MQ=255)
                                                                            aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTg                                                                    <  1:163208/151‑1 (MQ=255)
                                                                             aaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGG‑CTTTTTGGAACAAAAAAAAAATC‑AGCAACAATATAAAAAATGGCC‑AATTTCTTTTAAAATACTTAAAAAT‑attatcaaactccattt                                                                     >  1:1159499/1‑139 (MQ=255)
                                                                             aaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTTGTATGAAAAAAAAAAAAAGCAATAAAAT‑AAAAAATGGCTTATTTTTTTTTACATTATTAAACAT‑ATTATTATCCTCTtttttt                                                                   >  1:275674/1‑149 (MQ=255)
                                                                             aaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGG‑CTTTTTTGTATTAAAAAAAAAAAAAACAAAAAAAT‑AAAAAAAGGCTAAATTTTTTTTAAATCATTAAGCAA‑ATAATAATCATCTAttt                                                                     >  1:5907/1‑149 (MQ=255)
                                                                               tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTTGGAAGAAAAAAAAAGAAAGCAAAAAAAT‑AAAAAAAGCCTAATTTTTTTTAAAATCATTTAGCAT‑ATTATCATCGTCTAtttt                                                                    >  2:522936/1‑147 (MQ=255)
                                                                                 ttATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CGTTTTGTTGTGAAAAGAAATTCAAAGAATAAAAA‑AATAAAGGCATTATTTCTTTGAACATAATTAATAAT‑ATTATTATCATCTtttt                                                                     >  1:725808/1‑145 (MQ=255)
                                                                                   aTTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGG‑CTTTTTGGGATGAAAAAAAAAACAAACAAAAATAT‑AATAAAAGCCTTATTTTTTTTTAAATTATTAAAAAT‑ATCATCATCGTCTAtttttccgaa                                                              >  2:245941/1‑143 (MQ=255)
                                                                                     tGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑GTTTTTTGGATGAAAAAAAAAAAAAACAATAAAAT‑AATAAAAGGCTAATTTTTTTTTAAATCATTAAACAA‑ATTATTCTCTTCTTTTTTTCAGa                                                               >  1:1133771/1‑147 (MQ=255)
                                                                                     tGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTGGTATTAAAAAAAAAACAAACAATAATAT‑AATAAAAGGCTTTTTTTTTTTTAAAACATTAATAAT‑ATTATCATCGTTTAtttttcacaa                                                              >  1:425493/1‑141 (MQ=255)
                                                                                          aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑GTTTTTTGAATTAAAAAAAAAACAAGCAAAAAAAT‑AAAAAATGGCTTATTTTTTTTAACAACATTTAACAT‑AATATCATCGTCTATTTTTCAGAAACCGTc                                                        >  1:444851/1‑149 (MQ=255)
                                                                                           ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTTGAATTAAAAAAAAAAAAAACAAAAATAT‑AATAATAGTCTTATTTTTTTTTACATCATTAAACAT‑ATCATCATATTCTATTTTTCAGAATCCCTCGAc                                                     >  1:47137/1‑151 (MQ=255)
                                                                                           ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGAc                                                     <  1:614006/151‑1 (MQ=255)
                                                                                            cTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTGGAATGAAAAAAAATGCAAGCAATAATAT‑AATAATGGGCTAATTTTTTTTTACATCATTAAGCAT‑ATCATCATCGTCTATTTTCCACAATCCCTCGAc                                                     >  2:1121767/1‑150 (MQ=255)
                                                                                             tAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑TTTTTTTGGAAGAAAGAAAATATACACAAAAAAAT‑AATAAAAGCCACATTTTTTTGAACAAAAATAATAAT‑ACAATAATCGTTTATTTTCCCGCATcactcggact                                                   >  2:1143988/1‑142 (MQ=255)
                                                                                              aaTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTta                                                  <  1:317754/151‑1 (MQ=255)
                                                                                                    tttCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGG‑GTTTTTGGGATGAAAAAAAAAACAAGAAATAATAT‑AATAATAGCCTAATTTTTTTTAACATCATTTAAAAT‑ATTATCATCGTCTATTTTTCACAATCCGTCGACTTATaa                                               >  2:363300/1‑148 (MQ=255)
                                                                                                       ccTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCa                                         <  2:989081/151‑1 (MQ=255)
                                                                                                       ccTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCa                                         <  1:57945/151‑1 (MQ=255)
                                                                                                       ccTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCa                                         <  2:912210/151‑1 (MQ=255)
                                                                                                         tttGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTGGGATGAAAAAAAAAAAAAGCAATAAAAT‑AATAATAGCCTAATTTCTTTTTAAATCATTAAGAAA‑ATTATTATCATCTATTTTTCAGAATACCTCGACTTATAATAGcc                                          >  1:497704/1‑148 (MQ=255)
                                                                                                         tttGCATAATGCGCCTGGGATTAATGGC‑GGGGG‑CTTTTTGGGAAGAAAAAAAAAGCAAAAAATAAAAT‑AAAAAATGGCTAATTTTTTTTAAAAACATTAAAAAT‑ATCATCATCGTCTATTTTCCAGAAACCGTCGACTTATAATAGcc                                          >  2:767542/1‑148 (MQ=255)
                                                                                                           tGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGaa                                      <  1:932170/150‑1 (MQ=255)
                                                                                                             cATAATGCGCCTGGGATTAATGGC‑GGGGG‑CTTTTAGGGATGGAAAAAAAAGAAAACAATAAATT‑AAAAAAGGGCCAATTTTTTTTTACAACATTAAACATAATCA‑CATTGTCTTTTTTTCAGAAATCGTCGACTTTTaaattccacaaagg                                   >  1:300603/1‑139 (MQ=255)
                                                                                                              aTAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTTTGGATGAAAAAAAAAGCAAACAAAAATAT‑AAAAAAAGGCTAATTTTTTTTAAAATCATTAAACAT‑ATCATCATCGTCTATTTTTCAGAATCCGCCGACTTATAATTACCAGaa                                      >  2:340346/1‑147 (MQ=255)
                                                                                                               tAATGCGCCTGGGATTAATGGC‑GGGGG‑GTTTTTGGGATGAAAGAAAAAGAAAGGAAAAAAAT‑AAAAAAAGGCTAATTTTTTTTTAAATAATTAAGAAT‑ATAATCATCGTCTTTTTTTCAGAAAACGTCGGCTTATAAAATCCAGaa                                      >  1:827798/1‑146 (MQ=255)
                                                                                                                 aTGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCata                               <  1:527063/151‑1 (MQ=255)
                                                                                                                    cgcCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCa                            <  1:682344/151‑1 (MQ=255)
                                                                                                                         gggATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGa                         <  1:1059118/149‑1 (MQ=255)
                                                                                                                               aaTGGA‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgctg                 >  1:778982/1‑151 (MQ=255)
                                                                                                                                 tGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCg               >  1:1072072/1‑151 (MQ=255)
                                                                                                                                   gC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGaa             <  1:372293/151‑1 (MQ=255)
                                                                                                                                    c‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc            <  1:805919/151‑1 (MQ=255)
                                                                                                                                    c‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc            >  1:328410/1‑151 (MQ=255)
                                                                                                                                            cTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg      >  2:210516/1‑151 (MQ=255)
                                                                                                                                               ttAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGcc   >  1:195329/1‑151 (MQ=255)
                                                                                                                                                tAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGcc   <  1:30733/150‑1 (MQ=255)
                                                                                                                                                 aGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc  <  1:121553/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
CAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTG‑CTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCAT‑ATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536054‑4536346

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