Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,778 C→T 7.3% intergenic (+324/+13) 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,7780CT 7.3% 1168.5 / 51.4 374intergenic (+324/+13)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (215/126);  new base T (16/11);  total (235/139)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 377 total aligned reads displayed.

CCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG  >  NC_000913/357629‑357927
                                                                                                                                                     |                                                                                                                                                     
   aCGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtg                                                                                                                                                    >  2:5311892/1‑150 (MQ=255)
    cGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtgt                                                                                                                                                   <  2:1017273/150‑1 (MQ=255)
         tgtgCAAAATTGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGCGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCAACAGGCAATTCCAGCCGCAGACCTGTGTCAg                                                                                                                                              <  1:1556566/150‑1 (MQ=255)
           tgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCg                                                                                                                                            >  1:5015027/1‑150 (MQ=255)
             cAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGc                                                                                                                                          >  1:80806/1‑150 (MQ=255)
              aaaaaTGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCt                                                                                                                                         >  2:1886610/1‑150 (MQ=255)
              aaaaaTGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCt                                                                                                                                         >  2:4840265/1‑150 (MQ=255)
              aaaaaTGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCt                                                                                                                                         >  1:334981/1‑150 (MQ=255)
              aaaaaTGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCt                                                                                                                                         >  1:2943476/1‑150 (MQ=255)
              aaaaaTGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATACGGCCTACGTGGAGCACTGAAATCGTAGGCCTGATAATCGTAGCGCCTCAGGCTATTCCAGCCGCAGACCTGTGTCTGAGGCt                                                                                                                                         >  2:1167728/1‑150 (MQ=18)
               aaaaTGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTa                                                                                                                                        >  1:3839453/1‑150 (MQ=255)
               aaaaTGCCTTCGCAGCGGCTGTCGGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACCCGTTGGCCTGTTTTGGGTAGCGCATCTGGCTATTCCAGCCGCAGACCTGTGTCAGCGGCTa                                                                                                                                        <  1:2615698/150‑1 (MQ=18)
                aactGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTAc                                                                                                                                       <  2:2021125/147‑1 (MQ=255)
                aaaTGCCTTCGCAGCGGCGGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTAc                                                                                                                                       <  1:3015932/150‑1 (MQ=255)
                 aaTGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTAcc                                                                                                                                      >  1:3228859/1‑150 (MQ=255)
                     cgttCGCCGCGGCGGTCTGGCTGCTTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGa                                                                                                                                  <  2:3559262/148‑1 (MQ=14)
                     ccTTCGCAGCGGCTGTCTGCCTGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGa                                                                                                                                  <  2:93421/150‑1 (MQ=255)
                      cTTCTCAGCGGGTGTATGCCTGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAt                                                                                                                                 <  2:574441/150‑1 (MQ=255)
                       ttCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAtt                                                                                                                                >  2:4910368/1‑150 (MQ=255)
                          gCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCAt                                                                                                                             >  2:5096440/1‑150 (MQ=255)
                          gCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCAt                                                                                                                             >  2:2279498/1‑150 (MQ=255)
                            aGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCAttt                                                                                                                           >  2:3362467/1‑150 (MQ=255)
                             gCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACAGAACTCGTATGCCTGATAAGCGTAGCGCATCATGCAAATCCAGCCGCAGACCTGTGACAGCGGCTACCGTGATTCATTTc                                                                                                                          >  1:5242965/1‑150 (MQ=17)
                              cGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTcc                                                                                                                         <  1:2529897/150‑1 (MQ=255)
                                gcgggccgCCAGCTGTAGCTGCTGCCGGATGCGGCGTGAACGCCTTATCCGGGCTACGTAGAGCACTGAACTCGTAGGCCTGCTATGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGc                                                                                                                       <  2:4215120/143‑1 (MQ=14)
                                      tGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAAAGCCTTATCCGGCCTAAGTAGAGAACTGAAATGGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTCCAGCCGCAGACCTGTTTAAGCGGCTACCGTGATTCATTTCCGCcaacaa                                                                                                                 >  1:3735542/1‑150 (MQ=14)
                                          aGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcgc                                                                                                             >  2:4314555/1‑150 (MQ=40)
                                             cttAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCa                                                                                                          <  1:1797792/148‑1 (MQ=25)
                                                  tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCCATTCCAGCCGCAGACCTGTGTCAGCGGCAACCGTGATTCATTTCCGCCCACAACCGCGCATTTAt                                                                                                     >  1:3543543/1‑150 (MQ=255)
                                                  tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAt                                                                                                     <  2:4619295/150‑1 (MQ=39)
                                                  tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAt                                                                                                     >  2:1442568/1‑150 (MQ=39)
                                                   ttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATc                                                                                                    <  1:5106742/150‑1 (MQ=40)
                                                      ggCCGGATAAGGCGTTCACGCCGCATCCGGCATACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCaa                                                                                                 >  2:2650815/2‑150 (MQ=14)
                                                         cGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGc                                                                                              >  2:375254/1‑150 (MQ=255)
                                                           gATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCAGTGAACTCGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCa                                                                                            <  1:565992/150‑1 (MQ=255)
                                                           gATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTc                                                                                                                               <  1:107997/115‑1 (MQ=35)
                                                           gATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCa                                                                                            <  1:4194613/150‑1 (MQ=255)
                                                             tGCGGCGTGAACGCCTTCTCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAtg                                                                                          <  2:1636739/150‑1 (MQ=255)
                                                              gcggcgTGAACGCCTTTTCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAGTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAtgt                                                                                         <  1:159402/150‑1 (MQ=255)
                                                              gcggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAtgt                                                                                         <  2:4120400/150‑1 (MQ=255)
                                                                  cgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaa                                                                                     <  1:4111037/150‑1 (MQ=255)
                                                                    tGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaaaa                                                                                   <  2:4709514/150‑1 (MQ=255)
                                                                     gAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaaaaa                                                                                  <  2:728429/150‑1 (MQ=255)
                                                                       aCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTCGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATg                                                                                >  1:2421284/1‑150 (MQ=255)
                                                                       aCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATg                                                                                >  1:3987464/1‑150 (MQ=17)
                                                                        cGCCTTATCGGGCCTACGTAGAGCACTGAAGTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGTTGTGTGTCAGCGGGTAGCGTGATTCTTTTGCGCCTACAACCGCGCATTTATCCAACGCCTTGTGCAAAATTGc                                                                               >  2:3189289/1‑150 (MQ=11)
                                                                        cGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGc                                                                               <  2:1074961/150‑1 (MQ=255)
                                                                         gCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGcc                                                                              >  2:2423577/1‑150 (MQ=255)
                                                                         gCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGcc                                                                              >  2:1703716/1‑150 (MQ=255)
                                                                           cTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCtt                                                                            >  1:2221446/1‑150 (MQ=255)
                                                                              aTCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTcgc                                                                         >  1:1403894/1‑150 (MQ=255)
                                                                               tCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTcgcg                                                                        >  2:3251558/1‑150 (MQ=255)
                                                                                ccGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCgcgg                                                                       >  1:3465504/1‑150 (MQ=255)
                                                                                ccGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCAACAGGCAATTCCAGCCGCAGAACTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCgcgg                                                                       <  1:3189289/150‑1 (MQ=255)
                                                                                  ggCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCgcggcg                                                                     >  2:4372219/1‑150 (MQ=255)
                                                                                   gCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGcggcgg                                                                    <  1:5382848/150‑1 (MQ=255)
                                                                                   gCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCACATCAGGCAAGTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGcggcgg                                                                    <  2:2316609/150‑1 (MQ=255)
                                                                                     cTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCGTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCt                                                                  >  2:2456172/1‑150 (MQ=255)
                                                                                       aCGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGt                                                                <  1:4842990/150‑1 (MQ=255)
                                                                                        cGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTc                                                               >  1:651197/1‑150 (MQ=255)
                                                                                            gagCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGcc                                                           <  2:2995323/150‑1 (MQ=255)
                                                                                             agCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCa                                                          >  1:3040571/1‑150 (MQ=255)
                                                                                             agCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCa                                                          <  1:3437795/150‑1 (MQ=255)
                                                                                              gCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAg                                                         <  1:4750399/150‑1 (MQ=255)
                                                                                               cACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGAGCTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAACCGCGCATTAATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGc                                                        >  1:4668201/1‑150 (MQ=255)
                                                                                                  tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt                                                     >  1:1922141/1‑150 (MQ=255)
                                                                                                  tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt                                                     >  2:817941/1‑150 (MQ=255)
                                                                                                  tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt                                                     >  2:3214260/1‑150 (MQ=255)
                                                                                                       tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                <  1:101785/150‑1 (MQ=255)
                                                                                                       tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                <  1:5185691/150‑1 (MQ=255)
                                                                                                         gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGAGCTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                              >  2:2233081/1‑150 (MQ=11)
                                                                                                         gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGACTGCCAGCTCTGTTTccg                                              >  1:230156/1‑150 (MQ=255)
                                                                                                             gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                          >  2:1865356/1‑150 (MQ=255)
                                                                                                              ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                         >  1:3492532/1‑150 (MQ=255)
                                                                                                              ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                         >  2:3118072/1‑150 (MQ=255)
                                                                                                               cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa                                        >  1:4507378/1‑150 (MQ=255)
                                                                                                                tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                       <  1:4275078/150‑1 (MQ=255)
                                                                                                                  aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGATCAGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAACCGCGCATTAATCCAACGCCATGTGCAAAAATGCCTTCGCCGCGGCTGTCTGCCAGCTAATTTTCATCCCcaacaa                                     >  1:2914615/1‑150 (MQ=11)
                                                                                                                  aTAAGCGTAGCGCAACAGGAAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGAAAAAATGCCTTCGCTGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaccaa                                     >  1:1185031/1‑150 (MQ=17)
                                                                                                                    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                   >  2:4304907/1‑150 (MQ=255)
                                                                                                                      gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTAc                                                                                                                                       <  2:1749480/48‑1 (MQ=255)
                                                                                                                         tAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCTCAACAAAAccgcc                              <  2:1496357/150‑1 (MQ=255)
                                                                                                                           gcgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAACCGCGCATTAATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCAGCCGCAACAAAAccgcaga                            >  2:2731078/1‑147 (MQ=14)
                                                                                                                               aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                        >  2:672403/1‑150 (MQ=255)
                                                                                                                                tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                       <  2:4171661/150‑1 (MQ=255)
                                                                                                                                tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                       >  2:3451562/1‑150 (MQ=255)
                                                                                                                                  aGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                              <  2:5000491/141‑1 (MQ=14)
                                                                                                                                   ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCa                    >  2:3671938/1‑150 (MQ=255)
                                                                                                                                    gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCag                   >  2:4698055/1‑150 (MQ=255)
                                                                                                                                    gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCag                   >  1:3294426/1‑150 (MQ=255)
                                                                                                                                     cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                  >  1:4556201/1‑150 (MQ=255)
                                                                                                                                      aaTTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagta                 <  2:3492532/150‑1 (MQ=255)
                                                                                                                                      aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagta                 <  2:4539357/150‑1 (MQ=255)
                                                                                                                                           cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCCTTACCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAACAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc            <  1:4419233/150‑1 (MQ=255)
                                                                                                                                            aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg           >  2:129173/1‑150 (MQ=255)
                                                                                                                                             gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg          <  1:5315842/150‑1 (MQ=255)
                                                                                                                                             gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg          >  2:4445907/1‑150 (MQ=255)
                                                                                                                                              ccGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg         <  1:707481/150‑1 (MQ=17)
                                                                                                                                              ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTACCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGGAGTGGCggg         >  2:286426/1‑150 (MQ=255)
                                                                                                                                                    aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga   >  2:521866/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG  >  NC_000913/357629‑357927

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: