Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,536,1970AT42.5% 111.6 / 74.8 87intergenic (+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/46);  new base T (37/0);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 2.01e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536047‑4536346
                                                                                                                                                        |                                                                                                                                                      
cATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTa                                                                                                                                                        <  2:941285/151‑1 (MQ=255)
cATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTa                                                                                                                                                        <  1:117286/151‑1 (MQ=255)
   gATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATTATGCGCCTGGGATTAATGGCGAGTGCTTTTAgg                                                                                                                                                      <  1:103003/150‑1 (MQ=255)
       cAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATga                                                                                                                                                 <  2:283587/151‑1 (MQ=255)
        aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaa                                                                                                                                                <  1:524986/151‑1 (MQ=255)
         ggAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaa                                                                                                                                               <  2:349387/151‑1 (MQ=255)
              agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATgaaaaaaa                                                                                                                                           >  2:432258/1‑150 (MQ=255)
              agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaagaaa                                                                                                                                           <  2:247658/150‑1 (MQ=255)
              agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaagaaa                                                                                                                                           <  1:579806/150‑1 (MQ=255)
               gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGaaaaaaaa                                                                                                                                          >  2:892388/1‑150 (MQ=255)
                 aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATg                                                                                                                                        <  1:88174/150‑1 (MQ=255)
                  aaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgca                                                                                                                                      <  2:1365099/151‑1 (MQ=255)
                   aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATGaaaaaaaaatga                                                                                                                                      >  2:810306/1‑146 (MQ=255)
                   aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaa                                                                                                                                     <  2:371444/151‑1 (MQ=255)
                      gCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaagca                                                                                                                                  <  2:510421/151‑1 (MQ=255)
                       cATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGGTTTTTGAATGaaaaaaaaaggaaaca                                                                                                                                  >  2:712852/1‑142 (MQ=255)
                        atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGTATTTAAAAAAAAgcaagaa                                                                                                                                  >  1:303238/1‑147 (MQ=255)
                          atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaataa                                                                                                                              <  1:91360/151‑1 (MQ=255)
                          atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATGAAAAAAAAAgaaaacaa                                                                                                                                 >  2:354189/1‑148 (MQ=255)
                                 cAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAAtaatataata                                                                                                                       <  1:355353/151‑1 (MQ=255)
                                  aaTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAtaataa                                                                                                                      <  2:12137/151‑1 (MQ=255)
                                     gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATg                                                                                                                   <  2:527824/150‑1 (MQ=255)
                                     gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGCGCTTTTTTCGATTAGAGAAAAATCCAGGAATAATAtaataaat                                                                                                                    >  1:291892/1‑148 (MQ=255)
                                       tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCCTTTTTGAAATGAAAAAAATGGAAAAAAAAATAaaaaaa                                                                                                                      >  2:382407/1‑146 (MQ=255)
                                       tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGAAAAAAAAAACAAAAAAAAATAtaaaaaaaga                                                                                                                  >  2:1020036/1‑146 (MQ=255)
                                        aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCt                                                                                                               <  1:773475/151‑1 (MQ=255)
                                         aGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTa                                                                                                              <  2:1000909/151‑1 (MQ=255)
                                         aGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAAAAAAAACAAGCAATAATATAATAAAGGCCt                                                                                                               >  1:1195765/1‑150 (MQ=255)
                                             gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGTATTAAAGAAAAAGAAAACAATAATATAAAAAATGCCTTAttt                                                                                                          >  1:12321/1‑151 (MQ=255)
                                              aCTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAGAAAAAGCAAGAAATAATATAATAATAGGCTAAttt                                                                                                          >  2:302387/1‑150 (MQ=255)
                                              aCTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAAttt                                                                                                          <  1:483700/150‑1 (MQ=255)
                                              aCTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTc                                                                                                         <  1:1208773/151‑1 (MQ=255)
                                              aCTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAACAAAAAAAAATATAATAAAAGCCTAAttt                                                                                                          >  2:991996/1‑150 (MQ=255)
                                               cTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGTAATGTAAAAAAAAGCAAAAAATAATATAATAATGGGCTAAttttt                                                                                                        >  1:1091070/1‑149 (MQ=255)
                                                tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGATTTTAAAAAAAAACAAGCAATAAAATAATAATTGGCTAAtttttt                                                                                                       >  1:903607/1‑151 (MQ=255)
                                                 ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCtt                                                                                                       <  1:629884/150‑1 (MQ=255)
                                                  tCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCGTTTTGGGATGAAAAAAAAAACAAAGAAAAATATAAAAAAAGGCTAATtttttt                                                                                                      >  2:616985/1‑150 (MQ=255)
                                                   ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGa                                                                                                    <  1:1057164/151‑1 (MQ=255)
                                                       tttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCGTTTTTGGATGAAAGAAAAAGCAAGAAAAAATATAAAAAAAGCCTAAATTTTTTTaaaa                                                                                                 >  2:704318/1‑148 (MQ=255)
                                                         tATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcatc                                                                                               <  2:1251175/150‑1 (MQ=255)
                                                         tATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGTTTTTTTGAATGAAAAAAAAATAAAACAAAAATATAATAAATGTCTAATTTTTTTTaaaaac                                                                                               >  1:1249668/1‑146 (MQ=255)
                                                           ttGCGTGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGGATGGAAGAAAAAGAAAACAATAAAATAATAATGGGCTAATTTTTTTTTAAATCAtt                                                                                            >  2:1122938/1‑151 (MQ=255)
                                                           ttGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTAATGAAAGAAAAAAAAAGCAATAATATAATAATGACCTAATtttttttttaatcatt                                                                                            >  1:806795/1‑141 (MQ=255)
                                                            tGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCAtt                                                                                            <  2:561088/150‑1 (MQ=255)
                                                                     ccTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCata                                                                                     <  2:388444/148‑1 (MQ=255)
                                                                         tCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatca                                                                               <  1:932621/150‑1 (MQ=255)
                                                                            tGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGTATGAAAAAAAAAAAAAGCAAAAAAATAATAAATGGCTTATTTTTTTTAACATCCTTAAACa                                                                                       >  1:1211075/1‑139 (MQ=255)
                                                                             gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCg                                                                            <  1:513831/149‑1 (MQ=255)
                                                                              aTTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTc                                                                          <  1:1323642/150‑1 (MQ=255)
                                                                              aTTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCt                                                                         <  2:115059/151‑1 (MQ=255)
                                                                                tataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGTATTAAAAAAAACACAAACAATAAAATAAAAAAAGGCTTATTTTTTTTGACATAAttttaaa                                                                                       >  1:951749/1‑130 (MQ=255)
                                                                                 ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTTCTTTTTTGATTGAGAAAAAAAAGAAACAAGAAAAAAAAAAAGGCCAAATTTTTTTTAAAATCATTAAATATAACATCATCctctct                                                                       >  1:1314778/1‑148 (MQ=255)
                                                                                  taAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTTATGAAAAAAAAAACAAGCAATAATATAAAAATAGGCTTATTTTTTTTTACATTATTAAACATATAATCATCGTCTAttt                                                                     >  1:530034/1‑151 (MQ=255)
                                                                                   aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTAttt                                                                     <  1:1225137/150‑1 (MQ=255)
                                                                                   aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGCTTTTGGGATGAAAAAAAAAAGAAAAAATAAAATAAAAAAAGGCCAATTTTTTTCTAAATCATTAAACATATTATCATTGTATCttt                                                                     >  2:870371/1‑150 (MQ=255)
                                                                                     aTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAAAAAAAGCAAGCAAAAAAATAAAAAAAGCCTAATTTTTTTTAACATCATTAAACATAACATTATCGTCTAtttt                                                                    >  2:680549/1‑148 (MQ=255)
                                                                                      tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGAAAGAAAAAGCAAACAAAAATATAATAATAGGCCAATTTTTTTTAAAAACCTTTAAGATATCATCATCGTCTAtttttc                                                                  >  2:1410709/1‑147 (MQ=255)
                                                                                       cTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATGAAAGAAAAAACAAGCAAAAAAATAATAAAGGCCTAATTTTTTTTTACATCATTAAACAAATCATCATCCTCTAtttt                                                                    >  2:985658/1‑146 (MQ=255)
                                                                                         tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAg                                                                <  2:1161286/149‑1 (MQ=255)
                                                                                          aTTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAAt                                                             <  2:1039434/151‑1 (MQ=255)
                                                                                           ttGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGAAAAAAAAAAAAAGAAAAAATATAATAAATGCCTAATTTTTTTTTAAATCATTAAACATATAATCATCGTCTTTTTTCCCGaa                                                              >  2:277080/1‑149 (MQ=255)
                                                                                           ttGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGGAAAAAAAAGAAAGAAAAAAAATAATAATTAACTAAATTTTTTTTAAATCATTTAAAAAATAATCATCTTCTATTTTCCAGAAt                                                             >  2:889272/1‑150 (MQ=255)
                                                                                            tGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATc                                                            <  2:1307938/150‑1 (MQ=255)
                                                                                             gTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATcc                                                           <  2:202062/150‑1 (MQ=255)
                                                                                               ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGt                                                         <  2:200545/150‑1 (MQ=255)
                                                                                                 aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGGATGAAAAAAAAAGGAAGCAAAAAAATAATAAAAGGCTAATTTTTTTTTACAACATTAAAAATATAATCATCGTCTATTTTTCAGAAACCGTc                                                        >  2:486459/1‑149 (MQ=255)
                                                                                                      aTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAAAGAAAAAAAAAGCAAGCAAAAAAATAAAAAATGGCTAATTTTTTTTTACATTATTTAAACTATTATCCTCGTCTATTTTTCCGAATCCCTCGGCCtat                                                 >  2:1181727/1‑151 (MQ=255)
                                                                                                         aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAAt                                              <  2:34650/151‑1 (MQ=255)
                                                                                                           tttCCTTTGCATAATGCGCCTGGGATTAGTGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATc                                             <  1:727779/150‑1 (MQ=255)
                                                                                                           tttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCGTTTTTGGATTTAAGAAAAATTAAGGAAAAATATAAAAAAAGCCTAATTTTTTTTTAAATCATTTATCATATTATTAACCTCTATTTTGCAGAAACCCTCCGCTTATaa                                               >  2:1254797/1‑148 (MQ=255)
                                                                                                               cTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGCTTTTTGGAAGAAAAAAAAAACAAACAAAAATATAAAAAAGGCCCAAATTTTTTTAAAATCATTAAACATATAATCATCGTCTTTTTTTCAGAAACCCTCGACTTATAATCGcc                                          >  2:1205891/1‑149 (MQ=255)
                                                                                                               cTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAAAAAAAACAAAAAAAAAAATAAAAAAAGGCCAATTTTTTTTTACATCAATTAACATATTATCATCCTCTTTTTTTCAGAAACCCTCGACTTTTAATcccc                                          >  2:463444/1‑149 (MQ=255)
                                                                                                                 ttGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGa                                       <  1:436733/150‑1 (MQ=255)
                                                                                                                  tGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGaa                                      <  2:601528/150‑1 (MQ=255)
                                                                                                                   gCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAAc                                     <  1:1190275/150‑1 (MQ=255)
                                                                                                                    cATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAGGAAAAAGCGAGCAAAAATATAATAATGGCCTAATTTCTTTTTACATCATTAAGAATATTATCATCGTCTCTTTTTCAGCAACCGTCGACTTTTAAACAACAGAACgg                                   >  1:742207/1‑151 (MQ=255)
                                                                                                                       aaTGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCAGTAAGCATATCAACATCGTCTATTTGCCAGAATCCGGAGACTTATAATCGCCAGAATGGCa                                 <  1:1271436/150‑1 (MQ=255)
                                                                                                                       aaTGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAAAGAAAGAAAAAAAAAGCAAAAATATAATAATGGGCTCATTTTTTTTAAAATCATTAAACATATCATTATCCTCTATTTTCCCGAATCCGTCGACTTATAAAAGGCAGAAAGGCa                                 >  2:1066706/1‑150 (MQ=255)
                                                                                                                        aTGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCata                               <  2:896825/151‑1 (MQ=255)
                                                                                                                          gcgcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCatat                              <  2:610827/150‑1 (MQ=255)
                                                                                                                           cgcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCa                            <  1:870221/151‑1 (MQ=255)
                                                                                                                                gggATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACg                       <  2:74779/151‑1 (MQ=255)
                                                                                                                                 ggATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACgg                      <  1:1017255/151‑1 (MQ=255)
                                                                                                                                             gggTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACg           >  2:1219376/1‑151 (MQ=255)
                                                                                                                                                tGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCt         >  2:1254873/1‑150 (MQ=255)
                                                                                                                                                     tttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGc    >  2:1278328/1‑151 (MQ=255)
                                                                                                                                                      ttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGc    <  1:805666/150‑1 (MQ=255)
                                                                                                                                                       tAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc  >  2:347203/1‑151 (MQ=255)
                                                                                                                                                       tAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCAAACGCTGTGGcc   <  1:948615/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                      
CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCC  >  NC_000913/4536047‑4536346

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