Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,614,387 C→T intergenic (+127/+293) yjjV → / ← yjjW putative metal‑dependent hydrolase YjjV/putative glycyl‑radical enzyme activating enzyme YjjW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,614,3870CT91.5% 117.2 / 2.1 47intergenic (+127/+293)yjjV/yjjWputative metal‑dependent hydrolase YjjV/putative glycyl‑radical enzyme activating enzyme YjjW
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (19/24);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCA  >  NC_000913/4614095‑4614682
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                       
aaCCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCGGCCTAACCGCCCNGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCNTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGt                                                                                                                                                                                                                                                                                                  <  2:645122/300‑1 (MQ=255)
            ggATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATg                                                                                                                                                                                                                                                                                      >  2:201890/1‑300 (MQ=255)
             gATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGc                                                                                                                                                                                                                                                                                     <  1:108448/300‑1 (MQ=255)
                atGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTg                                                                                                                                                                                                                                                                                  <  2:418398/300‑1 (MQ=255)
                             gTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGCTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGt                                                                                                                                                                                                                                                                     >  1:586124/1‑300 (MQ=255)
                               tttCATGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACTGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCATTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTcg                                                                                                                                                                                                                                                                   <  2:265388/300‑1 (MQ=255)
                                ttCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTcgc                                                                                                                                                                                                                                                                  >  1:621861/1‑300 (MQ=255)
                                          gCCTAACCGCCCGGAGCAGGCTGCCCGTGTGNTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGNCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCAt                                                                                                                                                                                                                                                        >  2:495846/1‑300 (MQ=255)
                                          gCCTAACCGCCCGGAGCAGACTGCCCGTGTGTTCGCCGTGTTTTGCGCGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCAt                                                                                                                                                                                                                                                        >  1:668173/1‑300 (MQ=255)
                                            ccaaCCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGAGGAGATGGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGc                                                                                                                                                                                                                                                      <  2:50818/298‑1 (MQ=255)
                                                  gCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACaa                                                                                                                                                                                                                                                <  1:670806/300‑1 (MQ=255)
                                                     cGGAGCNGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGNGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTc                                                                                                                                                                                                                                             <  2:577716/300‑1 (MQ=255)
                                                     cGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGGTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTc                                                                                                                                                                                                                                             <  1:164288/300‑1 (MQ=255)
                                                      ggAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAGTCt                                                                                                                                                                                                                                            >  1:713318/1‑300 (MQ=255)
                                                       gAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACaggtttg                                                                                                                                                                                                                                           >  2:453448/1‑294 (MQ=255)
                                                          cAGGCTGCCCGTGGGTTCGCCATGCTTTGCGAGTGGCGCCGGGAACCGGCGGATGAGATGGAGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCTTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGc                                                                                                                                                                                                                                        <  1:405984/300‑1 (MQ=255)
                                                                       cgtTCTCCGTGCTTTGCTAGTTGTGCCGGGAACCGGCGGATGAGATTGCGCAATCGTTGCTTAATAACACGTATACTTTGTTTAACGTGCCGTAGGCCGGAGAAGGCGTTCACGCCGCATCCGGCAGTTGGCNCACAGTGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGcc                                                                                                                                                                                                                           <  2:169394/299‑1 (MQ=21)
                                                                                tGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGc                                                                                                                                                                                                                  <  1:701357/300‑1 (MQ=255)
                                                                                 gCTTTGCGAGTTGCGCCGGGAACCGGCGGATNAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCNCAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGGCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCgtatccggca                                                                                                                                                                                                                 >  2:502736/1‑291 (MQ=255)
                                                                                                  gggAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                                                                                <  1:205381/300‑1 (MQ=255)
                                                                                                  gggAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGGTGTTTAACGTGCCGTATGCCGGATAAGGCGTTCACTCCGCATTCTGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                                                                                <  2:675086/300‑1 (MQ=255)
                                                                                                                                        aacgcatataCGTTGTTTAATGTGCCGTAGGTCGGATAATGCGTTCACTCCGCATCCGGCAGTTGGCGCACAATGCCTGATTCAACGCTTAACGCGTCTTATCGTGCCTACAGTTTTGTGCCGAACCGTAGGCCGGATAATGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTc                                                                                                                                                          <  1:58772/294‑1 (MQ=17)
                                                                                                                                            cGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTNTTATCATGCCTAAAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCAGCATCCGGCAGTCGGCGCACAACGCCTGATGCGACGCTTGTCGCGTCTTAt                                                                                                                                                      >  2:618292/1‑300 (MQ=255)
                                                                                                                                              gataCGTTGTGTAACGTGCTGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTTGGCGCACAATGCCTGAGGCGACACGTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGTCGTTCACGCCGCACTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGGCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCa                                                                                                                                                    <  2:685843/299‑1 (MQ=21)
                                                                                                                                                 gcagtatttAATGTACCGTAGGCCGGATAAGGCGTTCACGTCGCATCCGGCAGTTGGCACACAATGCCTGGGGCGACGCTTAGCGCGTCTTATCATGTCTACAGGTTTGTGCCGAACCGTAGGCCGGATAATGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGc                                                                                                                                                 <  2:258089/294‑1 (MQ=14)
                                                                                                                                                                    aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTAGCGTCTTATCATGCCTACAGGTTTGTGCCGAAc                                                                                                                              >  1:304832/1‑300 (MQ=255)
                                                                                                                                                                         ggATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAAGCCACATCTGGCAGTCGGCGTCCAATGACTGATGCGACGTTCGTCTCGTCTTCTCGTGCATACAGGTTTGTGCCGAACCGTAg                                                                                                                         >  2:670148/1‑300 (MQ=17)
                                                                                                                                                                            tAAGGCGTTCACGCCGCATCCGGCAGTTTGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGcc                                                                                                                      <  1:574891/300‑1 (MQ=255)
                                                                                                                                                                            tAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGcc                                                                                                                      <  1:165755/300‑1 (MQ=255)
                                                                                                                                                                               ggCGTTCACACCGCATCCGGCAGTTGGCGCACAATGCATGATGCGACGCTTACCGCGTCTTATCACGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCAGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAATCGTAGGCCGGATAAGGCGTTCATGCCGCATCCGGTAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCGTGCCCAAAAGCTCGTGTCGAACCGTAGGCCGGa                                                                                                                   >  2:319929/1‑300 (MQ=17)
                                                                                                                                                                                gtgTTCACGCCGCCTCCGGCAGTGGGCACACAATGCCGGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTATGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCGGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGttctgtctcttatacacatctcgc                                                                                                                  <  1:247/298‑23 (MQ=17)
                                                                                                                                                                                      aCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTAGCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCg                                                                                                            <  1:106647/300‑1 (MQ=255)
                                                                                                                                                                                                        ggCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGNATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCa                                                                                          >  1:506134/1‑300 (MQ=255)
                                                                                                                                                                                                         gcgcACACTGCCTGATGCTACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTACCGAACCGTAGGCCGGATAAGGCGTTCACGTCGCATTCGGCAGTTGGCGCNCAATGCCTGGTGCGACGCTTGTCGCGTCGTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAg                                                                                         <  2:608330/300‑1 (MQ=255)
                                                                                                                                                                                                           gcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGtt                                                                                       >  1:390622/1‑300 (MQ=255)
                                                                                                                                                                                                            cccaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTATGCCGCACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCAGCGCATAATGTCTTATGCGAAGCTTGTTGCGTCTTATCATGCCTAAAGGTTTGTGCCGAACCGTAGGTAGGATAGGGCGGTCGAGCCGCATCCAGCAGtcg                                                                                      >  2:370615/3‑298 (MQ=17)
                                                                                                                                                                                                                      gATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGNGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATg                                                                            >  1:497334/1‑300 (MQ=255)
                                                                                                                                                                                                                                        tctccccattcctgcAGGTTTATGCCGAACCATATGCCGGATAAGTCGTTCACGCCGCATTCGGCAGTTGGAACATAATGCCCGGTGCGACGTTTGTCGCGTCTTATCATGCCGACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGAc                                                          <  2:63033/286‑1 (MQ=17)
                                                                                                                                                                                                                                                        ggTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATCAGGCGTTCGCACCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACgtctgatgcgtgttctcgggcat                                          >  1:455868/1‑278 (MQ=17)
                                                                                                                                                                                                                                                        ggTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCt                                          >  1:136041/1‑300 (MQ=255)
                                                                                                                                                                                                                                                             gtgCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAg                                     >  1:179270/1‑300 (MQ=255)
                                                                                                                                                                                                                                                               gCCGAACCGGAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTATGCCGGATAAGGCGTTCACGCCACATCCGGCAGTCGGCGCATAATGCCTGATTCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTCATCAGGCCTACAAGTc                                   <  2:361942/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccGTAGGCCGGATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGcc                             >  1:657275/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                              ggATAAGGCGTTCACGCCGCATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCCACGCTTGTCGCGTCTTATCATGCCTACCGGTTTGTGCCGAACNGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCCTATCAGGCCTACAAGTCTGTGTCGAACCGTAg                    >  2:616922/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                aTAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGc                  <  1:461409/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       gTTCGCGCCGCATTCGGCATTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATcc           <  2:272059/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                gcATTCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTca  >  1:74434/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                       
AACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCA  >  NC_000913/4614095‑4614682

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: