Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,624 T→G 5.3% intergenic (+170/+167) codA → / ← cynR cytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,6240TG 5.3% 1765.9 / 15.1 612intergenic (+170/+167)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (243/321);  new base G (3/44);  total (247/365)
Fisher's exact test for biased strand distribution p-value = 9.40e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 623 total aligned reads displayed.

TTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCG  >  NC_000913/357475‑357773
                                                                                                                                                     |                                                                                                                                                      
  aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATtc                                                                                                                                                      <  1:9143567/150‑1 (MQ=38)
  aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATtc                                                                                                                                                      >  2:4811746/1‑150 (MQ=38)
   gTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGAACTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCAttaa                                                                                                                                                     >  1:10678989/1‑148 (MQ=21)
      tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATtcctcc                                                                                                                                                  >  1:1340391/1‑150 (MQ=40)
      tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGAGGCTACCGTGATTCAATCCCGCCAACAACCGCGCAGtcctcc                                                                                                                                                  >  2:11050393/1‑150 (MQ=255)
       aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCa                                                                                                                                                 <  2:4385294/150‑1 (MQ=255)
        tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCaa                                                                                                                                                <  1:8230541/150‑1 (MQ=255)
         gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc                                                                                                                                               <  2:199855/150‑1 (MQ=255)
         gCCGGATGCGGCGTGAACGCCATATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc                                                                                                                                               >  1:1851104/1‑150 (MQ=255)
            ggATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGcc                                                                                                                                            <  2:9242671/150‑1 (MQ=255)
            ggATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGcc                                                                                                                                            <  2:10569133/150‑1 (MQ=255)
             gATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAACCGCGCATTACTCCAACGCCa                                                                                                                                           >  1:3802335/1‑150 (MQ=255)
                 cggagTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCAtgtg                                                                                                                                       <  2:8405493/150‑1 (MQ=255)
                   gcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCa                                                                                                                                     >  1:5637375/1‑150 (MQ=255)
                        aaCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAt                                                                                                                                <  1:8865446/150‑1 (MQ=17)
                        aaCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAAt                                                                                                                                <  2:3504821/150‑1 (MQ=255)
                        aaCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAAt                                                                                                                                >  1:6312025/1‑150 (MQ=255)
                         aCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCATCCAACGCCATGTGCAAAAATg                                                                                                                               >  1:5837454/1‑150 (MQ=37)
                              ttATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTc                                                                                                                          <  2:10882292/150‑1 (MQ=255)
                                   cGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTAATCCAACGCCATGTGCAAAAATGCCTTCGCTGc                                                                                                                     >  1:792771/1‑150 (MQ=255)
                                   cGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAAACGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAga                                                                                                                     >  1:10261667/1‑149 (MQ=255)
                                   cGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCTGc                                                                                                                     >  1:6071097/1‑150 (MQ=255)
                                     gCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCgg                                                                                                                   <  2:6374537/150‑1 (MQ=255)
                                         aCGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGt                                                                                                               <  1:3326290/150‑1 (MQ=255)
                                           gTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCt                                                                                                            <  2:10301219/150‑1 (MQ=255)
                                           gTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCt                                                                                                            <  2:629268/150‑1 (MQ=255)
                                           gTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCATCCAACGCCATGTGCAAAAATGCCTTCGCTGCGGCTGTCt                                                                                                            >  1:9892777/1‑150 (MQ=255)
                                                gCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAg                                                                                                       <  1:5368187/150‑1 (MQ=255)
                                                 cACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGc                                                                                                      >  2:3011534/1‑150 (MQ=255)
                                                   cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTg                                                                                                    <  1:4309601/150‑1 (MQ=255)
                                                      aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAg                                                                                                 <  1:9408783/150‑1 (MQ=255)
                                                           gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTAt                                                                                            <  1:2013678/150‑1 (MQ=255)
                                                           gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATCCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTAt                                                                                            <  2:5284953/150‑1 (MQ=255)
                                                             aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGc                                                                                          <  1:6605819/150‑1 (MQ=255)
                                                              ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGcc                                                                                         <  2:9030586/150‑1 (MQ=255)
                                                              ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGcc                                                                                         <  2:1851104/150‑1 (MQ=255)
                                                              ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGcc                                                                                         <  2:451154/150‑1 (MQ=255)
                                                              ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGcc                                                                                         <  2:4703792/150‑1 (MQ=255)
                                                               gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCg                                                                                        >  2:2423347/1‑150 (MQ=255)
                                                               gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCg                                                                                        <  2:1520378/150‑1 (MQ=255)
                                                               gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCCTGTGCAAAAATGCCGACCCGGCGGCGGTCGGCCAGCTGAAGTGTCTGCCg                                                                                        >  2:5559493/1‑150 (MQ=14)
                                                                 cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGa                                                                                      <  2:3982543/150‑1 (MQ=255)
                                                                 cTGATAAGCGGCGCGCAGCAGGCAATTCCAGCCGCGGATCTGGGTCAGCGGCTACCGGGCGTCATTCCCGCCAACAACCGCGCAGACCGCCAACGCCAGGTGCAAAAATGCCTGCGCCGCGGCTGTCTCCCAGCTGTCGGTGATGCCCGa                                                                                      >  2:5715766/1‑150 (MQ=255)
                                                                       aGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcggc                                                                                <  1:9514895/150‑1 (MQ=255)
                                                                        gCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcggcg                                                                               <  2:6914792/150‑1 (MQ=255)
                                                                         cGTAGCGCATCAGGCAATTCCCGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGt                                                                              <  2:8300826/150‑1 (MQ=255)
                                                                         cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCt                                                                                                                 >  2:7062214/1‑116 (MQ=255)
                                                                         cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGt                                                                              >  1:846735/1‑150 (MQ=255)
                                                                          gTAGCTCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGTTTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTg                                                                             <  2:8064344/150‑1 (MQ=255)
                                                                          gTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTg                                                                             <  1:4018988/150‑1 (MQ=255)
                                                                           tAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGa                                                                            <  2:5130294/150‑1 (MQ=255)
                                                                                  tCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                     <  1:6139939/150‑1 (MQ=255)
                                                                                  tCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                     >  1:3195001/1‑150 (MQ=255)
                                                                                      gCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                 >  2:3644190/1‑150 (MQ=255)
                                                                                      gCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTACTCCAACGCCATGTGCAAAAAAGCCTTAGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                 >  1:153401/1‑150 (MQ=255)
                                                                                      gCAAGTCCAGCCGCTGAGCTGTGGCAGCGGCTCCCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCCACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGGAGGGTAAGCGGGAGGCGCCGGGAGCGCCTTATcc                                                                 >  2:6562786/1‑150 (MQ=17)
                                                                                       cAATTCCAGCCGCTGATCTGTGTCAGCGGCGACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                <  2:3389883/150‑1 (MQ=255)
                                                                                        aaTTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCAGTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                               <  1:2113681/150‑1 (MQ=255)
                                                                                             cAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                          <  1:10002069/150‑1 (MQ=255)
                                                                                              aGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                                         >  1:5122872/1‑150 (MQ=255)
                                                                                               gCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGt                                                        >  2:3773872/1‑150 (MQ=255)
                                                                                                ccGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa                                                       <  1:4685120/150‑1 (MQ=255)
                                                                                                ccGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCCACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa                                                       <  2:2748928/150‑1 (MQ=255)
                                                                                                 cGCTGATCTGTGTCCGCGGCTACCGAGATTCATTCCCGCCAACAACCGCGCAGTCCGCCAACGCCATGTGCAAAAATGGCTTCGCCTCTGCTTTATTCCTGGCGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTag                                                      <  2:2315682/150‑1 (MQ=21)
                                                                                                 cGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTag                                                      <  1:10869396/150‑1 (MQ=255)
                                                                                                  gCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCTACGCCATGTGCAAAAATGCCTTCGCTGCGGCTGTCTGCCTGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTaga                                                     <  2:792483/150‑1 (MQ=255)
                                                                                                   cTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTagag                                                    >  1:10969689/1‑150 (MQ=255)
                                                                                                   cTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTagag                                                    <  2:6061843/150‑1 (MQ=255)
                                                                                                     gATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATCCCCGATCCTGCCTCAAAGCCATATCCGTCCTACGTAGACCa                                                  >  1:7011848/1‑150 (MQ=18)
                                                                                                       tCTGGGTCAGCGGCTACCGTGAGTCCTTCCCGCCAACACCCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                                                <  2:7678900/150‑1 (MQ=255)
                                                                                                          gtgcCAGCGGCTACCGTGATTCAGTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGaa                                             <  1:5214003/150‑1 (MQ=255)
                                                                                                                gCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTa                                       <  1:2559822/150‑1 (MQ=255)
                                                                                                                   gCTACCGTGATTCAGTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGc                                    <  1:1132869/150‑1 (MQ=255)
                                                                                                                     tACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAATATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCt                                  <  2:5883241/150‑1 (MQ=255)
                                                                                                                     tACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAAAAAGCCTTCGCAGCGGCGGTCCGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCt                                  <  2:1207751/150‑1 (MQ=255)
                                                                                                                     gaCCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCt                                  <  1:9761600/149‑1 (MQ=255)
                                                                                                                       ccGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCCACGCCATGTGCAACAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGa                                <  1:4481960/150‑1 (MQ=255)
                                                                                                                        cGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGAt                               >  2:5389360/1‑150 (MQ=255)
                                                                                                                        cGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATACGGCCTACGTAGAGCCCTGAACTCGTAGGCCTGAt                               >  1:5005846/1‑150 (MQ=255)
                                                                                                                          tGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATaa                             >  1:1465453/1‑150 (MQ=255)
                                                                                                                          tGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTAATCCCGCAAGCGGCGTGAACGCCATAGCCGGCCGACGTCGAGCACTGAACTCCTAGGCCACATaa                             >  2:10106712/1‑150 (MQ=21)
                                                                                                                           caaacattcCCGCCAACAACCGCGCATTCCTCCAACGCCATCTGCAAAAATGCCTTCCCAGCGGCTGTCAGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAg                            <  2:7375050/146‑1 (MQ=255)
                                                                                                                            cttccatcCCGACAACACCCGCGCATGCGCCCAACGCCATGTGCAAAAACGCGTTCGCAGCGGGGGTCGGCCTGCAGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGc                           <  2:874048/144‑1 (MQ=18)
                                                                                                                               cattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTa                        <  1:8744428/150‑1 (MQ=255)
                                                                                                                                 ttcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGGCTTCGCATGGGCTTTCTGCCTGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgc                      <  2:4132207/150‑1 (MQ=25)
                                                                                                                                 ttcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgc                      <  2:977487/150‑1 (MQ=255)
                                                                                                                                   cccGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgcgc                    <  1:7570072/150‑1 (MQ=255)
                                                                                                                                     cGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATg                                                                                    >  2:6785359/1‑84 (MQ=255)
                                                                                                                                     cGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATg                                                                                    <  1:6785359/84‑1 (MQ=255)
                                                                                                                                       ccATCAACCGCGCAGCCCCCACACGCCATGTGCAAAAATGGCTTCGCCGCGGGTTTCCCGCGCGCGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCa                <  2:1391867/150‑1 (MQ=14)
                                                                                                                                        caacaaCCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTATTTTATGACGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGACCTCGTAGGCCTGATCAGCGTGGCGCATCAg               >  2:3164044/1‑150 (MQ=25)
                                                                                                                                        caacaaCCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAg               <  1:3622471/150‑1 (MQ=255)
                                                                                                                                         aacaacCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAgg              <  1:2423347/150‑1 (MQ=255)
                                                                                                                                            aacCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTGCGTAGAGCACTGACCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa           >  2:9390206/1‑150 (MQ=255)
                                                                                                                                            aacCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa           >  2:5035111/1‑150 (MQ=255)
                                                                                                                                             acCGCGCAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAAt          <  2:9632988/150‑1 (MQ=255)
                                                                                                                                                 cgcATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa      >  1:295301/1‑150 (MQ=255)
                                                                                                                                                 cgcAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa      <  2:726944/150‑1 (MQ=255)
                                                                                                                                                    aTTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTAGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGaa                                             <  2:2838338/108‑1 (MQ=255)
                                                                                                                                                     ttCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCg  >  1:7356915/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCG  >  NC_000913/357475‑357773

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: