Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,169,836 A→G L180P (CTG→CCG)  ldtC ← L,D‑transpeptidase LdtC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,169,8360AG88.2% 72.0 / 1.0 34L180P (CTG→CCG) ldtCL,D‑transpeptidase LdtC
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (18/12);  total (20/14)
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.05e-01

GGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCG  >  NC_000913/1169590‑1170076
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                
ggCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTAtccggtcc                                                                                                                                                                                                                                              >  1:219679/1‑251 (MQ=255)
   ggattaagttCAACCAGCCGCGCACCGTTTGGTGCGGTGGAGTCTTTTATCGGATTGTTGATGATATTCTCTTTGTTGCCTGGGGTGACCTGGCTAACGAGTGTTTTGATATCGTCATCCCGTAGACGAATACAGCCAGAACTTACCCTCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTAtccggtccag                                                                                                                                                                                                                                            <  2:252975/242‑1 (MQ=255)
      tGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCgg                                                                                                                                                                                                                                        >  1:70285/1‑251 (MQ=255)
      tGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCgg                                                                                                                                                                                                                                        >  1:284165/1‑251 (MQ=255)
       gTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGc                                                                                                                                                                                                                                       <  1:140888/251‑1 (MQ=255)
       gTACTTCAACCAGACGCACACAGGTTGTTTCGGCTGAGACTTTTACCGGAGGGTTGATGATATTCACTTTGGTGCTTGGGGTGACCTGGCTAAATAGTGTTTGGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCATTGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGTCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGc                                                                                                                                                                                                                                       <  2:339152/251‑1 (MQ=255)
           ttcaacctatcgcgCCCCGTTTGGTCTGGCGTAGACTTTTATCGGAGTGTTGATGATATTCACTTGGGTGCCTGGGGTGCCCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACt                                                                                                                                                                                                                                    <  2:59122/240‑1 (MQ=255)
                   gACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGc                                                                                                                                                                                                                           >  1:157977/1‑251 (MQ=255)
                   gACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTAACCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCCGACGAATCGCATGATGGCCCATTGGGATATCCAGTCCAGCCGGCACTACCGCAGGc                                                                                                                                                                                                                           >  2:66286/1‑251 (MQ=255)
                      gcgcgctgttgggtccgtcAGATACTTGTATCGGAGTGTTTATGATATTCTCTTTGGTGCCTGGGGTGCCCTGGCTAAAGAGGTTATTGATTCCGTCATCCCGCAGCCGAAGACAGCCAGAACTTACCCGCATGCCAATGCCGAAACCGGCGTTCGTACCATGAAGCAAATAACCGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAAc                                                                                                                                                                                                                        <  2:147804/233‑1 (MQ=255)
                             gTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCaa                                                                                                                                                                                                                         >  1:246499/1‑243 (MQ=255)
                             gTTTGGTCGGGCAGAGACTTTTATCGGAGTGTTGATGATATGCACTTTGGTGCCTGGGTTGACCTGGCTAAAGAGTTTTTTGATAGCGTCAGACCGCAGACGAATACAGCCAGAATGTACCCGCATGCCAATGCCAAAATCGGCTTTCGTACCATGAAGCAAATAAACTCCGTCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTTTCGCAGGCaa                                                                                                                                                                                                                         <  2:246499/243‑1 (MQ=255)
                                     gggCAGAGACTTTTATCTGAGTGTTTATGATATTCCTTGTGGTGCCGGGGGTGACCGGGCTAAAGAGTGTTTTGATACCGTCAGCCCGTAGACAAATCCAGCCTGAACTTACCCGCATGCCAATGCCGAACTCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCt                                                                                                                                                                                                         <  2:211169/250‑1 (MQ=255)
                                       gCAGAGGCTTTTATCGGAGTGTTGTAGATATTAACGTTGGTGCCAGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCAGCCCGCATACGAATACAGCCAGAACTAACCCGTATGCCAATGCCGAAACCGGCGTTCGTACCATGAAGCAAAGAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCTAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCttt                                                                                                                                                                                                       <  2:120286/251‑1 (MQ=255)
                                           agaCTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGTTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCAAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATaa                                                                                                                                                                                                   <  1:96476/251‑1 (MQ=255)
                                                                     ttCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCAGCAGACGAATACAGCCAGACCTTACCCGCATGCCCATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGACATCGGCCGCCAGACGAATAGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCCATACCGCAGGCAACTCAATTCCCTGTGATTTATAACGGGCGCGGTTGGTTGACGTTAGCGt                                                                                                                                                                         >  1:360536/1‑251 (MQ=255)
                                                                     ttCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATAGTCATCCCGCAGACGAATACAGCCAGCACTTAGCCGCATGCCAATGCCGAAATCGGCGTTAGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGAATGATGGCCCATTGTGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTAAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGATGGCgg                                                                                                                                                                         >  2:251366/1‑250 (MQ=255)
                                                                                    tGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTCTAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGTTTTg                                                                                                                                                           >  2:176954/1‑251 (MQ=255)
                                                                                    tGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTg                                                                                                                                                           >  1:108033/1‑251 (MQ=255)
                                                                                    tGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCAGCCATAGGCCACCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCAAGGTTGGGTTTg                                                                                                                                                           >  1:95965/1‑251 (MQ=255)
                                                                                                               tGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGATGGGATTGCACGTTTGTCTGAAACGGTGGTCAcc                                                                                                                                 >  1:271585/1‑251 (MQ=255)
                                                                                                               tGATAGCGTCATCCCGCCTACGCATACAGCCCAAACTTACCCCCAGGCCTATCCCAAAATCGGCGTTCGCACCATGAAGCACATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCAc                                                                                                                                  <  2:108033/250‑1 (MQ=255)
                                                                                                                          tCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATACCGAAATCGGCGTTCGTACCATGAAGAAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGCCCCATTGGGTTATCCGGTCATGCCGACACTACCGCAGGCAACTCAATTCCACGTGCTTTATAACGGGCGCGGAAGTTGGCCGTTGGCGTCCAGGTTGGATTTGCAAGATTGTCGGAAACgaggggacacattttcgttgt                                                                                                                      >  2:317968/1‑231 (MQ=255)
                                                                                                                                    gAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTcc                                                                                                                                                                       >  1:251418/1‑192 (MQ=255)
                                                                                                                                       tACAGCCAGAACTTACCCGCATGCCAATGCCGAAATTGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTcacc                                                                                                          >  1:25966/1‑250 (MQ=255)
                                                                                                                                                  cTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGCTATCCGGTCCAGCCGGCCATACCGCAAGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTAGCGTCCAGGTTGGGATTGTACGTTTGTCTGAAACGGTGGCCACCATTGTCGGTGGCAGCGAGTCCACACATAACTgccc                                                                                              >  2:301516/1‑248 (MQ=255)
                                                                                                                                                      cccGCATGCCAAAGCCGAAATCGGCGTTCATACCCTGAAGCAAATACACGCCGACATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGACCTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGACCGCGGATGTTTGCCGTTGGCGTACAGGTTGGGCTTGCCAGTTTGTCTGAAACGGTGGTCCACATTGACGGTGTCAGCGTGTCACCACCTACCTgcacaatc                                                                                          >  2:359512/1‑244 (MQ=255)
                                                                                                                                                                 aTGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATa                                                                               <  1:66286/251‑1 (MQ=255)
                                                                                                                                                                                      ccATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTATCTGACCAATACCTATTGGATACACGGTTACCGAATTTTTAcc                                                          <  1:331190/251‑1 (MQ=255)
                                                                                                                                                                                                     aCGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAAt                                            >  1:198045/1‑250 (MQ=255)
                                                                                                                                                                                                      cgccgcCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCAcc                                                                                                                                 <  1:139384/164‑1 (MQ=255)
                                                                                                                                                                                                      cgcagccATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGGTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACggtggtaccc                                                                                                                                 >  2:139384/5‑160 (MQ=255)
                                                                                                                                                                                                        ccgccATAGGCCGCCAGACGAATGGAACGATGGCCTATTGGGTTATCCGGTCCAGCCGGCACTACCTCAGGCAACTCAATTTCCTGGGCTTTATAACGGGCGCGGATGTTTGCAGTTGGCGTCCAGGTTGGGTTTGCACGTTTTGCGGAAACGGTGGTCACCATTGGCGGTGTCAGCGTGTAAACACCTAACTGACCAATACCTATGAGATACACGGTTACCGAATTTTTTCCCGGCGGTTAGTAATagag                                        <  2:234284/251‑1 (MQ=255)
                                                                                                                                                                                                        ccgccATAGGCAGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGGGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATagcg                                        >  1:256078/1‑249 (MQ=255)
                                                                                                                                                                                                                                           cATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTcgc     <  1:176954/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCCCCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCCGCTCCGCAATGTTGCTCACAATGCCTTcgctcg  >  1:154574/1‑248 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                
GGCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCG  >  NC_000913/1169590‑1170076

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