Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,606 C→T 5.1% intergenic (+152/+185) 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,6060CT 5.1% 868.1 / 22.8 264intergenic (+152/+185)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (80/162);  new base T (11/2);  total (98/166)
Fisher's exact test for biased strand distribution p-value = 3.03e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

GACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGC  >  NC_000913/357457‑357753
                                                                                                                                                     |                                                                                                                                                    
  cTGGTTTACAGCTTGATTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGAACTGTGTCAGCGGCTACCGTGATTCATTtcc                                                                                                                                                    <  2:2583461/150‑1 (MQ=11)
  cTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTccc                                                                                                                                                    <  1:2459587/150‑1 (MQ=255)
  cTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTccc                                                                                                                                                    <  2:1818319/150‑1 (MQ=255)
   tGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCg                                                                                                                                                   >  1:1058719/1‑150 (MQ=255)
   tGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCg                                                                                                                                                   >  1:4123583/1‑150 (MQ=255)
   tGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTACAGCCGCAGAACTGTGTCAGCGGCGACCGTGATTAATTaccc                                                                                                                                                   >  1:802986/1‑149 (MQ=14)
   tGGGTTACAGCGAGCTTAGTATATGCCGGATGCGGCGTGAACGACTTATCCGGCCTACGAAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCg                                                                                                                                                   >  1:1251072/1‑150 (MQ=11)
     ggTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGcc                                                                                                                                                 >  1:355021/1‑150 (MQ=255)
     ggTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGcc                                                                                                                                                 >  1:3611440/1‑150 (MQ=255)
     ggTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTACCGcc                                                                                                                                                 >  1:528252/1‑150 (MQ=255)
     ggTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCGGATCTGTGGCAGCGGCTAGCGTGATTCATTCCCGcc                                                                                                                                                 >  1:859403/1‑150 (MQ=255)
          cAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCcaacaa                                                                                                                                            <  1:4047851/150‑1 (MQ=255)
                 cTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCa                                                                                                                                     <  1:1582162/150‑1 (MQ=38)
                           gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc                                                                                                                           <  1:1296611/150‑1 (MQ=255)
                             cGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGc                                                                                                                         >  2:4008730/1‑150 (MQ=255)
                                  gcggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCAtgt                                                                                                                    <  2:616677/150‑1 (MQ=255)
                                   cggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTCCTCCAACGCCAtgtg                                                                                                                   >  2:3035161/1‑150 (MQ=37)
                                   cggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGAGCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCAtgtg                                                                                                                   >  2:4069967/1‑150 (MQ=21)
                                     gcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCa                                                                                                                 <  1:3285607/150‑1 (MQ=255)
                                          aaCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTtcc                                                                                                                                                    >  1:3988753/1‑110 (MQ=17)
                                          aaCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCTACAACCGCGCATTCCTCCAACGCCATGTGCAAAAAt                                                                                                            <  2:3158099/150‑1 (MQ=255)
                                          aaCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAAt                                                                                                            <  1:51182/150‑1 (MQ=255)
                                             gCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGcc                                                                                                         >  1:2335663/1‑150 (MQ=255)
                                               ctgaTCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCtt                                                                                                       <  1:1440340/147‑1 (MQ=255)
                                               cTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCtt                                                                                                       <  1:2219115/150‑1 (MQ=255)
                                                ttATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTc                                                                                                      <  1:193446/150‑1 (MQ=255)
                                                ttATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTc                                                                                                      >  2:2406797/1‑150 (MQ=255)
                                                 tATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCg                                                                                                     <  2:459984/150‑1 (MQ=255)
                                                      ggCCTACGTAGAGCACTGAACTCGTAGGACTCATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCCTGATTCAATCCCCCCAACAACCGCGCAATCCTCCAACGCCATGAGCAAAAATTCCTTCGCAGCg                                                                                                >  1:3124387/1‑150 (MQ=14)
                                                        ccTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGc                                                                                              >  2:3616484/1‑150 (MQ=255)
                                                          tACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTg                                                                                            >  2:1890799/1‑150 (MQ=255)
                                                              tAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTAATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTg                                                                                       >  1:3659837/1‑150 (MQ=12)
                                                               agagCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAACCGCGCATTGCTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGc                                                                                      <  1:302301/150‑1 (MQ=255)
                                                                   cACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGc                                                                                  >  2:3687898/1‑150 (MQ=255)
                                                                    aCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAAGTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCt                                                                                 <  1:1402810/150‑1 (MQ=255)
                                                                      tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCcaa                                                                                                                                               >  1:604452/1‑87 (MQ=255)
                                                                      tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAg                                                                                   >  1:2775809/1‑146 (MQ=255)
                                                                      tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGt                                                                               >  1:3240244/1‑150 (MQ=255)
                                                                      tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAACAGCGCATTCCTCAAACGCCATGTGCAAAAATGCCGTCGCTGCGGCTGTCTGCCAGCTGt                                                                               >  1:3675622/1‑150 (MQ=14)
                                                                        aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAg                                                                             <  1:169437/150‑1 (MQ=255)
                                                                        aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAg                                                                             >  2:2138674/1‑150 (MQ=255)
                                                                         aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGt                                                                            >  2:3601689/1‑150 (MQ=255)
                                                                           tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGttt                                                                          <  1:779688/150‑1 (MQ=255)
                                                                           tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGttt                                                                          >  1:782233/1‑150 (MQ=255)
                                                                           tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCATCCAACGCCATGTGCAAAAATGCCTTCGCTGCGGCTGTCTGCCAGCtatttttc                                                                          <  1:2691213/150‑8 (MQ=255)
                                                                             gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTAt                                                                        <  2:709365/150‑1 (MQ=255)
                                                                              tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATg                                                                       <  1:2828605/150‑1 (MQ=255)
                                                                                 gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCg                                                                    <  1:1662530/150‑1 (MQ=255)
                                                                                  ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGc                                                                                                 >  2:3483401/1‑121 (MQ=255)
                                                                                  ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGc                                                                                                 <  1:1426712/121‑1 (MQ=255)
                                                                                  ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGc                                                                                                 <  1:3483401/121‑1 (MQ=255)
                                                                                   cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGa                                                                  <  2:2776410/150‑1 (MQ=255)
                                                                                       tAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcg                                                              <  2:3719414/150‑1 (MQ=255)
                                                                                        tagagTATCGCATCATGCAATTCCATCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcgg                                                             <  2:2228730/146‑1 (MQ=255)
                                                                                        aaGCGTAGCGCATCAGGCAATTCCAGCGGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCGGTAGTTTATGCCGGATgcgg                                                             <  2:2997274/150‑1 (MQ=255)
                                                                                        aaGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcgg                                                             <  1:3464252/150‑1 (MQ=255)
                                                                                        aaGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGAGTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcgg                                                             <  1:1507031/150‑1 (MQ=255)
                                                                                           cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGTCGGATGCGGCGt                                                          <  2:1145352/150‑1 (MQ=255)
                                                                                               gcgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTATATTTTATCCCCGAAGCAGCCTGaaa                                                      >  2:2351889/1‑149 (MQ=14)
                                                                                                  cATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGcc                                                   <  1:187237/150‑1 (MQ=255)
                                                                                                        gCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCa                                                                                    <  2:1982236/111‑1 (MQ=255)
                                                                                                        gCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGAAGCGGCGTGAACGCCTTATcc                                             >  1:1965682/1‑150 (MQ=255)
                                                                                                         cAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCg                                            <  1:2363019/150‑1 (MQ=255)
                                                                                                         cAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCg                                            <  1:1204458/150‑1 (MQ=255)
                                                                                                          aaTTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCAGTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCgg                                           <  2:3672824/150‑1 (MQ=255)
                                                                                                           aTTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                          <  1:103708/150‑1 (MQ=255)
                                                                                                           aTTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                          <  2:1239352/150‑1 (MQ=255)
                                                                                                            ttCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCACCGCCATGTGCAATAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                         <  1:690641/150‑1 (MQ=255)
                                                                                                            ttCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                         <  1:397357/150‑1 (MQ=255)
                                                                                                             tCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTACTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                        >  1:2455829/1‑150 (MQ=255)
                                                                                                              ccAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                       >  2:274249/1‑150 (MQ=255)
                                                                                                               cAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCCGAAGATTCCTCAAAGCATTATCCTGCCTAc                                      >  2:2527833/1‑150 (MQ=17)
                                                                                                                 gCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGt                                    <  2:355021/150‑1 (MQ=255)
                                                                                                                 gCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGt                                    <  2:24716/150‑1 (MQ=255)
                                                                                                                 gCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGt                                    <  1:3762059/150‑1 (MQ=255)
                                                                                                                  ccGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCACAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa                                   >  1:2440478/1‑150 (MQ=255)
                                                                                                                  ccGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa                                   >  1:1105703/1‑150 (MQ=255)
                                                                                                                    gCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTCTATTTGATCCCCGAGGCGGCCTCAACGCCATAGCCGGCCTACGTaga                                 >  2:2303826/1‑150 (MQ=18)
                                                                                                                     cTGATCTGTGTCAGCGGCTACCGTGATTCAGTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTGCGTagag                                <  2:1314821/150‑1 (MQ=255)
                                                                                                                        aTCTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAACCGCGCATTACTCCAACGCCATGTGCAAAAATGAATACCCACCTGCAGTCAGACAGCTTAACTAAATGACTGATTCTGCGTGAACGCCTTATCCGGCATACGTAGAGctc                             >  2:548674/1‑148 (MQ=12)
                                                                                                                         tCTGTGTCAGCGGCTACCGTGAGTCATTCCCGCCACCAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                            <  2:3826280/150‑1 (MQ=255)
                                                                                                                         tCTGTGTCAGCGGCTACCGTGACTCATTCCCGCCCACATCCGCGGAGTACTCCCACGCCATGTGAAAAAATGCCTTCGCATCGGCGGTCCTGCTTGTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                            <  2:3395302/150‑1 (MQ=21)
                                                                                                                               tCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTc                      <  1:336743/150‑1 (MQ=255)
                                                                                                                               tCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTc                      <  2:1687088/150‑1 (MQ=255)
                                                                                                                                   cGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAg                  <  2:2274438/150‑1 (MQ=255)
                                                                                                                                     gCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGc                <  2:3587428/150‑1 (MQ=255)
                                                                                                                                     gCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCGTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGc                >  2:3756767/1‑150 (MQ=255)
                                                                                                                                       tACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCt              >  1:3584161/1‑150 (MQ=255)
                                                                                                                                       tACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGAGGCGGCGTGACCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCt              >  2:2628741/1‑150 (MQ=255)
                                                                                                                                       tACCGTGATTCATTCCCGCCAACAACCGCGCATGCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCt              <  2:1945401/150‑1 (MQ=255)
                                                                                                                                          cGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGAAGCGGCGTGAACGCCTTATCCGGCCTTCGTAGAGCACGGAACTCGTAGGCCTGAt           >  2:892432/1‑150 (MQ=255)
                                                                                                                                          cGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATCCCGGATGCGGCGTCAACGCATTAACCGGGCTACGTGGAGCACTGAAGTCGTAGGCCGGAt           >  2:798620/1‑150 (MQ=25)
                                                                                                                                          cGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTCTATTTTATGCCGGATGCGGCGTGACCGCCCTCGCCGGCCGTCGTAGAGCAATGAACTCGTAggacagag           >  2:1198011/1‑144 (MQ=21)
                                                                                                                                          cGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTCTAGTTTATCCCGGAAGCGACCTCAACGCCATCTCCGGCCTACGTTGATCACTGAAAGCGAAGGccacat           >  2:3394325/1‑146 (MQ=18)
                                                                                                                                              attcattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCTGCGGCTGTCTGCCAGCTGTAGTTGATCCCGGAGGCGCCGTGACCGCGGTGTCCGGCCTACGTAGAGCACTGAACTCGTAGGCGGCATAAGc       >  1:2876810/1‑150 (MQ=21)
                                                                                                                                              attcattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTGATGCCGGATGCGGCGTGAACGCATTATCCGGCCTGCGGAGAGCACTGAAATCGTAGGCCTGATAAGc       >  1:256253/1‑150 (MQ=37)
                                                                                                                                               ttcattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCg      >  1:2158322/1‑150 (MQ=255)
                                                                                                                                                 cattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTa    <  1:2406797/150‑1 (MQ=255)
                                                                                                                                                 cattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTa    <  1:2288311/150‑1 (MQ=255)
                                                                                                                                                   ttcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCg                     <  1:2467024/131‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGC  >  NC_000913/357457‑357753

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: