Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 11,814 G→A intergenic (‑28/‑349) yaaI ← / → dnaK DUF2541 domain‑containing protein YaaI/chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091311,8140GA94.4% 94.0 / ‑1.9 36intergenic (‑28/‑349)yaaI/dnaKDUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (17/17);  total (18/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 = 4.79e-01

AAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCC  >  NC_000913/11573‑12035
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                             
aaaCCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAAt                                                                                                                                                                                                                      >  1:481609/1‑251 (MQ=255)
 aaCCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATa                                                                                                                                                                                                                     <  1:406806/251‑1 (MQ=255)
     gATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATAccc                                                                                                                                                                                                                  >  1:349503/1‑250 (MQ=255)
          aCCGCTTAACTGTAAATCTCCATGGGCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGGCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGt                                                                                                                                                                                                            <  2:375659/251‑1 (MQ=255)
               ttAAGTGTAAATCGCCATGGTCGGCAGAGAGTTTTATGCGTTTCACAATGCGACAAACAGGAATTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTAc                                                                                                                                                                                                        <  2:231937/250‑1 (MQ=255)
                           cGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAag                                                                                                                                                                                             >  1:168212/1‑249 (MQ=255)
                               aTGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGTTAAAGagaagg                                                                                                                                                                                       >  2:293793/1‑250 (MQ=255)
                               aTGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGagaagg                                                                                                                                                                                       <  2:30140/251‑2 (MQ=255)
                                    cggtcgagagTTGGATGCGTGTCACAATGCGTCAAACAGGAAGTTCCAGCGTCAGTACGTTGGTTTCGTTACGCTGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGTCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTTCTGCTGATAAAGAGAAGGTTCa                                                                                                                                                                                   <  2:93889/245‑1 (MQ=255)
                                    cGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCTCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAg                                                                                                                                                                                  <  1:239905/251‑1 (MQ=255)
                                    cGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAg                                                                                                                                                                                  <  1:292295/251‑1 (MQ=255)
                                              ttGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAACGAGAAGGTTCAGGTACGTAAtc                                                                                                                                                                        >  1:237765/1‑251 (MQ=255)
                                                        ttCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATcttttctttt                                                                                                                                                               <  1:387390/250‑1 (MQ=255)
                                                              aTGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTttattaa                                                                                                                                                         >  2:261247/1‑249 (MQ=255)
                                                                    cAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAAttttt                                                                                                                                                  <  2:428583/251‑1 (MQ=255)
                                                                     aaaCAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGACATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACCGCATCAGGCTAATCGCCAGGCTGGCGGACATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGACGAGGTAGAGGGATGAATACCCACTGGTTTACTGCTGTTAAAGAGAAGGTTCAGGGACGTAATCTTGTCTTTTTATTACAAtttttt                                                                                                                                                 >  2:323199/1‑251 (MQ=255)
                                                                           gAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTtgatga                                                                                                                                            >  1:128602/1‑250 (MQ=255)
                                                                                     cgcCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCt                                                                                                                                 >  1:258097/1‑251 (MQ=255)
                                                                                                                   tGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCGGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCGCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGGTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGAGTCATTCTTTATCTGAATaaaactg                                                                                                    <  2:238904/250‑4 (MQ=255)
                                                                                                                     cAATGACGCCGAGGAGATTATGGACGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGGTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGCTGAATGCCTTGGCTGCGATTCATTCTTTTTCTGAATAAAATTGCTg                                                                                                 >  2:376013/1‑251 (MQ=255)
                                                                                                                                          ggTCGGTTGCCCGCGCCGTGCAGCCAAGCATCAGGGTAATCGCCCGGCTGGGGGAAATCGGAAAAACGGATTTCATAAGGATTTTCTTAGTGGGAAGAGGTAGGGGGAAGAATTCCCAATAGTTTTCTGCTGATAAAGAGAAGATTCAGGGACCTAATCTTTTTTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTg                                                                            <  2:290778/251‑1 (MQ=255)
                                                                                                                                                           gTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAAAGGATTTCATAAGGATTCTCTTAGTGGGCAGAGGTAGAGGGATGCATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGTATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAAtttt                                                           >  1:282808/1‑250 (MQ=255)
                                                                                                                                                              cagcacagcaTCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAggta                                                                                                                                                                               <  1:410620/132‑3 (MQ=255)
                                                                                                                                                              cagcacagcaTCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAggta                                                                                                                                                                               >  2:410620/1‑130 (MQ=255)
                                                                                                                                                                 cacagcaTCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACATTTTTTTGATGAATGCCTTGGCTGCGTTTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTAAGTCTTGTCCTGCCATATCGCGCAATTTCTgcgc                                                      >  2:318996/1‑250 (MQ=255)
                                                                                                                                                                           ggCTAATCGCCAGGCTGGCGTAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGTTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCCGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACaaa                                           <  1:174847/251‑1 (MQ=255)
                                                                                                                                                                            gCTAATAGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACaaaa                                          <  1:198659/251‑1 (MQ=255)
                                                                                                                                                                               aaTCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACaaaaaa                                        <  1:261247/250‑1 (MQ=255)
                                                                                                                                                                                  cGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGACTAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAAtttt                                    >  1:292306/1‑251 (MQ=255)
                                                                                                                                                                                       ggCTGGCGGAAATCGTCAAACCGGATTTAATACGGATTCTCTTAGTGGGAAGAGGTAGAGGATTGAATACCCACTAGTTTACTGCTGCTAAAGAGAAGGTTCCGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGTTTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTATGCGAAACAGCACAAAAAATTTTTGCAt                               >  2:313007/1‑251 (MQ=255)
                                                                                                                                                                                                    agTAAAAACGGTGTTCATAGGGATTATCTTAGTGGGAAGAGGGAGAGGGATGAATACCCACTATATTTCTACAGATAAAGAGAAGGTTCAGGTACGTAATCTTTTTTTTTTATTACAATGATTTAATGAATGCCTTGGCGGCGATCCATTCTTTATCTGAATAAAATTGCTGTAAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTtgat                   <  2:80548/249‑1 (MQ=255)
                                                                                                                                                                                                          accGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAAt                                                                                            >  2:351719/3‑171 (MQ=255)
                                                                                                                                                                                                          aaCGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAAt                                                                                            <  1:351719/171‑1 (MQ=255)
                                                                                                                                                                                                              gATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGACTGCGATTCATTCTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTAc       >  1:16881/1‑251 (MQ=255)
                                                                                                                                                                                                                tttCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACg      >  1:198251/1‑251 (MQ=255)
                                                                                                                                                                                                                    aTAAGGATTTTCTTAGTGGGAAGAGGTAGAGGGAAGAATACCCACTAGTTTACTGCGGATAAAGAGAATTGTCAGTTACGTAATCTTTTCTTTTTATTACAATTTTTTGAGGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGCCAATTTTACTTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCGCCCCCTTGATGACGTGGGTTACGAccc  <  2:301696/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                             
AAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCC  >  NC_000913/11573‑12035

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: