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.

TTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCG  >  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)
                                           gTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑Ct                                                                                                            <  2:10301219/150‑1 (MQ=255)
                                           gTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑Ct                                                                                                            <  2:629268/150‑1 (MQ=255)
                                           gTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCATCCAACGCCATGTGCAAAAATGCCTTCGCTGCGGCTGT‑Ct                                                                                                            >  1:9892777/1‑150 (MQ=255)
                                                gCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAg                                                                                                       <  1:5368187/150‑1 (MQ=255)
                                                 cACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGc                                                                                                      >  2:3011534/1‑150 (MQ=255)
                                                   cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTg                                                                                                    <  1:4309601/150‑1 (MQ=255)
                                                      aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAg                                                                                                 <  1:9408783/150‑1 (MQ=255)
                                                           gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTAt                                                                                            <  1:2013678/150‑1 (MQ=255)
                                                           gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATCCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTAt                                                                                            <  2:5284953/150‑1 (MQ=255)
                                                             aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGc                                                                                          <  1:6605819/150‑1 (MQ=255)
                                                              ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGcc                                                                                         <  2:9030586/150‑1 (MQ=255)
                                                              ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGcc                                                                                         <  2:1851104/150‑1 (MQ=255)
                                                              ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGcc                                                                                         <  2:451154/150‑1 (MQ=255)
                                                              ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGcc                                                                                         <  2:4703792/150‑1 (MQ=255)
                                                               gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCg                                                                                        >  2:2423347/1‑150 (MQ=255)
                                                               gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCg                                                                                        <  2:1520378/150‑1 (MQ=255)
                                                               gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCCTGTGCAAAAATGCCGACCCGGCGGCGGT‑CGGCCAGCTGAAGTGTCTGCCg                                                                                        >  2:5559493/1‑150 (MQ=14)
                                                                 cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGa                                                                                      <  2:3982543/150‑1 (MQ=255)
                                                                 cTGATAAGCGGCGCGCAGCAGGCAATTCCAGCCGCGGATCTGGGTCAGCGGCTACCGGGCGTCATTCCCGCCAACAACCGCGCAGACCGCCAACGCCAGGTGCAAAAATGCCTGCGCCGCGGCTGT‑CTCCCAGCTGTCGGTGATGCCCGa                                                                                      >  2:5715766/1‑150 (MQ=255)
                                                                       aGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATgcggc                                                                                <  1:9514895/150‑1 (MQ=255)
                                                                        gCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATgcggcg                                                                               <  2:6914792/150‑1 (MQ=255)
                                                                         cGTAGCGCATCAGGCAATTCCCGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGt                                                                              <  2:8300826/150‑1 (MQ=255)
                                                                         cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCt                                                                                                                 >  2:7062214/1‑116 (MQ=255)
                                                                         cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGt                                                                              >  1:846735/1‑150 (MQ=255)
                                                                          gTAGCTCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGTTTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTg                                                                             <  2:8064344/150‑1 (MQ=255)
                                                                          gTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTg                                                                             <  1:4018988/150‑1 (MQ=255)
                                                                           tAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGa                                                                            <  2:5130294/150‑1 (MQ=255)
                                                                                  tCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                     <  1:6139939/150‑1 (MQ=255)
                                                                                  tCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                     >  1:3195001/1‑150 (MQ=255)
                                                                                      gCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                 >  2:3644190/1‑150 (MQ=255)
                                                                                      gCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTACTCCAACGCCATGTGCAAAAAAGCCTTAGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                 >  1:153401/1‑150 (MQ=255)
                                                                                      gCAAGTCCAGCCGCTGAGCTGTGGCAGCGGCTCCCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCCACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGGAGGGTAAGCGGGAGGCGCCGGGAGCGCCTTATcc                                                                 >  2:6562786/1‑150 (MQ=17)
                                                                                       cAATTCCAGCCGCTGATCTGTGTCAGCGGCGACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                <  2:3389883/150‑1 (MQ=255)
                                                                                        aaTTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCAGTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                               <  1:2113681/150‑1 (MQ=255)
                                                                                             cAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                          <  1:10002069/150‑1 (MQ=255)
                                                                                              aGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                                         >  1:5122872/1‑150 (MQ=255)
                                                                                               gCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGt                                                        >  2:3773872/1‑150 (MQ=255)
                                                                                                ccGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa                                                       <  1:4685120/150‑1 (MQ=255)
                                                                                                ccGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCCACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa                                                       <  2:2748928/150‑1 (MQ=255)
                                                                                                 cGCTGATCTGTGTCCGCGGCTACCGAGATTCATTCCCGCCAACAACCGCGCAGTCCGCCAACGCCATGTGCAAAAATGGCTTCGCCTCTGCTTT‑ATTCCTGGCGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTag                                                      <  2:2315682/150‑1 (MQ=21)
                                                                                                 cGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTag                                                      <  1:10869396/150‑1 (MQ=255)
                                                                                                  gCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCTACGCCATGTGCAAAAATGCCTTCGCTGCGGCTGT‑CTGCCTGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTaga                                                     <  2:792483/150‑1 (MQ=255)
                                                                                                   cTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTagag                                                    >  1:10969689/1‑150 (MQ=255)
                                                                                                   cTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTagag                                                    <  2:6061843/150‑1 (MQ=255)
                                                                                                     gATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATCCCCGATCCTGCCTCAAAGCCATATCCGTCCTACGTAGACCa                                                  >  1:7011848/1‑150 (MQ=18)
                                                                                                       tCTGGGTCAGCGGCTACCGTGAGTCCTTCCCGCCAACACCCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                                                <  2:7678900/150‑1 (MQ=255)
                                                                                                          gtgcCAGCGGCTACCGTGATTCAGTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGaa                                             <  1:5214003/150‑1 (MQ=255)
                                                                                                                gCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTa                                       <  1:2559822/150‑1 (MQ=255)
                                                                                                                   gCTACCGTGATTCAGTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGc                                    <  1:1132869/150‑1 (MQ=255)
                                                                                                                     tACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAATATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCt                                  <  2:5883241/150‑1 (MQ=255)
                                                                                                                     tACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAAAAAGCCTTCGCAGCGGCGGT‑CCGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCt                                  <  2:1207751/150‑1 (MQ=255)
                                                                                                                     gaCCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCt                                  <  1:9761600/149‑1 (MQ=255)
                                                                                                                       ccGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCCACGCCATGTGCAACAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGa                                <  1:4481960/150‑1 (MQ=255)
                                                                                                                        cGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGAt                               >  2:5389360/1‑150 (MQ=255)
                                                                                                                        cGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATACGGCCTACGTAGAGCCCTGAACTCGTAGGCCTGAt                               >  1:5005846/1‑150 (MQ=255)
                                                                                                                          tGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATaa                             >  1:1465453/1‑150 (MQ=255)
                                                                                                                          tGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTAATCCCGCAAGCGGCGTGAACGCCATAGCCGGCCGACGTCGAGCACTGAACTCCTAGGCCACATaa                             >  2:10106712/1‑150 (MQ=21)
                                                                                                                           caaacattcCCGCCAACAACCGCGCATTCCTCCAACGCCATCTGCAAAAATGCCTTCCCAGCGGCTGT‑CAGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAg                            <  2:7375050/146‑1 (MQ=255)
                                                                                                                            cttccatcCCGACAACACCCGCGCATGCGCCCAACGCCATGTGCAAAAACGCGTTCGCAGCGGGGGT‑CGGCCTGCAGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGc                           <  2:874048/144‑1 (MQ=18)
                                                                                                                               cattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTa                        <  1:8744428/150‑1 (MQ=255)
                                                                                                                                 ttcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGGCTTCGCATGGGCTTT‑CTGCCTGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgc                      <  2:4132207/150‑1 (MQ=25)
                                                                                                                                 ttcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgc                      <  2:977487/150‑1 (MQ=255)
                                                                                                                                   cccGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgcgc                    <  1:7570072/150‑1 (MQ=255)
                                                                                                                                     cGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATg                                                                                    >  2:6785359/1‑84 (MQ=255)
                                                                                                                                     cGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATg                                                                                    <  1:6785359/84‑1 (MQ=255)
                                                                                                                                       ccATCAACCGCGCAGCCCCCACACGCCATGTGCAAAAATGGCTTCGCCGCGGGTTTCCCGCGCGC‑GTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCa                <  2:1391867/150‑1 (MQ=14)
                                                                                                                                        caacaaCCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTATTTTATGACGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGACCTCGTAGGCCTGATCAGCGTGGCGCATCAg               >  2:3164044/1‑150 (MQ=25)
                                                                                                                                        caacaaCCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAg               <  1:3622471/150‑1 (MQ=255)
                                                                                                                                         aacaacCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAgg              <  1:2423347/150‑1 (MQ=255)
                                                                                                                                            aacCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTGCGTAGAGCACTGACCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa           >  2:9390206/1‑150 (MQ=255)
                                                                                                                                            aacCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa           >  2:5035111/1‑150 (MQ=255)
                                                                                                                                             acCGCGCAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAAt          <  2:9632988/150‑1 (MQ=255)
                                                                                                                                                 cgcATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa      >  1:295301/1‑150 (MQ=255)
                                                                                                                                                 cgcAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa      <  2:726944/150‑1 (MQ=255)
                                                                                                                                                    aTTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTAGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGaa                                             <  2:2838338/108‑1 (MQ=255)
                                                                                                                                                     ttCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCg  >  1:7356915/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGT‑CTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCG  >  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: ‑