Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,614,346 T→C intergenic (+86/+334) 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,3460TC88.9% 117.8 / 7.0 63intergenic (+86/+334)yjjV/yjjWputative metal‑dependent hydrolase YjjV/putative glycyl‑radical enzyme activating enzyme YjjW
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base C (35/21);  total (39/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

GCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAG  >  NC_000913/4614102‑4614573
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
gcgcCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGGGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGcc                                                                                                                                                                                                                               <  2:385827/251‑1 (MQ=255)
 cgcCGGATATGCCTCTCAACGGTTTACAGGGGCAGCCTACCCGCCCGGAGCAGGCGGCCTGTGTGTTTGCCGTGCTTTGCGATTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACTTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCg                                                                                                                                                                                                                              <  1:420888/251‑1 (MQ=25)
    cGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGGGGGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAAc                                                                                                                                                                                                                           <  2:580712/251‑1 (MQ=37)
      gatagcccGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCg                                                                                                                                                                                                                         <  1:66773/245‑1 (MQ=37)
      gATATGCCGCTAAACGGTTGTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGGGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCg                                                                                                                                                                                                                         <  1:6496/251‑1 (MQ=25)
                 caccGGTTTTCAGGGGCATCCTAACCGCCCGTGGCAGGCTGCTCTTGTGTTTGCCGTGCTTTGCGTGTTGCGCGGGGAACCGGCGGATGAGATTGTGTAAGTGTTGCTTAAGATCACGTATACGTTGTTTACCGTGCCGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGACCCGTAGGCCGGATa                                                                                                                                                                                                              <  2:630292/248‑1 (MQ=18)
                 cAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGa                                                                                                                                                                                                                <  1:222535/249‑1 (MQ=255)
                  aaCGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTGGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATaa                                                                                                                                                                                                             <  1:341533/251‑1 (MQ=255)
                   aCGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATa                                                                                                                                                                                                              >  1:586917/1‑249 (MQ=255)
                      gTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATAAGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCg                                                                                                                                                                                                         >  2:440193/1‑251 (MQ=255)
                       ttttCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTTTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGt                                                                                                                                                                                                        >  1:74620/1‑251 (MQ=255)
                         ttCAGGGGCCGCCTACCCGCCCGGACCAGGCTGCACGTGCGTTTGCCGTGCTTTGCGAGCTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATACCACGTATACGTTGTTTAACGTGCCGTAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGCTGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTc                                                                                                                                                                                                      >  2:95600/1‑251 (MQ=21)
                           cAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAc                                                                                                                                                                                                    >  1:397691/1‑251 (MQ=255)
                            aGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACCCGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                                                                                   >  2:523080/1‑251 (MQ=255)
                             ggggGAGCCTAACCGCCCGGATGAGGGCGCCCGGGTGTGTGCCGGGCGTTGGGAGGTGCGCCGGGAACCGGAGGATGAGAGTGCGCAAGCGTTGCTTAATAGCTCGGATACGTTTTTTAACGTGCCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTTTCGCGTCTTAGCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAc                                                                                                                                                                                                    <  2:381076/249‑1 (MQ=17)
                             ggggCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                                                                                  >  2:564291/1‑251 (MQ=255)
                                gCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGTGACGCTTGTCGCGTCTTATCATGCCTACAAATCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                                                                               >  2:335265/1‑251 (MQ=25)
                                 cAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGCTAAGGCGTTCAcgccgc                                                                                                                                                                                               >  2:136247/1‑250 (MQ=255)
                                    ccTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTACCGTGCCGTAGACCGGCTAAGGCGTTCACGCAGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCCTGCCTACAAGTCTGTGCCGAACCGTACGCCGGCTAAGGCGTTCACGCCGCCTcc                                                                                                                                                                                           >  1:623797/1‑251 (MQ=21)
                                     cTAACCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                                                                          >  1:239582/1‑251 (MQ=255)
                                       aaCCGCCCGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                                                                         <  1:170531/250‑1 (MQ=255)
                                         ccGCCCGGAGCAGGCTGCTCGTGTGTTGGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAg                                                                                                                                                                                      <  1:125193/251‑1 (MQ=21)
                                         ccGCCCGGAGCAGGCTGCCCGGGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGGAGGCCGGATCAGGCGTTCACGCCGCATCCGGccg                                                                                                                                                                                      >  2:126126/1‑249 (MQ=21)
                                             ccGGAGCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGtcgg                                                                                                                                                                                  >  1:15330/1‑248 (MQ=255)
                                                  gCAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATGGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGcat                                                                                                                                                                             <  1:335265/251‑2 (MQ=18)
                                                   cAGGCTGCCGGTGTGGTGGCCGGGCTTTGCGAGGTGCCCGGGGAACCGGCGGAGGGTATGGCGCACGCGTGGCTTAATAGCCCGTGTACGTTGTTGAACGGGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTTCCATCCCTACAATTCTGTGCCGAACCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAGtctgcgcata                                                                                                                                                                            <  2:33200/251‑10 (MQ=2)
                                                   cAGGCTGCCCGTGTGTTTGCCGTGTTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGcata                                                                                                                                                                            <  1:278243/251‑3 (MQ=18)
                                                   cAGGCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGCTGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGcat                                                                                                                                                                             >  1:216096/1‑249 (MQ=18)
                                                     ggCTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAAt                                                                                                                                                                          >  1:344176/1‑251 (MQ=255)
                                                       cTGCCCGTGTGTTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGc                                                                                                                                                                        >  1:604882/1‑251 (MQ=18)
                                                               gtgtTTGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCg                                                                                                                                                                >  1:427869/1‑251 (MQ=18)
                                                                     gCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACACTGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCtt                                                                                                                                                          >  1:164906/1‑251 (MQ=18)
                                                                        gTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGCTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGAAGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCGTCGGCGCATAATGCCTGATGCGACGCTTGTc                                                                                                                                                       >  1:294137/1‑251 (MQ=17)
                                                                         tGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCCTAATGCCTGATGCGACGCTTGTcg                                                                                                                                                      >  1:249541/1‑251 (MQ=17)
                                                                           cTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCTACGCTTGTcgcg                                                                                                                                                    >  2:231111/1‑251 (MQ=17)
                                                                           cTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTcgcg                                                                                                                                                    >  1:470984/1‑251 (MQ=255)
                                                                               gCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCtt                                                                                                                                                <  1:418525/251‑1 (MQ=255)
                                                                               gCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCtt                                                                                                                                                >  1:398735/1‑251 (MQ=255)
                                                                                 gAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACTCTTGTCGCGTCTTAt                                                                                                                                              <  1:176056/251‑1 (MQ=17)
                                                                                        gcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCt                                                                                                                                       >  1:119354/1‑251 (MQ=255)
                                                                                         ccGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTCGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCtc                                                                                                                                      >  2:428371/1‑250 (MQ=17)
                                                                                          cGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCCTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTAc                                                                                                                                     >  2:384056/1‑251 (MQ=14)
                                                                                           gggAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTGTCATGCCTACAAGTCTGTGCCGAACCGAAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATACCTGATGCGAAGCTTGTCGCGTCTTAACATGCCTAcc                                                                                                                                    >  2:604988/1‑250 (MQ=14)
                                                                                            gaaaCCGGCGGAGGGGATTGCGCAAGCGTTGTTTAATAACACGTACACGTTGTTTAAGGGGCCGTGGCCCGGATATGGGGTTCCCGCCGCATCCGGCGGTTGGCGCCCAATGCCGGAGGCGACGCTGGTCGGGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTGATCATGCCTACag                                                                                                                                   <  2:613116/249‑2 (MQ=2)
                                                                                             aaaCCGGCTGATGGGATTGCGCAAGCGTTGCTTATTAACACGTATACGTTGTTTAGCGGGCCGTAGCCCGGATAGGGTGTGCTCGCCGCATCCGGCAGTTGGCGCCCAATGCCGGGTTCGACGCTTGTCGCGTTTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACagg                                                                                                                                  <  2:227158/250‑3 (MQ=11)
                                                                                              aaCCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTATTATCATGCCTACaggt                                                                                                                                 >  1:325204/1‑248 (MQ=14)
                                                                                                   gcggATGAGATTGCGCAAGCGTTGCTTAATAACCCGTATACGTGGTTTAACGGGCCGTAGGCCGGATAAGGCGTTCACGCCGCACCCGGCAGTGGGCGCACAATGCCTGATGCGACGCTTGTCGCGTTTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTtgtg                                                                                                                            <  1:231111/251‑1 (MQ=11)
                                                                                                      gATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCg                                                                                                                         >  1:204733/1‑251 (MQ=11)
                                                                                                       aTGAGATAGCGCACGCGTTCCTTAATAACACGTCTACGTTCTTTAACGTGCCGTACGCCGCCTAAGGCGTTCACGCCGCATCCAACAATCGGCGCACAAAGCCTGATGCGACGATTGTCGCGGCTTATCCTGACTACAAGTCTATGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCACACAATGCCTGATGCGCCGCATATCACGTCTTaccacacctgcaggattgggccg                                                                                                                         >  2:359723/1‑228 (MQ=2)
                                                                                                         gagaTTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAAc                                                                                                                      >  1:473794/1‑251 (MQ=11)
                                                                                                          agaTTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTCACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAAcc                                                                                                                     >  1:506679/1‑251 (MQ=11)
                                                                                                              tGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGt                                                                                                                   >  1:171059/1‑249 (MQ=11)
                                                                                                                cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGc                                                                                                               <  1:40747/251‑1 (MQ=11)
                                                                                                                 gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGcc                                                                                                              <  1:334821/251‑1 (MQ=11)
                                                                                                                  cAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCg                                                                                                             >  1:276819/1‑251 (MQ=11)
                                                                                                                       gTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACACTGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAGGCGACGCATAATGCCTGATCCGACGCTTGTCGCGTCTTATCCTGCCTACAAGTTTGTCCCGaaacggcggccggttac                                                                                                        >  2:539640/1‑236 (MQ=11)
                                                                                                                        ttGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATa                                                                                                         >  1:257836/1‑249 (MQ=11)
                                                                                                                                        atccGTTGTTTACCGGGCCGTAGGCGGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                       <  1:169305/248‑1 (MQ=255)
                                                                                                                                        ataCGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                       >  1:500179/1‑251 (MQ=255)
                                                                                                                                          tcgtggtgcagcgggCCGGAGGCCGGGTAAGGCGTCCGCGCCGCATCCGGCAGTTGGCTCCAAAGGCCTGATGCGCCGCTTAACGCGTATTATAATGCCTACAGGTTGGGGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCTGCAGTTGGCGCACAATGCCTGATGCGACGCTTGGCGCGTCTTATCATGCCTACAGGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                      <  2:148976/236‑1 (MQ=14)
                                                                                                                                                                        ggCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTAGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGACTGTGCCGAACCGTAGGCCGGGTACGGAGTTCACGACGCATCCGGCCGTAGGCGCagaaggcccgagtcggccc                                                       >  2:426681/1‑233 (MQ=14)
                                                                                                                                                                                                       caATGCCTGATGCGACGCTTACCTCGTCTTATCATGCCTACAGGTTTGTGCCGAACCTAAGGCCGGATAAGGCTTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTtgt                        <  1:54304/251‑1 (MQ=255)
                                                                                                                                                                                                                             cgcgTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACCCCGCAGCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGCCTTATCATGCCTCCAGGTTTGTGCCGCACCGTCGGCCggctccc  >  2:457055/1‑246 (MQ=17)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
GCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAG  >  NC_000913/4614102‑4614573

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