Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,226 G→A R188R (CGC→CGT ytjB ← putative inner membrane protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,2260GA93.9% 82.0 / ‑5.8 34R188R (CGC→CGTytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/17);  minor base T (1/1);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

CTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCG  >  NC_000913/4623985‑4624469
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
ctacCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTTACAGGGGGGGCGATTGCTGCCAGCGGGGCCGTTTAccc                                                                                                                                                                                                                                            >  2:219822/1‑251 (MQ=255)
  acCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAGAGGCGCTTGCTGCCAGCGGGGACGTTTATCCTg                                                                                                                                                                                                                                          >  2:219282/1‑251 (MQ=255)
         cGGCATTGCGACAGAGAAACAGAACGCGCTGCGGGGCGGGCATCTGGCGAAAAATACACTACCCCACCGCCAGGGTAAAACACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTGACCCTGTagcag                                                                                                                                                                                                                                    <  2:226647/250‑1 (MQ=255)
                gcgcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTTTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCgg                                                                                                                                                                                                                             <  1:120080/250‑1 (MQ=255)
                gcgcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCgg                                                                                                                                                                                                                             <  1:267869/250‑1 (MQ=255)
                  gcCAGCGAAACAGAACGGCCTGCCTGGCGGGCATCTGGCGTAAAATACACTCCTCCCCCGCCCGGTTAAACCACGGGTCGTGAGAGTCAGAGATGCGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGCAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACCGGGAAGGCGCTTACTGCCCGCGGGTAAGTTTAccccatcgccgcaggcgcgtt                                                                                                                                                                                                                          >  2:20263/1‑233 (MQ=255)
                       agaAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGGGGGCGGGTccgccc                                                                                                                                                                                                                     >  2:266490/1‑246 (MQ=255)
                                             cGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGGTTGcccgggaa                                                                                                                                                                                                  >  2:178166/1‑242 (MQ=255)
                                                  atctGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTTCCAGTGAAAgcagcac                                                                                                                                                                                          >  1:245923/4‑250 (MQ=255)
                                                      ggCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCGTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATc                                                                                                                                                                                      <  1:178166/251‑1 (MQ=255)
                                                         gAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATccgg                                                                                                                                                                                   >  1:38959/1‑248 (MQ=255)
                                                            aaagacacTCCTCCTCCGCCAGGTTAAACCACGGTTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTCTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCATCCGGCAGCGGTGTGCTGGCGGTTAAAAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTTCCCTGTAGCAGCGTGCGGGGCAGCACTAAGCCGATTGCCAGTGATAGCAGCACCATTAGGCGt                                                                                                                                                                                <  2:230059/247‑1 (MQ=255)
                                                                     ccccccccGCCAGGTTAAACCACCGGTCGTAAGAGTCAGAGACGAGCAGGCGTAATGTGGACCTAACGATGTCCTTTCTTGTAGTAATGGGTCACTCGTTTTTCTCGCTTTCTTCCGCTTCCGGCACCGGTTTTCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCGGCGTGCGGGTCAGTACAACTCCGATTGGCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTa                                                                                                                                                                       <  2:49306/245‑1 (MQ=255)
                                                                        ccACCGCCAGGTTAAACCACGGGTCGTAAGAGGCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAgtg                                                                                                                                                                    <  2:184823/251‑1 (MQ=255)
                                                                            cGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCg                                                                                                                                                                                                                              <  1:182770/189‑1 (MQ=255)
                                                                            cGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGTAGCAGCGGGCg                                                                                                                                                                                                                              >  2:182770/1‑189 (MQ=255)
                                                                                    tAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTAAGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTCAAATGTTAGTGGGGTTATCCAcc                                                                                                                                                        >  1:95098/1‑251 (MQ=255)
                                                                                          aCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGGTATCCACCTGTTgg                                                                                                                                                  <  1:219822/251‑1 (MQ=255)
                                                                                             ggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTTCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCt                                                                                                                                               <  2:247931/251‑1 (MQ=255)
                                                                                                                    aGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTc                                                                                                                        <  1:189357/251‑1 (MQ=255)
                                                                                                                               ggACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTGTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGGATGGGGGTCGAGCGTCAgg                                                                                                              <  2:38959/250‑1 (MQ=255)
                                                                                                                                       cGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTTACAGGAAGGGCGATTGCTGCCAGCGGGTAAGTTTTCCCTGTAGCAGCGGGCGGGGTAGCACTACGCCGATTTGCAGTGAAAACCGCAACATCAGGCGGAAAATGTTAGGGGGGTTATTCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGGCGAGCGTCAggggacgata                                                                                                      >  2:240094/1‑242 (MQ=255)
                                                                                                                                               tttCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTACCAGGAAGGGCGATTTCTGGCAGCGGGTTCGTTTACCCTGTCGCAGGGGGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCCGGCGGAAAATGTTTGTGGTGTTATCCACCTGTTGGGCTTTGGTGGCGAGCGTATGGGTGTCGAGCTTCAGGCGCAGATAGCCGAGCgg                                                                                             >  2:98004/1‑251 (MQ=255)
                                                                                                                                                     gtagtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGt                                                                                        <  1:266490/250‑1 (MQ=255)
                                                                                                                                                           aTGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttccc                                                                                 <  1:150506/251‑4 (MQ=255)
                                                                                                                                                            tGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGGTTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGGGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAACCGAGCGGTCCGtttttcccg                                                                                >  1:68577/1‑247 (MQ=255)
                                                                                                                                                                           tctcGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGTTCAGCACTACGCCGATTGCCATTGAAAGCAGCAACATCATGCGTAAAATGTTAGTGGTGTTATCCACTTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGTCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGAc                                                                  <  1:172193/250‑1 (MQ=255)
                                                                                                                                                                           tctcGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGAGGTTCTCAGGAAGTTCCTTTGCTTCCAGCGGGTACGTTTACCCTGTAGCCGCGTGCGGGTCCGCACTACCCCGCTTGCCAGTGAAAGCAGCCACATACGGGGAAAAATGGTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGGGCGTCCGGCGCAGATAACCGAGCGGTCCGtttttcacggctattggttcgtc                                                                  >  2:185861/1‑232 (MQ=255)
                                                                                                                                                                                             ccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATa                                                <  1:59681/250‑1 (MQ=255)
                                                                                                                                                                                                                       aGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTCAAATAGCCGCCCGCTTTTTTACCGTCGAgcgc                     >  1:248956/1‑251 (MQ=255)
                                                                                                                                                                                                                             ggCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCAGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCGGTGGGGCCACGGTGGCGAGCGTATGGGTGTCGAGCGGCAGGCGCAGATTGTCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGTCAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCg               <  1:20263/251‑1 (MQ=255)
                                                                                                                                                                                                                              gCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAACAGCAAAATCAGGCGTAAAATGTTAGTGGTGTTATCCAACTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCAGCCCGCTTTTTTACCGTCGAGCGACAGCCgg              >  1:16028/1‑251 (MQ=255)
                                                                                                                                                                                                                                gagtgttgcCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTTCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGAAATTGGTTCGACAATCTGCTTGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCGGCCGGTc            <  1:86111/245‑1 (MQ=255)
                                                                                                                                                                                                                                          aGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCg  >  1:209745/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
CTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCG  >  NC_000913/4623985‑4624469

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