Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,509,8800AC42.7% 123.5 / 13.0 76D84A (GAC→GCC) yfeOputative ion channel protein
Reads supporting (aligned to +/- strand):  ref base A (36/7);  new base C (0/32);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.93e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTG  >  NC_000913/2509730‑2510026
                                                                                                                                                       |                                                                                                                                                  
ccTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGa                                                                                                                                                    <  2:680702/151‑1 (MQ=255)
   cGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGAccc                                                                                                                                                 <  2:119089/151‑1 (MQ=255)
     gTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCg                                                                                                                                               >  2:167938/1‑151 (MQ=255)
       aTTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGcc                                                                                                                                             >  2:917186/1‑151 (MQ=255)
          aCAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCtg                                                                                                                                           <  2:582850/150‑1 (MQ=255)
            aGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCtgtg                                                                                                                                         >  1:305224/1‑150 (MQ=255)
             gAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGa                                                                                                                                        >  2:944267/1‑150 (MQ=255)
               aTTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCTTTTCAGCCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAAcc                                                                                                                                     <  2:137342/151‑1 (MQ=255)
                 ttGCTCGGGCAACGTCTGCCGGCAACTCCGGGGATGGCCCAGGGTTCTCCCCTCTGGATCATCGGTGTATAACCGATAACGGGTATTGCGGTGGGTTTGGTTCCCTGTTCCAGCCAGGTCAATGCCGGACCAGCCCCCGCCtgtg                                                                                                                                         <  1:1196961/145‑1 (MQ=255)
                  tGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCt                                                                                                                                  >  1:303633/1‑151 (MQ=255)
                  tGCCCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGc                                                                                                                                   >  1:1275635/1‑150 (MQ=255)
                    ctctGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTg                                                                                                                                 >  1:1292036/1‑150 (MQ=255)
                         gCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGgc                                                                                                                           >  1:616155/1‑151 (MQ=255)
                          cAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCCCCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTTTTGCGGTGGGGTTGGTTATCCGTTTCAGCCGGGGTCAGGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGgc                                                                                                                           <  1:883635/150‑1 (MQ=255)
                           aaCGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGgc                                                                                                                           <  2:249222/149‑1 (MQ=255)
                           aaCGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGgcgc                                                                                                                         >  1:802773/1‑151 (MQ=255)
                           aaCGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGgcgc                                                                                                                         >  2:1060607/1‑151 (MQ=255)
                             cGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCAc                                                                                                                       >  1:154569/1‑151 (MQ=255)
                              gACTGCGGGGACCTCGGGGGATAGCCCAGGATTCCCCCCTCTGGATCATCGGTTTATTACCGCTAACGGGTTTTGGGGTGGGGTTGGTTATCCGTTTCAGCCGGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCa                                                                                                                        <  1:1088817/149‑1 (MQ=255)
                                   cgggtaaCGCTGGGAAAGGCCCGGAATCCACCCTCCTGGATAAGGGGTGTATTAACGCAAACGGGTATGGGGGGGGGGTTGGTTATCCGTTCCGGCCGGGGTAATGCCGAACCGGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGt                                                                                                                   <  1:630809/144‑1 (MQ=255)
                                     gggACCTCTGGGGATACCCCAGGATCCCCCCCCCTGGATAATCGTTGTATTAACGCAAACGGGTATTGCGGTGGGTTGGTTTATCCTTTTCACCCGGGGTCATCCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTcc                                                                                                                <  2:176417/150‑1 (MQ=255)
                                     gggACCTCGGGGGATAGCCCGGGATTCCCCCCTTTGGATCATCGGTGTATTACCGCTACCGGGTTTTGGGGGGGGTTTGGTTCCCCGTTTCAGCCAGGGTCTTGCCGCACAAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTcc                                                                                                                <  1:243434/150‑1 (MQ=255)
                                     gggAACTCTGGGGATAGCCCAGGATTCCCCCCTCTGGATCATCGGTGTATTAACGCTAAGGGTTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCGGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGt                                                                                                                   <  2:802805/147‑1 (MQ=255)
                                     gggAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCg                                                                                                               >  2:940106/1‑151 (MQ=255)
                                     gggAACTCTCGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGAATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCATCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCCCCGGCTc                                                                                                                 >  1:1013728/1‑149 (MQ=255)
                                         aCTCTGGGGATAGCCCAGGATTCACCCCTCTGGTTCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGccc                                                                                                            >  2:674793/1‑150 (MQ=255)
                                           tctGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCctct                                                                                                         >  1:1271516/1‑151 (MQ=255)
                                               gggATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTgcgc                                                                                                     >  1:499706/1‑151 (MQ=255)
                                                    aGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCt                                                                                                >  2:336799/1‑151 (MQ=255)
                                                    aGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCt                                                                                                >  2:894493/1‑151 (MQ=255)
                                                      cccaggtgccACCCCACTGGATCATCGGTTTATTAACGCTAAGGGTTATTGGGGTGGGTTTGGTTTCCCTTTTCACCCAGGGTCATGCGGGACCGGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCt                                                                                                <  1:617643/140‑1 (MQ=255)
                                                      cccAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGTTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCt                                                                                                <  2:306401/149‑1 (MQ=255)
                                                          ggttgccccccTGTGGATATTCGGTGTTTTAACCCTAACGTTTTTTGGGTGGGGGTTGTTTCTCTGTTTCGCCAAGGGTAATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGAc                                                                                            <  1:1084807/142‑1 (MQ=255)
                                                          ggATTCCCCCCTCTGGATCATGGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGTTTATCCGTTTCAGCCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACt                                                                                           <  2:1142947/150‑1 (MQ=255)
                                                          ggATTCCCCCCTCTGGATCATCGTTGTATTAACGCTAACGGGTATGGCGGGGGTTTTGGTTATCCGTTTCACCCGGGGTCAGGCCGCCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACt                                                                                           <  1:1119816/150‑1 (MQ=255)
                                                          ggATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACtt                                                                                          <  2:315569/151‑1 (MQ=255)
                                                            aTTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTa                                                                                         <  2:199706/150‑1 (MQ=255)
                                                              tCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCg                                                                                      >  2:1209410/1‑151 (MQ=255)
                                                               ccccccTCTGGATCATCGGTGTTTAAACGCTAACGGGTATTGGGGTGGGGTTGGTTATCCGTTTCACCCAGGGTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATc                                                                                       <  1:705181/147‑1 (MQ=255)
                                                                 ccccTTGGGATCATCGGTTTATTAGCGCTAAGGGTTTTTGGGGTGGGTTTGTTTATCTGTTTCAGCCGGGTTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAg                                                                                   <  1:657658/151‑1 (MQ=255)
                                                                   ccTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCTTTTCAGCCAGGTTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCa                                                                                 <  1:1071942/151‑1 (MQ=255)
                                                                     tctGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCAtt                                                                               >  2:969920/1‑151 (MQ=255)
                                                                      ctGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCAtt                                                                               >  2:887744/1‑150 (MQ=255)
                                                                           tcatcGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATtct                                                                         >  2:398410/1‑151 (MQ=255)
                                                                             atcGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATtctc                                                                        <  2:1053814/150‑1 (MQ=255)
                                                                               cGGTGTTTAACCGCTAGCGTGTATTGCGGTGGGTTTGGTTTCCCGTTTCACCCGGGGTCTTGCCGCCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTa                                                                              <  1:735600/142‑1 (MQ=255)
                                                                                ggTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTc                                                                    >  2:927343/1‑151 (MQ=255)
                                                                                ggTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTc                                                                    >  1:637019/1‑151 (MQ=255)
                                                                                    tATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGc                                                                >  2:83550/1‑151 (MQ=255)
                                                                                       tAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTgg                                                             >  1:347537/1‑151 (MQ=255)
                                                                                            cTAACGGGTATTGCGGTGGGGTTGGTTATCCTTTTCAGCCGGGGTCATGCCGAACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTggcggc                                                         <  2:1136834/150‑1 (MQ=255)
                                                                                             tAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGt                                                       >  2:1071783/1‑151 (MQ=255)
                                                                                               aCGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTc                                                      >  1:955244/1‑150 (MQ=255)
                                                                                               aCGGATATGGGGGTGGGTTTGTTTATCCATGTCGGCGGGGCTCACGCCGGCCAATCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTggcggc                                                         <  2:876036/147‑1 (MQ=255)
                                                                                                 gggTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGc                                                   >  1:215703/1‑151 (MQ=255)
                                                                                                    tATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTg                                                >  2:1290136/1‑151 (MQ=255)
                                                                                                       tGCGGTGGGTTTGTTTATCCTTTTCAGCCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTggg                                              <  2:191818/150‑1 (MQ=255)
                                                                                                        tgcggggggTTGGTTTATCCTTTTCAGCCGGGGTCATGCCGACCAAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTg                                                <  1:1075361/143‑1 (MQ=255)
                                                                                                        gggggggggTTGTTTACCCCTTTCACCCGGGCCCATGCCGTACCAGGCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgggg                                             <  1:1151035/145‑1 (MQ=255)
                                                                                                        gCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgggg                                             <  2:80417/150‑1 (MQ=255)
                                                                                                          ggTGGGTTTGTTTATCCGTTTCAGCCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGcc                                           <  1:835812/150‑1 (MQ=255)
                                                                                                                 ttgtttATCCTTTTCAGCCGGGTTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACAt                                    <  1:1111283/150‑1 (MQ=255)
                                                                                                                    gtttttcgttttCACCCGGGGTTATCCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCg                                 <  1:1216972/142‑1 (MQ=255)
                                                                                                                    gtttATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATcc                                  <  2:521230/148‑1 (MQ=255)
                                                                                                                      ttATCCTTTTCAGCCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGa                                <  2:96917/150‑1 (MQ=255)
                                                                                                                       tATCCGTTTCAGCCGGGGTCATGCCGCACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGAt                               <  1:887708/150‑1 (MQ=255)
                                                                                                                         tcgttttCAGCCAGAGTCAGCCCGGACAAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCa                             <  1:601956/146‑1 (MQ=255)
                                                                                                                          ccGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCAt                            >  2:197007/1‑150 (MQ=255)
                                                                                                                                   ccAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAat                  >  2:180277/1‑151 (MQ=255)
                                                                                                                                   ccAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAat                  >  2:1195790/1‑151 (MQ=255)
                                                                                                                                      gggTCATGCCGGACCGTCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCg                       <  1:830278/143‑1 (MQ=255)
                                                                                                                                      gggTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATc               >  1:553810/1‑151 (MQ=255)
                                                                                                                                      gggTCATGCCGACCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAat                  <  2:499727/148‑1 (MQ=255)
                                                                                                                                       ggTCAGGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCg              <  2:1292106/151‑1 (MQ=255)
                                                                                                                                              ccGGACCGCCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCt          <  1:976498/148‑1 (MQ=255)
                                                                                                                                                ggACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCgg     >  1:359086/1‑151 (MQ=255)
                                                                                                                                                   ccAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTg  >  2:1278612/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
CCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTG  >  NC_000913/2509730‑2510026

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