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

GATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTG  >  NC_000913/4536050‑4536340
                                                                                                                                                     |                                                                                                                                                
gATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGa                                                                                                                                               <  1:667305/151‑1 (MQ=255)
  tACAGTAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGTTTTTTGGTAt                                                                                                                                             >  2:1190509/1‑151 (MQ=255)
   aCAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAAt                                                                                                                                             <  1:654730/150‑1 (MQ=255)
     aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaa                                                                                                                                          <  1:67470/151‑1 (MQ=255)
     aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaa                                                                                                                                          <  1:500862/151‑1 (MQ=255)
        agtagtAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaaga                                                                                                                                       <  1:1350727/151‑1 (MQ=255)
            gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAAt                                                                                                                                   <  2:106743/151‑1 (MQ=255)
             tcaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGTATTaaaaaaaa                                                                                                                                    >  2:904173/3‑149 (MQ=255)
             taAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAATaaaa                                                                                                                                    >  2:439759/1‑149 (MQ=255)
                aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaa                                                                                                                               <  1:1209008/151‑1 (MQ=255)
                aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGCGCGTTTTTGAATGGAAGAAAAAgca                                                                                                                                >  1:933170/1‑150 (MQ=255)
                aGGGCATATGTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTACTGCTTATAAATCTTCTTGTTTACCTAATCATTATCCTTTGCATAATGCGCCTGGGATTAATGTCTTGTGGTTTTAGGACTGATAGTAAATgaac                                                                                                                               <  2:1097741/151‑4 (MQ=255)
                  ggCATATTTTGCCAATTTTAAGGCGACTTTCCGTTTTAGTGCGGGTCGGCTTTCCTCATCATAAAGCTTATTGTGGACCTAATCTTGTGCTTTTGCGTAATGCTCCTGCAATTAATGACGGGTGCTTTTAGGAATGAAAGAGAATgcaa                                                                                                                               <  2:1091875/149‑1 (MQ=255)
                  ggCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaagc                                                                                                                             <  2:399941/151‑1 (MQ=255)
                    cATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAgcaaaca                                                                                                                            >  1:1286945/1‑150 (MQ=255)
                     atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaagcaa                                                                                                                           <  1:25492/150‑1 (MQ=255)
                          ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataata                                                                                                                      <  1:1027010/150‑1 (MQ=255)
                          ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataata                                                                                                                      <  2:133894/150‑1 (MQ=255)
                          ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataata                                                                                                                      <  1:178986/150‑1 (MQ=255)
                           tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatat                                                                                                                     <  1:1121407/150‑1 (MQ=255)
                           tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatat                                                                                                                     <  1:1261033/150‑1 (MQ=255)
                           tACCAATGTTAAGGAGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTGTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaataattt                                                                                                                     <  2:1108541/150‑3 (MQ=255)
                            aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAACCAataatata                                                                                                                    <  2:469213/150‑1 (MQ=255)
                             ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatataat                                                                                                                  <  2:38279/151‑1 (MQ=255)
                             ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATTGAAAAAAAAgcacaaaaaaaaataa                                                                                                                   >  2:1149040/1‑137 (MQ=255)
                              cAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAAtaatataata                                                                                                                <  1:1018165/151‑1 (MQ=255)
                               aaTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAtaataa                                                                                                               <  1:792900/151‑1 (MQ=255)
                                  gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATgg                                                                                                            <  1:904270/151‑1 (MQ=255)
                                  gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATgg                                                                                                            <  1:445317/151‑1 (MQ=255)
                                   ttAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGAGCTTTTTAGAATGAAAAAAAAAGCAAGCAATAAAAtaataa                                                                                                               >  2:1232374/1‑147 (MQ=255)
                                     aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGTATTAAAGAAAATACAAGCAATAATAtaataa                                                                                                               >  2:19182/1‑145 (MQ=255)
                                     aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAATGAAAAAAAAAACAAACAAAAATATAATAAAggtc                                                                                                          >  1:315024/1‑148 (MQ=255)
                                          gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAAttt                                                                                                    <  1:930379/151‑1 (MQ=255)
                                              ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCt                                                                                                  <  1:661798/149‑1 (MQ=255)
                                              ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATTAAAAAAAAAATAAACAATAAAATAAAAATAGCCTCATtttttt                                                                                                >  2:110439/1‑151 (MQ=255)
                                                     ttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcat                                                                                          <  2:1349327/150‑1 (MQ=255)
                                                       aTTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcatca                                                                                        <  1:170558/150‑1 (MQ=255)
                                                           cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAg                                                                                   <  1:1143179/151‑1 (MQ=255)
                                                                  ccTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGTTATTTTGGATTTAAGAAGAAGCAAACAATAATATAATAAATTCCTTATTTTTTTTAACattatttaaca                                                                                 >  2:324184/1‑141 (MQ=255)
                                                                   cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAGAAAAAGCAAGCAAAAAAAAAATAAAGGGCTTATTTTTTTTAAAATCATTAAGAata                                                                               >  2:1362505/1‑147 (MQ=255)
                                                                   cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCCTTTTTGAATGAAAGAAAATGAAAACAATAAAATAATAATAGCCTAATTTTTTTGAAAATTATTaaaaa                                                                                 >  1:1356853/1‑142 (MQ=255)
                                                                      tCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGCGGTTTTTGGAAAGAAAGAAAATAAAAAAAAAAAAAAAAAAAAGGACTTATTTCTTTgtataaaattaaaca                                                                                 >  1:692195/1‑128 (MQ=255)
                                                                        cTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAACAAACAATAATATAATAAATGGCTAATTTTTTTTAAAAACATTAAACata                                                                               >  1:1219127/1‑142 (MQ=255)
                                                                          gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCAGATCATCATCGTc                                                                    <  2:1101609/151‑1 (MQ=255)
                                                                              ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATGAAAAAAAAAGCAAACAAAAATATAATAAATGGCTAATTTTTTTTTACATCATTTAAAATATCATCATTGTCTAtt                                                                >  1:665815/1‑151 (MQ=255)
                                                                              ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTAtt                                                                <  2:1115101/151‑1 (MQ=255)
                                                                              ataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTAt                                                                 <  1:326083/150‑1 (MQ=255)
                                                                                aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAAAAAAAGCAAACAATAAAATAAAAAAAGGCTAATTTTTTTTAAAATCATTAAACATATCATCATCGTCTAtttt                                                              >  1:833629/1‑150 (MQ=255)
                                                                                aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTg                                                              <  2:683857/151‑1 (MQ=255)
                                                                                   tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTAGGAAGAAAAAAAATGCAAGCAAAAAAAAAAAAAAAGGCTAATTTTTTTTTAAATTATTAAACATATTAATATTGTCTTTTTGcc                                                            >  1:549043/1‑150 (MQ=255)
                                                                                   tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCCTTTTTGAATTTAAAAAAAAACAAGCAAATATTTAATAATTGCCTCATTTTTTTTTAAATCATTTATCATattttcatcacccactttccc                                                           >  2:1339000/1‑139 (MQ=255)
                                                                                   tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGAAAGAAAAAGCAAAAAAAAATATAATAATAGGCTAATTTTTTTTTAAATCAATAAACATAATATCATCCTCTtttt                                                               >  1:404076/1‑147 (MQ=255)
                                                                                      tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGTCCTAATTTCTTTGAACGTCATTAAGCATATCATCAGCGTCTATGTGCCa                                                           <  2:1314894/148‑1 (MQ=255)
                                                                                      tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGa                                                         <  1:640408/150‑1 (MQ=255)
                                                                                       aTTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGaa                                                        <  2:932545/150‑1 (MQ=255)
                                                                                       aTTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGTTTTTTGGGATTAAAGAAAATGCAAGCAAAAATATAATAAAGGCCTAAATTTTTTTAACATCATTAAAAATATCATCATCGTCTAtttttcag                                                          >  1:548839/1‑143 (MQ=255)
                                                                                          gTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGTTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCACTAAGCATATCATCATCGTCTATTTGCCAGAATc                                                      <  1:724506/149‑1 (MQ=255)
                                                                                          gTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATcc                                                     <  1:122733/150‑1 (MQ=255)
                                                                                          gTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAAAAAAAAGAAGCAAAAATAAAATAAATGCCTAATTTTTTTTAACAACATTAAACATAAAATCATTGTCTTTTTTTCAGaa                                                        >  1:142690/1‑147 (MQ=255)
                                                                                           tttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGAATGAAAAAAAAAACAAGCAAAAAAATAATAAAAGCCCAATTTTTTTTAACATCATTTAGCATATCATCATCGTCTATTTTTCAGAAttc                                                     >  1:1212906/1‑147 (MQ=255)
                                                                                            ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGt                                                   <  2:875882/150‑1 (MQ=255)
                                                                                             tACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTc                                                  <  2:308257/150‑1 (MQ=255)
                                                                                               ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATGAAAGAAAAAGCAAGCAAAAAAATAAAAAAGGCCTAATTTTTTTTTACATCATTAAGAATATCATCATCGTCTTTTTTCCCGAAACCGTCGAc                                               >  1:1103314/1‑151 (MQ=255)
                                                                                                cTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCAGTCg                                                 <  2:329139/148‑1 (MQ=255)
                                                                                                cTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTATTTTTTTGATTTTAAAAAAAAACAAGCAAAAAAATAATAAAAGGCTTATTTTTTTTTAAATCATTAAAAATATAATCATTGTCTCTTTTCTAGAATccctaca                                                >  2:252507/1‑144 (MQ=255)
                                                                                                cTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCCTTTTGGGATGAAAAAAAAAGCAAACCAAAATATAAAAAAGGGCTAATTTTTTTTAAAATCATTAAAAATATCATCATCCTCTATTTTTCAGAATACGTCGAc                                               >  1:458421/1‑150 (MQ=255)
                                                                                                 tAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGCTTTTTGGATGAAAAAAAAAAGAAGCAAAAAAATAATAAATGGCCAATTTTTTTTAACAACAATTAACAGATCATCATCGTCTCTTTTCCAGGAtaccgctac                                               >  1:810536/1‑141 (MQ=255)
                                                                                                  aaTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACt                                              <  2:45822/149‑1 (MQ=255)
                                                                                                     cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGGATGAAAAAAAAAGGAAGCAATAATATAATAATAGGCTAATTTTTTTTAACAATATTAAAAATATAATAATCTTCTATTTTTCAGAAAACGCCGACTta                                            >  1:29883/1‑148 (MQ=255)
                                                                                                     cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATaa                                         <  1:65206/151‑1 (MQ=255)
                                                                                                      aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGGATTAAAAAAAAAGCAAACAATAAAATAAAAATAGCCTTATTTTTTTTTAAATCATTAAACATATCATCATCGTCTTTTTTTCAGAAACCGCCGACTTTTaa                                         >  1:970641/1‑150 (MQ=255)
                                                                                                      aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGTTGAAAGAAAAAGCAAGGAAAAATATAATAAAGGCCTAATTTCTTTTTAAAACATTTAAAATATAATCATCGTCTTTTTGTCAGAAACCGTCGACTtat                                           >  2:221961/1‑148 (MQ=255)
                                                                                                            cTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATTAAAAAAAAAGCAAGGAATTAAATAATAAAGGGCTTATTTTTTTTAAAAATATTAAAAAAATCATAATCGTCTATTTTTCAGAATTCGTCGCCTTATAATCTcc                                    >  2:158979/1‑149 (MQ=255)
                                                                                                            cTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAATGAAAAAAAAAGCAAGCAAAAATATAAAAAAAGCCTAATTTTTTTTTACATCATTAAACATATCATCATCGTCTATTTTTCAAAATCCGTCCACTTATAAttccc                                    >  1:994601/1‑145 (MQ=255)
                                                                                                            cTTTGCATAATGCGCCTGGGATTAATGGCGGGCGCTTTTTGGCATTAAAGAAAAATGAAGGAAAAATAAAATAATGGGATCATTTTTTTTAAAAACATTTAAAATATCATCATAGTCTATTTTTCCGAATTCGTCGACTTATAATTGcc                                    >  2:184703/1‑149 (MQ=255)
                                                                                                             tttGCATAATGCGCCTGGGATCAATGGCGGGTGCTTTTATGTATGGAAGAAAATGCCAGCAATAATATAATAAAGGGCTTATTTTTTTGAAAAAAATTAAACAAATCATCAAAGTCTATTTTCCAGAATACGTTGAATTATAAACGc                                     >  2:892678/1‑147 (MQ=255)
                                                                                                               tGCACAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGaa                                <  2:1315425/150‑1 (MQ=255)
                                                                                                                gCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATGAAAAAAAAAGCAAACAAAAATATAATAAAGGCCTAATTTTTTTTAAAATCATTAAGCATATTATCATCGTCTATTTTTCAGAAACCCGCGACCTATAATCGCCAGaa                                >  1:1073876/1‑149 (MQ=255)
                                                                                                                  aTAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAAAGAAAAAAAAAGGAAGCAAAAATATAATAATAGGCTAATTTTTTTTTAAATCATTTAGCATTTCATCATCGTCTATTTTCCAGAAACCGTCGACTTTTAAATGCCCGaa                                >  1:1143377/1‑147 (MQ=255)
                                                                                                                   tAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGAAAAAAAAAAAAACAAACAATAAAATAATAATAGCCTAATTTTTTTTTACAATATTAAACATATCATCATCGTCTTTTTTTCAGAATCCGTCGGCTTATAATTTCCAGAACGGc                            >  1:1103679/1‑150 (MQ=255)
                                                                                                                    aaTGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCat                          <  1:770462/151‑1 (MQ=255)
                                                                                                                          ccTGGGATTAAAGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAg                     <  2:99498/150‑1 (MQ=255)
                                                                                                                             gggATTAATGGCGGGTGGTTTTTTGAAATAAAAAAAAAACCAGGAAATAAATAATAATGGGCTAATTTTTTTTAAAATCATTAAAAAAATAAACATCATCTATTTTCTAGAATCCGTCGACTTATAATTACCAGAAAGGCctataaag                    >  2:59619/1‑144 (MQ=255)
                                                                                                                               gATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGt               <  1:399820/151‑1 (MQ=255)
                                                                                                                                aTTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtg              <  1:588511/151‑1 (MQ=255)
                                                                                                                                aTTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGt               <  2:748699/150‑1 (MQ=255)
                                                                                                                                     tGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCg         >  1:91807/1‑151 (MQ=255)
                                                                                                                                        cGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc      <  1:1233093/151‑1 (MQ=255)
                                                                                                                                        cGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc      <  1:430662/151‑1 (MQ=255)
                                                                                                                                              gCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtg  <  2:38530/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTG  >  NC_000913/4536050‑4536340

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