Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,138,533 C→T G52G (GGC→GGT)  flgK → flagellar hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,138,5330CT97.0% 91.9 / ‑4.3 34G52G (GGC→GGT) flgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/13);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACG‑CAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCA‑CCAGCCTGC‑AAACGCTGGTGAGTAACGCGGAAG‑ACCCGGCAGCGCGC  >  NC_000913/1138290‑1138767
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
tGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCggttg‑ggttg                                                                                                                                                                                                                                            >  1:303568/1‑251 (MQ=255)
                     gAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTTTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGTATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCt                                                                                                                                                                                                                       <  1:176691/250‑1 (MQ=255)
                        tAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGGTGGCAATGGTGTC‑TACGTTTTTGgtg                                                                                                                                                                                                                   >  1:326269/1‑251 (MQ=255)
                                               cccGATAATACTTTGTAATTGAAGGCTTATAAGGAACCTCCATGGCCATCTTTATTAATAACGACATGAGCGGACTGAGCGCGGCCCAGGCGGCGTAAAATACGGAAAGAAATATTATCTTCAGCTAAAACTTTGCCGGAGATCCCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGttt                                                                                                                                                                                          <  2:271803/250‑1 (MQ=255)
                                                   ataataCTCTGTAATTGAAGTTTTAAGAGGACCCTCCATGTCAAGCTTGATTAATAACGCCATGAGCGGACTGACCGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATGATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTttatta                                                                                                                                                                                      <  2:161799/251‑1 (MQ=255)
                                                    taataCTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GCTTGGCAATGGGGTC‑TACTTTTCTGGTGGGCAGCGGGA‑GTATGAGTG‑GTttatta                                                                                                                                                                                      >  2:111277/1‑250 (MQ=255)
                                                       taCTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGATTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGTGGTGGGCA‑TGGTGTA‑TACGTTTCTGGTGTGCGGCGAGG‑GGATTA‑TGCGTttattaaca                                                                                                                                                                                   >  1:199378/1‑247 (MQ=255)
                                                             atAATGGAAGGATTATAAGGACCCCCCATGTCCAGCGTGATTAATAACGCCATGAGCGGAATGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTGGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGGGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCAGtt                                                                                                                                                                            <  2:326269/250‑1 (MQ=255)
                                                                              aGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGGTGGCAATTGGGTC‑CACGTTTCTGTTGTGTACAGTGG‑GTGTTA‑TGCGTTTATTACCAACCCGTCTCTTGCGGCGCAGACG‑c                                                                                                                                                          >  2:215545/1‑251 (MQ=255)
                                                                                   ccTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGGGTC‑CACGTTTCTGGTGTGCCGCGGGA‑GTATGG‑TGGGTTTATTACCAACCAGTTACGTgcggcggcggcgcca                                                                                                                                                         >  2:499534/1‑239 (MQ=255)
                                                                                        aggTCCAGCGTGATTAATAACTCCTTGATCGGATTGAACGCGCCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCAATAACGTTGCAGTATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACG‑CAAAGTAGCgg                                                                                                                                                <  1:272966/249‑1 (MQ=255)
                                                                                                       aataaCGCCATGAGCGGACTGAACGCGGCCCAGCCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTGGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACG‑CAAAGTAGCGGTCTGACTGCCCGCt                                                                                                                                  <  1:355095/250‑1 (MQ=255)
                                                                                                              ccATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATACTATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGATGGCAATTGGGTC‑TACGGTTCTGGTGTGCATCG‑GACGTATGA‑TGAGTTTTTTACCTACCAGTTACCTGCGGCGCCAACA‑CAACGTAGCGGTATGACCGCCCGCttagcgc                                                                                                                            >  1:61415/1‑243 (MQ=255)
                                                                                                               cATGAGAGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTG‑GGTTGGCAATGGAGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCAGTTACGTGCGGCGCCGACG‑CAAAGTTGCGGTCTTACTGCCCGcatagcgcca                                                                                                                          >  2:333849/1‑241 (MQ=255)
                                                                                                                                    ccAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAAACACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGCGTCATA‑GTTTCTGGTAGGCAGAGTGA‑GTATGG‑TGAGTTTATTACCAACCAGTTACGTGCGGCGCAGACGACAAA‑TAACGGTCTGACTGCACGCGCTGAGAAAGTGTCGAAAATCGACAAAATGc                                                                                                    >  1:334221/1‑251 (MQ=255)
                                                                                                                                      aGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCa                                                                                                                                                                                                               <  1:437286/145‑1 (MQ=255)
                                                                                                                                            cgTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGGTATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGACGCTGACGGGTG‑GGGTGGCTAGTGTGTC‑TCCGTTTGTGGGGTGCGGCGTGG‑GTAAGA‑TATATTTATTGCAAAACCGTTAAGTGCGGCGCCGGCG‑CCAAGTAGCGGGCTGTCTGCCCGCTTAGAGAAGATGGCGAACATCgtagaaaaggcagtcacc                                                                                            >  2:104467/1‑234 (MQ=255)
                                                                                                                                            cgTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCAGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGGGTC‑TACGGTTCTGGTGTGCAGCGTGA‑GTATGG‑TGCGTTTATTACCAACCCGTTACGTGGGGGGCACACG‑CAAAATAGCGATCTGGCTGCCCGCCATAAGCAGATGTCGAAAAtaggaaataagctccccccc                                                                                            >  2:305259/1‑232 (MQ=255)
                                                                                                                                                     cGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGGTGGCAAAGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCCGTTACGTGCGGCGCCGACG‑CAAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACctct                                                                                   >  1:486025/1‑251 (MQ=255)
                                                                                                                                                             aataatATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGG‑GTATGA‑TGCGTTTATTACCAACCCGTTACGTGCGGCGCAGACG‑CAAAGTAGCGGTCTTACTGCCCGCTCTGAGCAGATGTCGAAAATTGAACATATGCTCTCCACCCATACCTCTTCGCTGGc                                                                           >  1:395120/1‑251 (MQ=255)
                                                                                                                                                                       ccAGCTATAACGTTGCCGGATATCCCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACG‑CAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATg                                                                 <  1:147666/251‑1 (MQ=255)
                                                                                                                                                                         aGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGGGA‑Ggatga‑t                                                                                                                                                                                                >  1:457886/1‑123 (MQ=255)
                                                                                                                                                                                aCGTTGCCGAATATACCCGCCAAACCACTATGATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACG‑CAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATttc                                                        <  1:104547/251‑1 (MQ=255)
                                                                                                                                                                                      ccGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGAAAAAGTGTC‑TGCGTTTCTGGGGTGCAGCGGGA‑GTATGA‑TGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACG‑CAAAGTAGCGGTCTGAATGCCCGCTATGAGCAGATGTCGTAAATCGACAATATGCTCTCCACCAGAACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTTA‑cc                                                 >  1:325665/1‑251 (MQ=255)
                                                                                                                                                                                          atatACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGGGCAGCGGGA‑GGATGA‑TGCGTTTATTACCACACAGTTACGTGCGGCGCAGACG‑CACAGTAGAGGTCTGACTGCCCGCCATGAGCAGATGTCGAAAATCAACACTACTATCTCCCCCAGTAACTCTTCCCTGGCAACACAGATGCAGGCTTTCTTTA‑CCaccc                                             >  2:109591/1‑248 (MQ=255)
                                                                                                                                                                                               cccGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCAGt                                                                                                                                                                             <  1:158226/120‑1 (MQ=255)
                                                                                                                                                                                               cccGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGGTGGCAATGGTGTC‑TACCGTTCTGGGGTGCCGCGTGA‑GTATGG‑TGGGGTTATTACCAACCAGt                                                                                                                                                                             >  2:158226/1‑120 (MQ=255)
                                                                                                                                                                                                 cGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACG‑CAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCA‑CCAGCCTGC‑AAAc                                     <  1:408216/251‑1 (MQ=255)
                                                                                                                                                                                                  gtcAAACCTCTTTTATGGCGCAGGCCAATAGCACGTGGGGCGCTGGCGGTTG‑GGGTGGCAATGGTGGC‑GACGGTTCTGGTGGGCAGCGGGA‑GTAGGA‑TGCGTTTATTAACATCCAGTTACGTGCGGCGCAGACG‑CATAGTAGCGGTCTGACTGCACGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTACCCCATTTCCTCTTCGCTGGCAACACGGATGCAGGATTTCTCCA‑CCAGCCTGC‑aaccg                                    <  2:489243/249‑4 (MQ=255)
                                                                                                                                                                                                       aCCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGGTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACG‑CAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCA‑CCAGCCTGC‑AAACGCTggg                               >  1:361648/1‑250 (MQ=255)
                                                                                                                                                                                                                atGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTATCAACCAGTTACGTGCGGCGCAGACG‑CAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCA‑CCAGCCTGC‑AAACGCTGGTGAGTAAcg                       <  1:7982/250‑1 (MQ=255)
                                                                                                                                                                                                                          ccAATAGCACGTTGGGCGCTGGCGGTTG‑GGGTGGCACTGGCTTC‑TACGCTTATGGTGTGAAGCGGGA‑GTATAA‑TGCGTTTATTACCACCCAGGTACGTTCGGCGCAGGCG‑CAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAAATGACAATATGCTCTCCACCAGTCACTCTTCGCTGGCACCACAGATTCAGGATTTTTTCA‑CCAGCCTGC‑AAAAGATGGAGCGTAACGCGGGAA‑ACCCg           >  1:417519/1‑251 (MQ=255)
                                                                                                                                                                                                                           cAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTC‑TACCTTTCTGGTGGGAAGCGGGG‑GTATGA‑TGGGGTTATTAACAACCAGTTACGTGCGGCGCAGAAG‑CAAAGGAGCGGTCTGACTGCCCCCTACGAGCAGGTGTCGGAAAACGACAATATGCTCTACACCCGTACCTCTTCGCTGGCAACACAAATGCAGGATTTCTTCA‑CCAGCCTGC‑AAACGCTGATGAGATACGCGGTAG‑ACTCgg          >  1:266904/1‑251 (MQ=255)
                                                                                                                                                                                                                            aaTAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAAAGGTGTC‑TAAGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTTTTACCAACCAGTTACGTGCGGCGCAGACG‑CAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGGCAATATGCTCTCCACCAGTACCTCTTCGCTGGCACCACAGATGCAGGATTTCTTCA‑CCCGCCTGC‑AAACGCTGGGGGGGAACACGGAAG‑ACCCGGc         >  2:270109/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cGTTGGGCGCTGGCGGTTG‑GGCTGGCAATGGTGTC‑TTCGTTTCTGGTGTGCCGCGTGA‑GGATGA‑TGCGTTTATTACCAACCAGTTACGCGCGGCGCAGACG‑CAAAGTAGCGGTCTGACTGCCCGCTATGAGCCGATGTCGAAAATCGACAATATGACCTCCACCAGTACCTCTTCGCTGGCCACACAGTTGCAGGTTTTCTTCAACCA‑CCTGCCAAAC‑CTGGTGGGTAACCAGGAAGCACCCG‑CAgcgccc  >  2:167382/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
TGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTG‑GGTTGGCAATGGTGTC‑TACGTTTCTGGTGTGCAGCGTGA‑GTATGA‑TGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACG‑CAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCA‑CCAGCCTGC‑AAACGCTGGTGAGTAACGCGGAAG‑ACCCGGCAGCGCGC  >  NC_000913/1138290‑1138767

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