Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,505 T→G 5.9% intergenic (+51/+286) 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,5050TG 5.9% 351.0 / 9.1 115intergenic (+51/+286)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (53/54);  new base G (4/3);  total (57/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

Only 100 of 118 total aligned reads displayed.

CGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCT  >  NC_000913/357356‑357652
                                                                                                                                                     |                                                                                                                                                   
cGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                     <  2:1064585/150‑1 (MQ=255)
cGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTGGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                     <  2:1521955/150‑1 (MQ=255)
 gTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                    >  2:356879/1‑150 (MQ=255)
  ttATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                   >  1:33270/1‑150 (MQ=255)
     ttCGGTACGTGGCGGCAAGGGGATTGCCAGCACACAACCGGCACAAACCACAGTATATCTGGAGCATCCAGAAGCCACCGATTACAAACGTTGAACGTCTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCATATcc                                                                                                                                                <  2:147581/150‑1 (MQ=255)
         gTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                            <  2:1120981/150‑2 (MQ=255)
              tGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGCGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGt                                                                                                                                       >  1:818289/1‑150 (MQ=255)
              tGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATaaac                                                                                                                                       <  1:1101469/150‑4 (MQ=38)
                  ggcAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTagag                                                                                                                                   <  2:711058/150‑1 (MQ=255)
                    cAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCa                                                                                                                                 >  1:684084/1‑150 (MQ=255)
                      aGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                                                                                                                               <  1:1063711/150‑1 (MQ=255)
                      aGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                                                                                                                               >  1:994051/1‑150 (MQ=255)
                      aGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                                                                                                                               >  1:741756/1‑150 (MQ=255)
                       ggTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTg                                                                                                                              <  1:1307002/150‑1 (MQ=255)
                       ggTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGACATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTg                                                                                                                              <  1:398762/150‑1 (MQ=255)
                        gTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGa                                                                                                                             >  2:1402819/1‑150 (MQ=255)
                         tGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGaa                                                                                                                            >  1:795626/1‑150 (MQ=255)
                          taTTTCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAc                                                                                                                           <  1:286407/149‑1 (MQ=255)
                          gtgtGCCATCACACAACCGGCACAAACCCCCGTATATCTGGAGCAGCCAGTAGCCATCGATTACACCCGTCGAACCACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGATCACTGAAc                                                                                                                           <  1:123752/147‑1 (MQ=37)
                          gATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCTGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAc                                                                                                                           <  1:991043/150‑1 (MQ=255)
                          gATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAc                                                                                                                           <  1:356879/150‑1 (MQ=255)
                           aTTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACt                                                                                                                          <  1:1237036/150‑1 (MQ=255)
                             tGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCg                                                                                                                        <  2:1184493/150‑1 (MQ=255)
                              gCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGt                                                                                                                       >  2:1201228/1‑150 (MQ=255)
                              gCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGt                                                                                                                       >  2:597295/1‑150 (MQ=255)
                                cAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAg                                                                                                                     >  2:363209/1‑150 (MQ=255)
                                cAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAg                                                                                                                     <  1:922332/150‑1 (MQ=255)
                                    acacaACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGCGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACAGAACTCGTAGGCCt                                                                                                                 >  1:1605867/1‑150 (MQ=255)
                                       caACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGAt                                                                                                              >  2:1595984/1‑150 (MQ=255)
                                        aaCCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATa                                                                                                             <  1:942911/150‑1 (MQ=255)
                                        aaCCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATa                                                                                                             <  1:1472455/150‑1 (MQ=255)
                                         aCCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATaa                                                                                                            >  1:1061562/1‑150 (MQ=255)
                                             gCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGt                                                                                                        >  2:47933/1‑150 (MQ=255)
                                               acaAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAg                                                                                                      <  1:613260/150‑1 (MQ=255)
                                                caAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgc                                                                                                     >  1:1670669/1‑150 (MQ=255)
                                                  aaCCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgcgc                                                                                                   >  1:1403652/1‑150 (MQ=255)
                                                  aaCCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgcgc                                                                                                   <  2:1076960/150‑1 (MQ=255)
                                                    ccaccGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCAt                                                                                                 <  2:907466/150‑1 (MQ=255)
                                                      accGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCa                                                                                               <  1:1402819/150‑1 (MQ=255)
                                                          tataTCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                           <  2:33270/150‑1 (MQ=255)
                                                          tataTCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                           >  2:1704618/1‑150 (MQ=255)
                                                            tatCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAAt                                                                                         >  1:1576500/1‑150 (MQ=255)
                                                             atCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAAtt                                                                                        >  2:8562/1‑150 (MQ=255)
                                                             atCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAAtt                                                                                        >  1:833935/1‑150 (MQ=255)
                                                              tCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTc                                                                                       >  2:314532/1‑150 (MQ=255)
                                                               cTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTcc                                                                                      <  1:1201228/150‑1 (MQ=255)
                                                               cTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTcc                                                                                      <  2:121607/150‑1 (MQ=255)
                                                                tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACTCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                                     >  1:580997/1‑150 (MQ=255)
                                                                tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                                     >  1:674817/1‑150 (MQ=255)
                                                                tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                                     >  2:504/1‑150 (MQ=255)
                                                                   cgcagcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGcc                                                                                  <  1:524172/149‑1 (MQ=255)
                                                                      agcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCt                                                                               >  1:1302704/1‑150 (MQ=255)
                                                                       gcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTg                                                                              >  1:7311/1‑150 (MQ=255)
                                                                       gcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTg                                                                              >  1:1615707/1‑150 (MQ=255)
                                                                       gcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCATGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTg                                                                              >  1:523917/1‑150 (MQ=255)
                                                                         cAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGAt                                                                            >  2:794780/1‑150 (MQ=255)
                                                                            aaGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCtg                                                                         <  2:693839/150‑1 (MQ=255)
                                                                             aGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCtgt                                                                        >  1:1137050/1‑150 (MQ=255)
                                                                              gCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCtgtg                                                                       <  2:488746/150‑1 (MQ=255)
                                                                                  tCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAg                                                                   <  2:707564/150‑1 (MQ=255)
                                                                                   cGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGATTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCAATGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCAGCAGCGCTGTGTCAGc                                                                  >  1:71145/1‑150 (MQ=17)
                                                                                    gATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCg                                                                 <  1:1595984/150‑1 (MQ=255)
                                                                                    gATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCg                                                                 <  2:795626/150‑1 (MQ=255)
                                                                                    gATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCg                                                                 >  1:64477/1‑150 (MQ=255)
                                                                                      ttACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGGCAGCGGc                                                               <  1:641551/150‑1 (MQ=255)
                                                                                       tACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGATCTGTGTCAGCGGCt                                                              >  2:468452/1‑150 (MQ=255)
                                                                                         cAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTAc                                                            <  2:1576500/150‑1 (MQ=255)
                                                                                         cAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTAc                                                            >  1:1482543/1‑150 (MQ=255)
                                                                                         cAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTAc                                                            >  1:838566/1‑150 (MQ=255)
                                                                                           aaCGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCg                                                          >  1:491307/1‑150 (MQ=255)
                                                                                            aCGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGt                                                         <  1:597295/150‑1 (MQ=255)
                                                                                            aCGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGt                                                         <  1:281964/150‑1 (MQ=255)
                                                                                                tGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGctt                                                     >  2:533018/1‑150 (MQ=255)
                                                                                                 gAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattc                                                    <  2:574923/150‑1 (MQ=255)
                                                                                                 gAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattc                                                    >  1:1586345/1‑150 (MQ=255)
                                                                                                 gAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattc                                                    <  2:442796/150‑1 (MQ=255)
                                                                                                    cGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattcatt                                                 >  1:800401/1‑150 (MQ=255)
                                                                                                       cTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTccc                                              <  1:1331339/150‑1 (MQ=255)
                                                                                                       cTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTccc                                              >  1:415970/1‑150 (MQ=255)
                                                                                                        tGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCg                                             <  2:610557/150‑1 (MQ=255)
                                                                                                         gggTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTGATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGc                                            <  2:1061562/150‑1 (MQ=255)
                                                                                                                aGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCcaacacc                                     >  1:1442150/1‑148 (MQ=255)
                                                                                                                 gCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGAACTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAAcc                                    >  2:590992/1‑150 (MQ=255)
                                                                                                                     gCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCgcgc                                <  2:86631/150‑1 (MQ=39)
                                                                                                                          gTTTATGCCGGATGCGGCGTGAACGCCTTATACGGACTACGTAGAGCATTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATtcc                           <  2:833935/150‑1 (MQ=255)
                                                                                                                           tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTACAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATtcct                          <  1:129399/150‑1 (MQ=255)
                                                                                                                           tttATGCCGGATGCGGCGTGAACGACTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATtcct                          <  2:1005689/150‑1 (MQ=255)
                                                                                                                             tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATtcctcc                        <  2:7311/150‑1 (MQ=40)
                                                                                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCaa                      >  2:1617923/2‑150 (MQ=14)
                                                                                                                                gCCGGATGCGGCGTGAACGCCTCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACACCCGCGCATTCCTCCAAc                     >  1:748528/1‑150 (MQ=255)
                                                                                                                                gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc                     >  2:402092/1‑150 (MQ=14)
                                                                                                                                gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc                     >  1:967617/1‑150 (MQ=14)
                                                                                                                                 ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg                    >  2:1227500/1‑150 (MQ=255)
                                                                                                                                 ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg                    >  2:84719/1‑150 (MQ=255)
                                                                                                                                 ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg                    <  1:1054316/150‑1 (MQ=255)
                                                                                                                                 ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg                    <  2:1615707/150‑1 (MQ=255)
                                                                                                                                        cggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCAtgtg             >  1:52622/1‑150 (MQ=255)
                                                                                                                                          gcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCa           <  2:205356/150‑1 (MQ=255)
                                                                                                                                              tcaCGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCaaaaa       <  1:1617923/148‑1 (MQ=255)
                                                                                                                                                   ccTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCt  <  2:1243804/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCT  >  NC_000913/357356‑357652

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