Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,492 C→A 5.9% intergenic (+38/+299) 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,4920CA 5.9% 366.5 / 9.2 118intergenic (+38/+299)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (54/57);  new base A (4/3);  total (58/60)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

Only 100 of 121 total aligned reads displayed.

CAGGTTCCGGTACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTG  >  NC_000913/357344‑357641
                                                                                                                                                    |                                                                                                                                                     
cAGGTTCCGGTACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAACCGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATgcg                                                                                                                                                      >  2:666335/1‑150 (MQ=255)
  ggTTCCGGTACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATgcggc                                                                                                                                                    >  2:1641946/1‑150 (MQ=255)
          tACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAACCCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                                            <  2:786343/150‑1 (MQ=255)
          tACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                                            >  1:707564/1‑150 (MQ=255)
            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)
                     gTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                 <  2:179572/150‑1 (MQ=255)
                          tGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGCGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGt                                                                                                                            >  1:818289/1‑150 (MQ=255)
                          tGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATaaac                                                                                                                            <  1:1101469/150‑4 (MQ=38)
                           ggcggcAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTg                                                                                                                                                 >  2:488704/1‑128 (MQ=255)
                           ggcggcAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa                                                                                                                           >  1:530817/1‑150 (MQ=255)
                                cAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCa                                                                                                                      >  1:684084/1‑150 (MQ=255)
                                  aGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                                                                                                                    <  1:1063711/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)
                                      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)
                                       aTTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACt                                                                                                               <  1:1237036/150‑1 (MQ=255)
                                         tGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCg                                                                                                             <  2:1184493/150‑1 (MQ=255)
                                          gCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGt                                                                                                            >  2:1201228/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)
                                                    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)
                                                         gCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGt                                                                                             <  2:1022557/150‑1 (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:684084/150‑1 (MQ=255)
                                                                           cTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTcc                                                                           <  2:121607/150‑1 (MQ=255)
                                                                            tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                          >  2:504/1‑150 (MQ=255)
                                                                            tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                          >  1:1371571/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:7311/1‑150 (MQ=255)
                                                                                   gcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTg                                                                   >  1:1615707/1‑150 (MQ=255)
                                                                                   gcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCATGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTg                                                                   >  1:523917/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                                                      <  2:795626/150‑1 (MQ=255)
                                                                                                gATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCg                                                      <  1:1595984/150‑1 (MQ=255)
                                                                                                  ttACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGGCAGCGGc                                                    <  1:641551/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:281964/150‑1 (MQ=255)
                                                                                                        aCGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGt                                              <  1:597295/150‑1 (MQ=255)
                                                                                                           ttGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCATGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGat                                           <  2:523917/150‑1 (MQ=255)
                                                                                                            tGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGctt                                          >  2:533018/1‑150 (MQ=255)
                                                                                                             gAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattc                                         >  1:1586345/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)
                                                                                                                    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)
                                                                                                                              cGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACcg                        >  2:1204096/1‑150 (MQ=255)
                                                                                                                                 gCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCgcgc                     <  2:86631/150‑1 (MQ=39)
                                                                                                                                 gCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGGGGCAGCGGCTACCGGGATTCCGTCACGCCAACAACCgcgc                     >  1:1206045/1‑150 (MQ=14)
                                                                                                                                       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:395943/1‑150 (MQ=14)
                                                                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc          >  2:402092/1‑150 (MQ=14)
                                                                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc          >  1:967617/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)
                                                                                                                                                    cggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCAtgtg  >  1:52622/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGGTTCCGGTACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTG  >  NC_000913/357344‑357641

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