Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 10,973 T→A K128N (AAA→AAT
I48I (ATT→ATA
yaaW ←
mbiA →
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091310,9730TA97.1% 89.0 / NA 34K128N (AAA→AAT
I48I (ATT→ATA
yaaW
mbiA
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (14/19);  minor base C (1/0);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTC  >  NC_000913/10726‑11218
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
tGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATAttc                                                                                                                                                                                                                                                     >  1:271220/1‑251 (MQ=255)
       tttAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTAAGTTCCTCGTCCATATTCTTCCAgg                                                                                                                                                                                                                                              >  2:337061/1‑251 (MQ=255)
                            ggccaaaggggcaccctggcccGCGGCGCGTAGTAAATCATGCCCAAGCACGCTCATGGCTGAGGGGGTGCGTAAAATGCTGGTAAATTGGCGTGAAAGCAAATGCGACACGTCTTTTGCCAATAATGTGTCTTTCATCAGCAGCGGCATCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACGCCTGCTGATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTcc                                                                                                                                                                                                                          <  2:360044/231‑1 (MQ=255)
                             cccaCAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAg                                                                                                                                                                                                                        <  1:37449/248‑1 (MQ=255)
                             ccaaCAGGGCCTCCCAGACCCGCGCCGCACAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTACAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTATTACAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCATCGTCCATATTCTTCAAGGTGTTACGCAGAAATTGTTcc                                                                                                                                                                                                                          >  2:316429/4‑249 (MQ=255)
                                    ggCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGt                                                                                                                                                                                                                 >  1:73749/1‑251 (MQ=255)
                                        ccccAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCt                                                                                                                                                                                                             <  1:63033/250‑1 (MQ=255)
                                                                       ccAAGCACGCTCATTGATGCGAGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCAAATAATTTGTCTTTCATCCGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGAACGCTGGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTAAAACGTAGACATCTCTTTGTCggt                                                                                                                                                                              >  2:213282/1‑250 (MQ=255)
                                                                        cAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCt                                                                                                                                                                             <  1:38405/251‑1 (MQ=255)
                                                                                 ccatggctgcGTGGGTGCGTAAAATGCGGGTCAATTTGCGGGAAAGCAAATGCGACACGACTTTTGCCAATAATTTGTCTTTGAGCAGCAGCGGCAGCAGGTCTGCCAGCTCAGGCGCCCGGGCATGGACCGCGTGCAGAAACGCCTGCTGATGTTCCGCGATCATATGCGTCCAGGTGTTACGCAGAAATTGTTCGAGCAACTGTTGCACAAGCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCtgc                                                                                                                                                                    <  1:358813/246‑3 (MQ=255)
                                                                                  cATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATGTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCtt                                                                                                                                                                     <  2:195030/249‑1 (MQ=255)
                                                                                                                    tGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATCATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGCCATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACCTCACACAAAATGACCCGAt                                                                                                                                 >  2:64789/1‑250 (MQ=255)
                                                                                                                                         acgccatttGCCAATAATTTGTCTTTCATCAGCAGCGGAAGCAGCTCTTCCAGCTCAATCACACTGGCATAGACCGCGTGCAGAAACTCCTGCTTAGGTTCCTCCTCCATATTCTTCCAGGGGTTACGCAGAAATTGTTCAAGAAACTGTTGCTCAATCTCAAACGTATACATCTCTTTGTCGGCTTTCAGCTTAAATCGGTTTACACATTCGAGCAAATTGGCCAGATAAAATTTAcaggatcatcgc                                                                                                             >  2:100202/7‑238 (MQ=255)
                                                                                                                                                 tGCCAATAATTTGTCTTTCATCAGAAGCGGCAGCAGCTCTTCCAGCCAATTCACCCTGGCATCGACCGCTGGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGAGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTCGACATCTCTTTGTAGGCTTTCAGCTTCAATCGTTTTGAAACATAGAGCAACATAGCCCGCTACAATTTTCCGTGTCCGCGCAGATTGTTgg                                                                                                   >  2:282430/1‑251 (MQ=255)
                                                                                                                                                  gCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTgg                                                                                                   <  1:169739/250‑1 (MQ=255)
                                                                                                                                                         atTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACt                                                                                            >  1:240287/1‑250 (MQ=255)
                                                                                                                                                                    atcaGCAGCGGCAACAGGTCTTCCAGCTCATTCACCCTGGCCTCGACCGCGTGCAGAAACTCATGCTTATGTTCCTCGTCCCTCTTCTTCCAGGTGTTACGCAGACATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACAGTGTCAGCGCAGTTTGTTGGCGATACTATCGCCACCaa                                                                                 >  2:50711/1‑250 (MQ=255)
                                                                                                                                                                      cagcagCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAAt                                                                              >  1:75118/1‑251 (MQ=255)
                                                                                                                                                                                       ctTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTccagc                                                              >  2:275591/1‑247 (MQ=255)
                                                                                                                                                                                               tCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTg                                                     >  1:163471/1‑251 (MQ=255)
                                                                                                                                                                                                cATTCACCCTGGCATCGACCGCGGGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTg                                                     <  1:286077/250‑1 (MQ=255)
                                                                                                                                                                                                     aCCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGtt                                               >  1:256519/1‑251 (MQ=255)
                                                                                                                                                                                                       ccTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTgc                                             <  1:265133/251‑1 (MQ=255)
                                                                                                                                                                                                         tGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTgcgg                                           >  1:291824/1‑251 (MQ=255)
                                                                                                                                                                                                          ggCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTgcggc                                          <  1:323050/251‑1 (MQ=255)
                                                                                                                                                                                                           gCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTgcggcg                                         <  1:321619/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cgtcttctgattcTCGTCCATACTCTTCCATGTGTTTAGCATAAATGGTTCCAGTACCTGTTTATTAATCTCAAACTTAGACAGCTCTTTGTCGTCGTTTAGCGTCAATCGTTTAGAAACATCGAGCAAAATTGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGTCACCAAAATGCTGTAATGCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTGCCa                  <  2:202864/239‑1 (MQ=255)
                                                                                                                                                                                                                                    gCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTTCGGCGATGTTGCTCGGGATGCCCTTCCAt                 <  1:122317/250‑1 (MQ=255)
                                                                                                                                                                                                                                    gCTTATGTTCCTCGTCCATATTCTTCCAGGGGTTACTCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCTATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCAt                 <  1:213282/250‑1 (MQ=255)
                                                                                                                                                                                                                                       tttGTTCCTCGTCCCTATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTTAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCATCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATTCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGa              <  1:149992/248‑1 (MQ=255)
                                                                                                                                                                                                                                             ccTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAAc       <  1:378606/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACa      <  1:337061/249‑1 (MQ=255)
                                                                                                                                                                                                                                                cGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAAGTTACCGTGTTCGCGCAGTTTGTTGGCGATACTATCGCCTCCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAg     <  1:353339/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTc  <  1:275591/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
TGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTC  >  NC_000913/10726‑11218

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: