Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,506 T→C 7.0% intergenic (+52/+285) 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,5060TC 7.0% 335.9 / 12.3 115intergenic (+52/+285)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (52/54);  new base C (5/3);  total (57/58)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01

Only 100 of 116 total aligned reads displayed.

GTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCT  >  NC_000913/357357‑357652
                                                                                                                                                     |                                                                                                                                                  
gTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                    >  2:356879/1‑150 (MQ=255)
 ttATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                   >  1:33270/1‑150 (MQ=255)
    ttCGGTACGTGGCGGCAAGGGGATTGCCAGCACACAACCGGCACAAACCACAGTATATCTGGAGCATCCAGAAGCCACCGATTACAAACGTTGAACGTCTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCATATcc                                                                                                                                                <  2:147581/150‑1 (MQ=255)
        gTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                            <  2:179572/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)
              ggcggcAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa                                                                                                                                      >  1:530817/1‑150 (MQ=255)
                 ggcAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTagag                                                                                                                                   <  2:711058/150‑1 (MQ=255)
                   cAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCa                                                                                                                                 >  1:684084/1‑150 (MQ=255)
                     aGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                                                                                                                               >  1:994051/1‑150 (MQ=255)
                     aGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                                                                                                                               <  1:1063711/150‑1 (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)
                         gATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAc                                                                                                                           <  2:818289/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)
                              ccAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTa                                                                                                                      >  1:1483590/1‑150 (MQ=255)
                               cAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAg                                                                                                                     <  1:922332/150‑1 (MQ=255)
                               cAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAg                                                                                                                     >  2:363209/1‑150 (MQ=255)
                                   acacaACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGCGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACAGAACTCGTAGGCCt                                                                                                                 >  1:1605867/1‑150 (MQ=255)
                                       aaCCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATa                                                                                                             <  1:1472455/150‑1 (MQ=255)
                                       aaCCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATa                                                                                                             <  1:942911/150‑1 (MQ=255)
                                        aCCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATaa                                                                                                            >  1:1061562/1‑150 (MQ=255)
                                        aCCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATaa                                                                                                            <  1:363209/150‑1 (MQ=255)
                                            gCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGt                                                                                                        >  2:47933/1‑150 (MQ=255)
                                              acaAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAg                                                                                                      <  1:613260/150‑1 (MQ=255)
                                               caAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgc                                                                                                     >  1:1670669/1‑150 (MQ=255)
                                                 aaCCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgcgc                                                                                                   <  2:1076960/150‑1 (MQ=255)
                                                 aaCCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgcgc                                                                                                   >  1:1403652/1‑150 (MQ=255)
                                                   ccaccGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCAt                                                                                                 <  2:907466/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)
                                                             tCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTc                                                                                       <  2:684084/150‑1 (MQ=255)
                                                              cTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTcc                                                                                      <  1:1201228/150‑1 (MQ=255)
                                                               tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACTCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                                     >  1:580997/1‑150 (MQ=255)
                                                               tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                                     >  1:674817/1‑150 (MQ=255)
                                                               tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                                     >  1:1371571/1‑150 (MQ=255)
                                                               tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                                     >  2:504/1‑150 (MQ=255)
                                                                  cgcagcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGcc                                                                                  <  1:524172/149‑1 (MQ=255)
                                                                     cgcCAGAAGCCTTCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGCTTATGCCGGCTGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCt                                                                               <  1:666335/149‑1 (MQ=255)
                                                                     agcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCt                                                                               >  1:1302704/1‑150 (MQ=255)
                                                                      gcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTg                                                                              >  1:1615707/1‑150 (MQ=255)
                                                                      gcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTg                                                                              >  1:7311/1‑150 (MQ=255)
                                                                        cAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGAt                                                                            >  1:86631/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)
                                                                                 tCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAg                                                                   <  2:707564/150‑1 (MQ=255)
                                                                                  cGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGATTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCAATGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCAGCAGCGCTGTGTCAGc                                                                  >  1:71145/1‑150 (MQ=17)
                                                                                   gATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCg                                                                 >  1:64477/1‑150 (MQ=255)
                                                                                   gATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCg                                                                 <  1:1595984/150‑1 (MQ=255)
                                                                                      tACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGATCTGTGTCAGCGGCt                                                              >  2:468452/1‑150 (MQ=255)
                                                                                       aCAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTa                                                             <  1:1641946/150‑1 (MQ=255)
                                                                                        cAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTAc                                                            >  1:1482543/1‑150 (MQ=255)
                                                                                        cAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTAc                                                            <  2:1576500/150‑1 (MQ=255)
                                                                                          aaCGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCg                                                          >  1:491307/1‑150 (MQ=255)
                                                                                           aCGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGt                                                         <  1:597295/150‑1 (MQ=255)
                                                                                              ttGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCATGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGat                                                      <  2:523917/150‑1 (MQ=255)
                                                                                               tGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGctt                                                     >  2:533018/1‑150 (MQ=255)
                                                                                                gAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattc                                                    <  2:574923/150‑1 (MQ=255)
                                                                                                gAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattc                                                    <  2:442796/150‑1 (MQ=255)
                                                                                                  aCGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattcat                                                  <  2:1483590/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)
                                                                                                        gggTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTGATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGc                                            <  2:1061562/150‑1 (MQ=255)
                                                                                                               aGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCcaacacc                                     >  1:1442150/1‑148 (MQ=255)
                                                                                                                gCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGAACTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAAcc                                    >  2:590992/1‑150 (MQ=255)
                                                                                                                 cGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACcg                                   >  2:1204096/1‑150 (MQ=255)
                                                                                                                    gCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCgcgc                                <  2:86631/150‑1 (MQ=39)
                                                                                                                    gCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGGGGCAGCGGCTACCGGGATTCCGTCACGCCAACAACCgcgc                                >  1:1206045/1‑150 (MQ=14)
                                                                                                                         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                     >  1:967617/1‑150 (MQ=14)
                                                                                                                               gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc                     >  2:402092/1‑150 (MQ=14)
                                                                                                                               gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc                     >  2:395943/1‑150 (MQ=14)
                                                                                                                                ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg                    <  1:1054316/150‑1 (MQ=255)
                                                                                                                                ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg                    >  2:84719/1‑150 (MQ=255)
                                                                                                                                ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg                    <  2:1615707/150‑1 (MQ=255)
                                                                                                                                ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg                    >  2:1227500/1‑150 (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)
                                                                                                                                                   cTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattca                                                   >  1:788072/1‑100 (MQ=37)
                                                                                                                                                     |                                                                                                                                                  
GTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCT  >  NC_000913/357357‑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: