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

CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536047‑4536346
                                                                                                                                                        |                                                                                                                                                      
cATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTa                                                                                                                                                        <  1:367532/151‑1 (MQ=255)
   gATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGa                                                                                                                                                     <  1:305880/151‑1 (MQ=255)
    aTACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGaa                                                                                                                                                    <  1:619202/151‑1 (MQ=255)
    aTACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTggg                                                                                                                                                     >  2:490564/1‑149 (MQ=255)
       cAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATga                                                                                                                                                 <  2:90449/151‑1 (MQ=255)
        aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaa                                                                                                                                                <  2:346392/151‑1 (MQ=255)
            gtagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCGTTTTTGAAt                                                                                                                                                   >  1:288142/1‑144 (MQ=255)
               gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAAt                                                                                                                                         <  2:416140/151‑1 (MQ=255)
                   aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaa                                                                                                                                     <  1:211957/151‑1 (MQ=255)
                   aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaa                                                                                                                                     <  2:873546/151‑1 (MQ=255)
                       cATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGGTTTTGGAATTAAAGAAAAAgcaa                                                                                                                                     >  2:27028/1‑147 (MQ=255)
                          atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGAATAaaaaaaaaacaaacaaaa                                                                                                                                >  1:713890/1‑140 (MQ=255)
                          atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGGATGAAAGAAAAAgcaaacaaaa                                                                                                                               >  1:144543/1‑148 (MQ=255)
                          atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaata                                                                                                                               <  2:293458/150‑1 (MQ=255)
                            tttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATTTAAGAAAAAGCCAGCaataatt                                                                                                                            >  2:543513/1‑150 (MQ=255)
                             ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataata                                                                                                                            <  1:481923/150‑1 (MQ=255)
                                   aTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATaataat                                                                                                                    <  2:484385/151‑1 (MQ=255)
                                     gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATgg                                                                                                                  <  2:518112/151‑1 (MQ=255)
                                     gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATgg                                                                                                                  <  1:667825/151‑1 (MQ=255)
                                        aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCt                                                                                                               <  2:978205/151‑1 (MQ=255)
                                        aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAAAAAATGGAAAAAAAAATAtaaaaa                                                                                                                     >  2:120271/1‑145 (MQ=255)
                                           gCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATTAAAAAAAAAACAAGCAATAATAtaataaatggctta                                                                                                             >  1:404013/1‑142 (MQ=255)
                                             gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAGAAAAAGGAAACAATAAAATAATAAATGGCTAAttt                                                                                                          >  2:591765/1‑151 (MQ=255)
                                             gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAAttt                                                                                                          <  1:427697/151‑1 (MQ=255)
                                                tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGCTTTTTGTATTGAAAAAAATAAAAAGCAATATATAATAAAggggcaatttttt                                                                                                       >  1:503420/1‑140 (MQ=255)
                                                tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAGGAAGCAATAATATAAAAAATGCCTAAtttttt                                                                                                       >  2:17629/1‑151 (MQ=255)
                                                   ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCATAATATAATAATGGCCTAATTTCTTTGa                                                                                                    <  2:991103/150‑1 (MQ=255)
                                                   ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAAACTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGa                                                                                                    <  2:288271/151‑1 (MQ=255)
                                                    cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGaa                                                                                                   <  1:473590/151‑1 (MQ=255)
                                                     gTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAGCAAAGAATAATATAATAATAGGCTAATTTCttttta                                                                                                   >  2:882580/1‑147 (MQ=255)
                                                       tttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAAAAAGCAAAAATATAAAAAATGCCCAATttttttttaca                                                                                                 >  2:708148/1‑145 (MQ=255)
                                                           ttGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAGAAAAAACAAGGAAAAATATAATAAAAGCCCAATttttttttaca                                                                                                 >  2:76843/1‑141 (MQ=255)
                                                            tGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCAtt                                                                                            <  2:208328/150‑1 (MQ=255)
                                                             gCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTa                                                                                           <  2:220362/150‑1 (MQ=255)
                                                              cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGTGTTTTTGGATGAAAAAAAAAGCAAGCAAAAAAATAAAAAAAGCCCTAATtttttttacatta                                                                                               >  2:840281/1‑139 (MQ=255)
                                                                aGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCa                                                                                       <  2:239356/151‑1 (MQ=255)
                                                                aGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCa                                                                                       <  2:582503/151‑1 (MQ=255)
                                                                 gTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCat                                                                                      <  1:58756/151‑1 (MQ=255)
                                                                 gTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAGAAAAAGCAAAAAATAATATAATAAAAGCCTAATTTTTTTTTAAATCATTAAACa                                                                                       >  2:941589/1‑150 (MQ=255)
                                                                   cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCata                                                                                     <  1:814121/150‑1 (MQ=255)
                                                                     ccTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatca                                                                                  <  2:377368/151‑1 (MQ=255)
                                                                     ccTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAAAGAAAAAAAAAGCAAACAATAATATAATAAATGCCTTATTTTTTTTTAAATCATTAAACa                                                                                       >  2:881608/1‑146 (MQ=255)
                                                                      cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTAGGGATGAAAAAAAAAACAAACAAAAAAATTATAATTGGTTTATTTCTTTTTACATCCTTa                                                                                           >  2:338100/1‑141 (MQ=255)
                                                                        ttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcat                                                                                 <  1:162475/149‑1 (MQ=255)
                                                                          ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAAAGAAAAAAAATGAAAGAAATAATATAATAAATGGCTAATTTTTTTTTACATCATTaaaaaaat                                                                                    >  2:501764/1‑138 (MQ=255)
                                                                          ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATTAAAAAAAAAGCAAGGAAAAATAAAATAATAGGCTAATTTTTTTTTAAATCCTTAAACaaattaaca                                                                               >  2:951841/1‑144 (MQ=255)
                                                                                 ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAAATAAAAAAAAAAAAAGCAATAATAAAAAAAAAGGCTAAATTCTTTTTACATTCTTAAACATAAAATCATCCTCtttt                                                                      >  2:656840/1‑151 (MQ=255)
                                                                                   aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAGAAAGAAAAAAAATAAAAATTGACCAAATTTTTTTAAAATCATTAAACATATCATCATCGTCTAttt                                                                     >  2:175542/1‑150 (MQ=255)
                                                                                     aTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATGAAAAAAAAAGCAAAACAAAAAATAATAAAGGGCTAATTTTTTTTTAAATCATTTAAAATATTATTATCCTCTCttt                                                                     >  2:1018880/1‑148 (MQ=255)
                                                                                         tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGa                                                               <  1:45767/150‑1 (MQ=255)
                                                                                         tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGa                                                               <  2:604851/150‑1 (MQ=255)
                                                                                            tGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATGAAAGAAAAAGCAAGCAATAAAATAAAAAAGGCCCAATTTTTTTTTAAATCAATTAACATATCATCATCGTCTATTTTTCCGAATcc                                                           >  2:868329/1‑151 (MQ=255)
                                                                                               ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGt                                                         <  2:796832/150‑1 (MQ=255)
                                                                                                tACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGATTGTAAGAAAAAGCAAAGAAAAAAATAATAAAGGGGTAATTTTTTTTTAAAACATTTAAGATATCATTTATGTCTATTTGCCaaaaattgtcg                                                       >  1:329270/1‑144 (MQ=255)
                                                                                                   cTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAGCAAACAAAAAAATAAAAAAAGGCTAATTTTTTTTTAAATCATTTAAAAAATAATCATCGTCTATTTTTCAGAAACCg                                                          >  2:672566/1‑145 (MQ=255)
                                                                                                        cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATaa                                               <  1:672607/151‑1 (MQ=255)
                                                                                                        cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAGAAAGCAAAAATATAATAATAGGCTAATTTTTTTTTAAATCATTAAGCATAATATCATCGTCTTTTTTTCAGAATCCGTCGACTtat                                                 >  2:550303/1‑149 (MQ=255)
                                                                                                             tCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGGATGAAAAAAAAACCAAGCAATAAAATAATAAAAGCCCAATTTTTTTTAATATTAATAAACATATCATCCTTGTCTATTTTTCAAAATCCGTCGACCTTTaa                                               >  2:628111/1‑146 (MQ=255)
                                                                                                             tCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATGAAAAAAAAAGAAAACAAAAAAATTAAAAAAGCCTAATTTTTTTTTAAAACATTAAACATATTATCCTCGGCTATTTTCCAGAATCCCGCGACtttttaacacc                                          >  2:705758/1‑142 (MQ=255)
                                                                                                              ccTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGcc                                          <  1:529974/150‑1 (MQ=255)
                                                                                                              ccTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCa                                         <  2:714236/151‑1 (MQ=255)
                                                                                                               cTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGTGATTAAAGAAAACAGCAGCAAAAAAATAAAAATTGTCTAAATTTTTTTTTAATAATTAAACATATAATCATCGTTTATTTTTCACAAACCGCCGTCTTTTaaacacc                                          >  1:580833/1‑144 (MQ=255)
                                                                                                               cTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGAAGAAAGAAAAAACAAGCAAAAATATAATAAAGGGCTAATTTTTTTTTACATCATTAAACATATCATCCTCCTCTATTTTTCCGAATACCTCGACTTATaaaaccc                                          >  2:802831/1‑144 (MQ=255)
                                                                                                                tttGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAATGAAAAAAAAAACAAAAAAAAAAATAAAAAAAGCCTAATTTCTTTTAACATCATTAAACATATCATCATCGTCTATTTTTCCGAATACGTCGACTTATAATTGcc                                          >  2:454265/1‑148 (MQ=255)
                                                                                                                 ttGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGa                                       <  1:179223/150‑1 (MQ=255)
                                                                                                                 ttGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAAAAAAAGCAAACAAAAATATAAAAAAGGCCTAATTTTTTTTAAAATCATTAAAAATATCATCATCGTCTATTTTTCACAAACCGTCGACTTATAATCGCCAGa                                       >  2:632105/1‑150 (MQ=255)
                                                                                                                  tGCATAATGCGCCTGGTATTAATTTCGGGTGCTTTTAGTAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGTCTTATAATCTCCAGat                                      <  1:366492/150‑2 (MQ=255)
                                                                                                                   gCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAGAAAAAGAAAGCAATAATATAAAAAATGGCTCATTTTTTTTTAAATCATTAAGCATATCCTCATCCGCTCTTTTTCCGAAACCGTCCACCTTTAATCCCCaaa                                       >  2:395363/1‑146 (MQ=255)
                                                                                                                     aTAATGCGCCTGGGATTAATGGCGGGGGGGTTTTTGGGAGAAAAAAAAAGCAAAAAAAAAAATAATAATAGTCTAATTTTTTTTAAAATAATTAAAAATATAATCAACGTCTATTTTTTCAAATCAGTCCACTTATAATTACCAGaaa                                     >  1:808172/1‑147 (MQ=255)
                                                                                                                       aaTGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGGAAAAAAAAACAAACAAAAAAATAATAATGGCCCAATTTTTTTTTACATCATTAAACATATTATCATCGTCTATTTTTCCGAATACCTCGACTTATAATATCCACAACgg                                   >  2:617223/1‑148 (MQ=255)
                                                                                                                            gcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAg                           <  1:672233/151‑1 (MQ=255)
                                                                                                                                   aTTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtg                    <  1:488310/151‑1 (MQ=255)
                                                                                                                                                   cTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg      >  1:331422/1‑151 (MQ=255)
                                                                                                                                                    ttttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTgg     >  2:747855/1‑151 (MQ=255)
                                                                                                                                                      ttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGc    <  1:120209/150‑1 (MQ=255)
                                                                                                                                                       tAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc  >  1:397489/1‑151 (MQ=255)
                                                                                                                                                        aGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc  <  1:490309/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                      
CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  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: