Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 71,877 A→G 24.9% E176G (GAG→GGG)  yabI → DedA family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,8770AG24.9% 59.7 / 3.6 37E176G (GAG→GGG) yabIDedA family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (4/22);  new base G (7/3);  total (12/25)
Fisher's exact test for biased strand distribution p-value = 3.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAA  >  NC_000913/71660‑72099
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 
gCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGTATGCTGGATCTGCCGGTGGCCCAATTTATTCCGCCGACTATTTACGGCTGCCTGCTGTGGCCGCCGTTTTTCTTCTTGCCAGGGATTTTGGCGGGCGCGGCCGTCGGATTTTCTTCCCGGAATCAGAGCGGTGGGGTTTaaag                                                                                                                                                                                                                           >  1:209057/1‑226 (MQ=255)
 cGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGGTCGCTATTCCTGCCGGAATGCCGAGCGGGGGGTTTAAATggc                                                                                                                                                                                                                        >  2:155575/1‑227 (MQ=255)
     gCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAAGATTATCGGCTGCCTGCTGTTGCCGCCGTTTTACTTCCTGCCCGGGATTCTGGCGGGGGCGGCGGTCCATATTCCTGCCGGAATGCAGAGGGGTGGGGGTAAATggggggg                                                                                                                                                                                                                    >  1:293514/1‑223 (MQ=255)
                gCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCCATGGTGGCGGGAATGCTGGATCTGCCGGTGGCCAAATTTATTACCCCGAATATTATCCGCTGCCTGCTGTTGCCGCCGGTTTACTTCCCTTCCGGGGTTCTGCCGGGCGCGGGGATCCGACATCCGTCCGGAATGCCGAGAGGGGAGttttacagggtggcgctggcccacg                                                                                                                                                                                                         >  1:432211/1‑206 (MQ=255)
                gCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGGTTTTGGCGGGGGCGGCGGCCGATATTTCTTCCCGCAAGAGGAGCGGGGGGGTTAAAGGGTGGTGGGTGGCAAAAg                                                                                                                                                                                                         >  1:304889/1‑228 (MQ=255)
                  ttGTTCCCCGTTTTGGGGGGTTTTTTTGTTGGCCGGGGGGGTCGGCTGGTGCCATGGGTGGCGGGAATGCTGGACCTCCGGGGGGTTGAGTTTATTCCGCCGAATATTACCGGGTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGAGTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCg                                                                                                                                                                                                       <  1:325177/227‑1 (MQ=255)
                              cTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGCTTCTGGCGGGCGCGGGGATCGATATTCCTGCCGGAATGCAGGGCGGGTGGTTTGAATGGTTGCTGCTGGCAACAGCGGGGGTTTTgggg                                                                                                                                                                                           >  2:104262/1‑228 (MQ=255)
                               tGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTTTGGGCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGGGATCGCAATTCCTTCCGGCATGCAGAGCCGGGCGGTTACATGGTGGCTGCTGGTAACAGCGGGGtgttttgggg                                                                                                                                                                                          >  1:189722/1‑223 (MQ=255)
                                 gtcggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTTGCGGGCGCGGGGGTCCGTTATCCTCCCCGACAGCAGCGGGGTGGGGTTAAAGGGGTGTCGCGTGCAACACCGGGGttttttgggggg                                                                                                                                                                                        >  1:379966/1‑221 (MQ=255)
                                    ggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGGGAGTTTAAATGGTTGCTGCTGGCAACAGCGGGGTTTTTGGGGGTtggt                                                                                                                                                                                     >  2:448242/1‑228 (MQ=255)
                                      tcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTTGATCTGCCGGTGGCTAAATTTATTACGCCGGATATTATCGGCTGCCCGCTGTGGCCCCCGTTTTTATTCCTGCCCGGGATTCTGGCGGGCCCGGGGATTGCTATTTCTGCCCGAATCGCGAACGGGGGGTTTTAATGGTTGTTTCATGCACCCACGGGGTTTTTTTGGGTtggggg                                                                                                                                                                                   >  1:175273/1‑228 (MQ=255)
                                                  cccGACGGGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCt                                                                                                                                                                        <  2:432211/227‑1 (MQ=255)
                                                              gtgtgggccatgggtggGGGGAATGGTGGATCGGCCGGTGGCTAATTTTATTACGCCGAATATTGTCGGCTGCCGGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                           <  2:175273/216‑1 (MQ=255)
                                                              gtgggTGCAAATGGGGGCGGGAAGGCTGGATCGGCCGGTGGCTAAATTTATTACGCCGAATATTGTCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                           <  1:488268/225‑1 (MQ=255)
                                                              ggtGGTGCCATGGTGGGGGGTATGTCTGGTCCTGCGGGTGCGAAAATTTTTTAGCCCAATATTTATGGGGTGCCGTTGTGGCCCGCGGTTTTTCTTCCGGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                           <  1:345358/226‑1 (MQ=255)
                                                                                        tgtataagggccagggCTAAATTTATTCCGCCGAATATTATCGGGTGCCTGCGGGGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGGTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGttt                                                                                                                                  <  2:531130/213‑1 (MQ=255)
                                                                                        ggggtacTGCCGGGGGCTAAATTTTTTGCGCCGAATATTATCGGCTGCCTGCTGGGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGttt                                                                                                                                  <  2:209057/221‑1 (MQ=255)
                                                                                        cTGGATCTGCCGGTGGCTAAATTTATTACGGCGAATATTAGCGGCTGCCTGGTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGttt                                                                                                                                  <  2:452875/227‑1 (MQ=255)
                                                                                                         tatttttATTACGCAGAATTGTAGCGGCTGCCTGCTGGGGCCGCAGTTTTACTTTCTGCCAGTGGTTCTGGCGGGCGCGGCGATCGATACTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGATGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcg                                                                                                                <  1:330623/224‑1 (MQ=255)
                                                                                                           aatttcagccggcGAATTTTATCGGCTGACGGGTGGGGCCGCCGTTTTACTGCCTGCCCGGGATGCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGGTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcgcg                                                                                                             <  2:555668/216‑1 (MQ=255)
                                                                                                            aTTTATTACGCCGAATTTTACCGGCTGCCGGCGGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcgcg                                                                                                             <  2:341547/227‑1 (MQ=255)
                                                                                                                 ttaCGCCGAATATTAGCGGCTGCCTGCTGTGGCCGGCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGt                                                                                                        <  2:586567/227‑1 (MQ=255)
                                                                                                                                        cctgcTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTTTGGGGCAGCGGGTAAGCGGCGGACCGGTTTAGTCATTATTTTTGCCCCGGGGGGTTGTTGTGGGTGGCGCCCtttgt                                                                                >  2:520821/1‑225 (MQ=255)
                                                                                                                                                           ttACCTCCTGCCGGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCgtggt                                                              <  2:201832/228‑1 (MQ=255)
                                                                                                                                                            ttctgccggCCAGGGATTTTGGGGGGCGGGGGTATCGATTGTCCTGCCGGAATGTAGAGGGGGGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTGGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCgtggt                                                              <  1:380458/219‑1 (MQ=255)
                                                                                                                                                                  cTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtg                                                       <  2:806/228‑1 (MQ=255)
                                                                                                                                                                       aGGGATTTTTGGTGGCGCGGGGGGTGATTTTCCGGCCGGAATGCAGGGGGGTGAGGTTAAATGGTTGCTGCTGGCAAGAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTggtggt                                                  <  1:104262/228‑1 (MQ=255)
                                                                                                                                                                                        cggcgATCGATATTGCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGt                                 <  2:379966/228‑1 (MQ=255)
                                                                                                                                                                                          cggtgatggtttttCCTGCCGGAATGCAGAGCGGTTGGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTGTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGt                                 <  1:566251/216‑1 (MQ=255)
                                                                                                                                                                                              tcgatATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTTGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTAATTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATg                            <  2:428020/227‑1 (MQ=255)
                                                                                                                                                                                                        ccTGCGGGAATGGAGGGTGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTTGTTTGGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGta                    <  2:608665/227‑1 (MQ=255)
                                                                                                                                                                                                          tGCCGGAATGCAGGGCGGTGATTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGtata                  <  2:304889/227‑1 (MQ=255)
                                                                                                                                                                                                            ccGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGGGGTGTTTTTGTGGGTTGGGTGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCg               <  2:189722/228‑1 (MQ=255)
                                                                                                                                                                                                                 aTGGAGAGTGGTGGGTTTAAATGGTTGCGTCGGGCAACAGGGGGTTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATAtt          <  1:493559/228‑1 (MQ=255)
                                                                                                                                                                                                                      ggggggTGGGTTTAAATGGTTGGTGCTGGCAACAGGGGGGTTTTTGGGGGTTGGTGGCTGGTTGTGCGGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCg     <  1:590194/224‑1 (MQ=255)
                                                                                                                                                                                                                       agCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTCTTTGGGCCGCGGGAAAGGCGTGGGACGGTTTAGTCCTTTTTTTTTCCGCCGCGGGTTGGTGGTGGCTGCGCGCGGGGGTTTTCTCCTCTGGCGCGGGGGGGGGGGTGGGGTTAATTTGCCACCCCTTGATGTCgggtgttagagcgtgtttgtg     >  1:185722/1‑209 (MQ=255)
                                                                                                                                                                                                                         cGGTGAGTTTAAATGGGTGCTGCTGGCAACAGCGGTGGTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTaa  <  2:293514/228‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 
GCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAA  >  NC_000913/71660‑72099

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