Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,536,1970AT46.7% 67.1 / 68.0 75intergenic (+166/+417)yjhR/nanSpseudogene, helicase family; putative frameshift suppressor/9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base A (4/36);  new base T (35/0);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 3.63e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536047‑4536346
                                                                                                                                                        |                                                                                                                                                      
cATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTa                                                                                                                                                        <  1:367532/151‑1 (MQ=255)
   gATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGa                                                                                                                                                     <  1:305880/151‑1 (MQ=255)
    aTACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGaa                                                                                                                                                    <  1:619202/151‑1 (MQ=255)
    aTACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTggg                                                                                                                                                     >  2:490564/1‑149 (MQ=255)
       cAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATga                                                                                                                                                 <  2:90449/151‑1 (MQ=255)
        aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATgaa                                                                                                                                                <  2:346392/151‑1 (MQ=255)
            gtagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCGTTTTTGAAt                                                                                                                                                   >  1:288142/1‑144 (MQ=255)
               gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAAt                                                                                                                                         <  2:416140/151‑1 (MQ=255)
                   aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATgcaa                                                                                                                                     <  1:211957/151‑1 (MQ=255)
                   aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATgcaa                                                                                                                                     <  2:873546/151‑1 (MQ=255)
                       cATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGGTTTTGGAATTAAAGAAAAAgcaa                                                                                                                                     >  2:27028/1‑147 (MQ=255)
                          atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTTGGAATAaaaaaaaaacaaacaaaa                                                                                                                                >  1:713890/1‑140 (MQ=255)
                          atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTTGGGATGAAAGAAAAAgcaaacaaaa                                                                                                                               >  1:144543/1‑148 (MQ=255)
                          atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCaata                                                                                                                               <  2:293458/150‑1 (MQ=255)
                            tttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GGTTTTTTGGATTTAAGAAAAAGCCAGCaataatt                                                                                                                            >  2:543513/1‑150 (MQ=255)
                             ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAataata                                                                                                                            <  1:481923/150‑1 (MQ=255)
                                   aTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑aataat                                                                                                                    <  2:484385/151‑1 (MQ=255)
                                     gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATgg                                                                                                                  <  2:518112/151‑1 (MQ=255)
                                     gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATgg                                                                                                                  <  1:667825/151‑1 (MQ=255)
                                        aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCt                                                                                                               <  2:978205/151‑1 (MQ=255)
                                        aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTTGGATGAAAAAAAATGGAAAAAAAAATAt‑aaaaa                                                                                                                     >  2:120271/1‑145 (MQ=255)
                                           gCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTTGAATTAAAAAAAAAACAAGCAATAATAt‑aataaatggctta                                                                                                             >  1:404013/1‑142 (MQ=255)
                                             gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GGTTTTTTGGATGAAAGAAAAAGGAAACAATAAAAT‑AATAAATGGCTAAttt                                                                                                          >  2:591765/1‑151 (MQ=255)
                                             gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAAttt                                                                                                          <  1:427697/151‑1 (MQ=255)
                                                tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGCTTTTTGTATTGAAAAAAATAAAAAGCAAT‑ATAT‑AATAAAggggcaatttttt                                                                                                       >  1:503420/1‑140 (MQ=255)
                                                tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTTGAATGAAAAAAAAAGGAAGCAATAATAT‑AAAAAATGCCTAAtttttt                                                                                                       >  2:17629/1‑151 (MQ=255)
                                                   ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGC‑ATAATAT‑AATAATGGCCTAATTTCTTTGa                                                                                                    <  2:991103/150‑1 (MQ=255)
                                                   ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAAACTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGa                                                                                                    <  2:288271/151‑1 (MQ=255)
                                                    cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGaa                                                                                                   <  1:473590/151‑1 (MQ=255)
                                                     gTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTGGGATGAAAAAAAAAGCAAAGAATAATAT‑AATAATAGGCTAATTTCttttta                                                                                                   >  2:882580/1‑147 (MQ=255)
                                                       tttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTGGGATGAAAAAAAAAAAAAGCAAAAATAT‑AAAAAATGCCCAATttttttttaca                                                                                                 >  2:708148/1‑145 (MQ=255)
                                                           ttGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTGGGATGAAAGAAAAAACAAGGAAAAATAT‑AATAAAAGCCCAATttttttttaca                                                                                                 >  2:76843/1‑141 (MQ=255)
                                                            tGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCAtt                                                                                            <  2:208328/150‑1 (MQ=255)
                                                             gCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTa                                                                                           <  2:220362/150‑1 (MQ=255)
                                                              cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GTGTTTTTGGATGAAAAAAAAAGCAAGCAAAAAAATAAAAAAAGCCCTAATtttttttacatta                                                                                               >  2:840281/1‑139 (MQ=255)
                                                                aGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCa                                                                                       <  2:239356/151‑1 (MQ=255)
                                                                aGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCa                                                                                       <  2:582503/151‑1 (MQ=255)
                                                                 gTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCat                                                                                      <  1:58756/151‑1 (MQ=255)
                                                                 gTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTGGGATGAAAGAAAAAGCAAAAAATAATAT‑AATAAAAGCCTAATTTTTTTTTAAATCATTAAACa                                                                                       >  2:941589/1‑150 (MQ=255)
                                                                   cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCata                                                                                     <  1:814121/150‑1 (MQ=255)
                                                                     ccTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATatca                                                                                  <  2:377368/151‑1 (MQ=255)
                                                                     ccTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTTGAAAGAAAAAAAAAGCAAACAATAATAT‑AATAAATGCCTTATTTTTTTTTAAATCATTAAACa                                                                                       >  2:881608/1‑146 (MQ=255)
                                                                      cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTAGGGATGAAAAAAAAAACAAACAAAAAAAT‑TATAATTGGTTTATTTCTTTTTACATCCTTa                                                                                           >  2:338100/1‑141 (MQ=255)
                                                                        ttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcat                                                                                 <  1:162475/149‑1 (MQ=255)
                                                                          ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTTGAAAGAAAAAAAATGAAAGAAATAATAT‑AATAAATGGCTAATTTTTTTTTACATCATTaaaaaaat                                                                                    >  2:501764/1‑138 (MQ=255)
                                                                          ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTGGGATTAAAAAAAAAGCAAGGAAAAATAA‑AATAATAGGCTAATTTTTTTTTAAATCCTTAAACaaattaaca                                                                               >  2:951841/1‑144 (MQ=255)
                                                                                 ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTGGAAATAAAAAAAAAAAAAGCAATAATAA‑AAAAAAAGGCTAAATTCTTTTTACATTCTTAAACATAAAATCATCCTCtttt                                                                      >  2:656840/1‑151 (MQ=255)
                                                                                   aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTTGAATGAAAAAAAAAGAAAGAAAAAAAAT‑AAAAATTGACCAAATTTTTTTAAAATCATTAAACATATCATCATCGTCTAttt                                                                     >  2:175542/1‑150 (MQ=255)
                                                                                     aTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTTGAATGAAAAAAAAAGCAAAACAAAAAAT‑AATAAAGGGCTAATTTTTTTTTAAATCATTTAAAATATTATTATCCTCTCttt                                                                     >  2:1018880/1‑148 (MQ=255)
                                                                                         tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGa                                                               <  1:45767/150‑1 (MQ=255)
                                                                                         tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGa                                                               <  2:604851/150‑1 (MQ=255)
                                                                                            tGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTTTGGATGAAAGAAAAAGCAAGCAATAAAAT‑AAAAAAGGCCCAATTTTTTTTTAAATCAATTAACATATCATCATCGTCTATTTTTCCGAATcc                                                           >  2:868329/1‑151 (MQ=255)
                                                                                               ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGt                                                         <  2:796832/150‑1 (MQ=255)
                                                                                                tACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTTGATTGTAAGAAAAAGCAAAGAAAAAAAT‑AATAAAGGGGTAATTTTTTTTTAAAACATTTAAGATATCATTTATGTCTATTTGCCaaaaattgtcg                                                       >  1:329270/1‑144 (MQ=255)
                                                                                                   cTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTGGGATGAAAAAAAAAGCAAACAAAAAAAT‑AAAAAAAGGCTAATTTTTTTTTAAATCATTTAAAAAATAATCATCGTCTATTTTTCAGAAACCg                                                          >  2:672566/1‑145 (MQ=255)
                                                                                                        cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATaa                                               <  1:672607/151‑1 (MQ=255)
                                                                                                        cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTTGAATGAAAAAAAAAGAAAGCAAAAATAT‑AATAATAGGCTAATTTTTTTTTAAATCATTAAGCATAATATCATCGTCTTTTTTTCAGAATCCGTCGACTtat                                                 >  2:550303/1‑149 (MQ=255)
                                                                                                             tCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTTGGGATGAAAAAAAAACCAAGCAATAAAAT‑AATAAAAGCCCAATTTTTTTTAATATTAATAAACATATCATCCTTGTCTATTTTTCAAAATCCGTCGACCTTTaa                                               >  2:628111/1‑146 (MQ=255)
                                                                                                             tCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTTGAATGAAAAAAAAAGAAAACAAAAAAAT‑TAAAAAAGCCTAATTTTTTTTTAAAACATTAAACATATTATCCTCGGCTATTTTCCAGAATCCCGCGACtttttaacacc                                          >  2:705758/1‑142 (MQ=255)
                                                                                                              ccTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGcc                                          <  1:529974/150‑1 (MQ=255)
                                                                                                              ccTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCa                                         <  2:714236/151‑1 (MQ=255)
                                                                                                               cTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTGTGATTAAAGAAAACAGCAGCAAAAAAAT‑AAAAATTGTCTAAATTTTTTTTTAATAATTAAACATATAATCATCGTTTATTTTTCACAAACCGCCGTCTTTTaaacacc                                          >  1:580833/1‑144 (MQ=255)
                                                                                                               cTTTGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTGGGAAGAAAGAAAAAACAAGCAAAAATAT‑AATAAAGGGCTAATTTTTTTTTACATCATTAAACATATCATCCTCCTCTATTTTTCCGAATACCTCGACTTATaaaaccc                                          >  2:802831/1‑144 (MQ=255)
                                                                                                                tttGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTGGAATGAAAAAAAAAACAAAAAAAAAAAT‑AAAAAAAGCCTAATTTCTTTTAACATCATTAAACATATCATCATCGTCTATTTTTCCGAATACGTCGACTTATAATTGcc                                          >  2:454265/1‑148 (MQ=255)
                                                                                                                 ttGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGa                                       <  1:179223/150‑1 (MQ=255)
                                                                                                                 ttGCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTGGGATGAAAAAAAAAGCAAACAAAAATAT‑AAAAAAGGCCTAATTTTTTTTAAAATCATTAAAAATATCATCATCGTCTATTTTTCACAAACCGTCGACTTATAATCGCCAGa                                       >  2:632105/1‑150 (MQ=255)
                                                                                                                  tGCATAATGCGCCTGGTATTAATTTC‑GGGT‑GCTTTTAGTAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGTCTTATAATCTCCAGat                                      <  1:366492/150‑2 (MQ=255)
                                                                                                                   gCATAATGCGCCTGGGATTAATGGC‑GGGG‑GCTTTTTTGGATGAAAGAAAAAGAAAGCAATAATAT‑AAAAAATGGCTCATTTTTTTTTAAATCATTAAGCATATCCTCATCCGCTCTTTTTCCGAAACCGTCCACCTTTAATCCCCaaa                                       >  2:395363/1‑146 (MQ=255)
                                                                                                                     aTAATGCGCCTGGGATTAATGGCGGGGG‑GGTTTTTGGGA‑GAAAAAAAAAGCAAAAAAAAAAAT‑AATAATAGTCTAATTTTTTTTAAAATAATTAAAAATATAATCAACGTCTATTTTTTCAAATCAGTCCACTTATAATTACCAGaaa                                     >  1:808172/1‑147 (MQ=255)
                                                                                                                       aaTGCGCCTGGGATTAATGGC‑GGGG‑GGTTTTTTGGATGGAAAAAAAAACAAACAAAAAAAT‑AATAATGGCCCAATTTTTTTTTACATCATTAAACATATTATCATCGTCTATTTTTCCGAATACCTCGACTTATAATATCCACAACgg                                   >  2:617223/1‑148 (MQ=255)
                                                                                                                            gcCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAg                           <  1:672233/151‑1 (MQ=255)
                                                                                                                                   aTTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtg                    <  1:488310/151‑1 (MQ=255)
                                                                                                                                                   cTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg      >  1:331422/1‑151 (MQ=255)
                                                                                                                                                    ttttAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTgg     >  2:747855/1‑151 (MQ=255)
                                                                                                                                                      ttAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGc    <  1:120209/150‑1 (MQ=255)
                                                                                                                                                       tAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc  >  1:397489/1‑151 (MQ=255)
                                                                                                                                                        aGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc  <  1:490309/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                      
CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGT‑GCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536047‑4536346

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