Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 71,877 A→G 28.7% 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,8770AG28.7% 117.5 / 9.7 73E176G (GAG→GGG) yabIDedA family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (14/35);  new base G (20/3);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 3.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAG  >  NC_000913/71645‑72115
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
gATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGCAATTCCTGCCGGAATGCGGAGCGGGGGGTTTaa                                                                                                                                                                                                                                               >  1:581658/1‑239 (MQ=255)
       ctgacctggCGTTGCATCAACCCAGCCGGTTCCCCATTCTGGGCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTtg                                                                                                                                                                                                                                        <  1:219131/231‑1 (MQ=255)
       cTGAACATGCGTTGCATTAAAACAGCATGTTCCCCATTCTGTCCGGTCGTTTTTTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTtg                                                                                                                                                                                                                                        <  1:310680/239‑1 (MQ=255)
       cTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCTAAATTTATTACGCCGGATATTATCGGCTGGCTGCTGTGGCCGCCGTTTTACTTCCCGCCAGGGATTCTGGCGGGCGCGGGGGTCGGGATTTCTTCCGGAATGGAGCGCGGGGGGGTTaaaaaggtg                                                                                                                                                                                                                                        >  2:398388/1‑233 (MQ=255)
          aaCATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGGAATTCCTGCCGGAATGCGAGGCGGGGGGTTTAAATGGTtgcgt                                                                                                                                                                                                                                     >  1:957111/1‑237 (MQ=255)
          aaCATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGGGTTTAAATGGTtgct                                                                                                                                                                                                                                      >  1:840812/1‑238 (MQ=255)
           aCATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGCTATTCCTGCCGGAATGCAGACCGGGGAGTTTAAATGGTtgttgc                                                                                                                                                                                                                                    >  1:518468/1‑239 (MQ=255)
               gCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACCCCGAATATTATCGGCTGCCTGCTTTGGCCGCCGTTTTACTTCCTGCCAGGGATTTTGGCGGGCGCGGCGAGCGGTATTCCTGCCGGAATGCAGGGCGGGGGGTTTAAATTGTtgctgggggg                                                                                                                                                                                                                                >  2:1077152/1‑234 (MQ=255)
                  ttGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGGTCTGCCGGTTGCTAAATTTATTACGCCGAATATTATCGGCTGCCCGCTGTGGCCGCCGTTTTACTTCCTGCCCGGGGTTTTGACGGGGGCGGCGCTCGCTATTCCTGCCCGAATGCAGAGCGGGGGGGTTAAAAGGTTGCGGGTGGGaag                                                                                                                                                                                                                             >  2:897071/1‑238 (MQ=255)
                    gCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTTGCGGGCGCGGCGGTCGATATTCCTTCCGGGATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAg                                                                                                                                                                                                                          >  1:229850/1‑239 (MQ=255)
                    gCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGGGCGGCGATCCATATTCCTACCGGAAAGCAAGGCGGGGAGGTTAAATGGTtgtggcggggaaaag                                                                                                                                                                                                                           >  2:922965/1‑226 (MQ=255)
                    ccatcacagcggcgtGTTCCCCTTTGTGGCCGGTTGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAg                                                                                                                                                                                                                          <  1:417835/225‑1 (MQ=255)
                         aacccgacttgtccacatttCTGGCGGGTCTTTTTGTTGGCCCGACGGGTCCGGTGGTGCCAATGGTGGGGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGgtg                                                                                                                                                                                                                    <  1:398388/221‑1 (MQ=255)
                           cacaGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGGGGGGTTAAATGGTTGCTGCTGGGAACAGCGGGGtt                                                                                                                                                                                                                  >  1:385275/1‑239 (MQ=255)
                             caGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGGGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGGGTTTAAATGGTTGCTGCTGGCAACAGCGGGGtttt                                                                                                                                                                                                                >  1:55706/1‑239 (MQ=255)
                               gCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGGATATTATCGGCGGCCTTCTGTGGCCGCCGGTTTACTTCCCGCCAGGGATTCTGGCGGGCGCGGCGCTCGTCATTCCATCCGGACTGCAGAGAGGGGGAGTTTAAAGGGTCCTGCTCGCAAAGGCGGGTGttttg                                                                                                                                                                                                               >  2:389074/1‑238 (MQ=255)
                               gCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTTCCGGAAAGCAGAGCGGGGGGTTTCAAATGTTTCGGCTGGGAACAGCGGGGTTTTtg                                                                                                                                                                                                              >  1:712150/1‑239 (MQ=255)
                               gCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGGTTTACTTCCTGCCCGGGATTCTGGCGGGCGCGGGGACCGGCATGCCTGCCCGAATGGAGGGGGGGGGGGTTAAAAGGGTGCATGTGGGCACAGGGGGGttttgg                                                                                                                                                                                                              >  2:587021/1‑237 (MQ=255)
                                cATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGGGTTCGTAATTCCTGCCGGAATGAAAAGCGGGGGGTTTAAATGGTTGCTGCTGGCAACAGCGGGGGTTTtgt                                                                                                                                                                                                             >  1:359668/1‑239 (MQ=255)
                                cATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGGGCGATATTCCTGCCGGCATGCAGAGCGGTGGGTTTAAATTGTTGCTGCTGGCAACAGCGGGGGTTTtgt                                                                                                                                                                                                             >  1:1071324/1‑239 (MQ=255)
                                cATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCTAAATTTATTACGCGGAATATTATCGGCTGCCTGCTTGGGCCGCCGGTTTATTTCCCGCCCGGGATCTTGGCGGGGCCGGGCGTCGCAATTTCTGCCCGGATTGAGAGGGGGGGGGTTAAAAGGTGGGCGATGGCACAAGGGGGGttttttg                                                                                                                                                                                                             >  2:350604/1‑237 (MQ=255)
                                 attgtccccttgtggggtggtggtTTTGTTGGCCCGACGGGTCCGCTGGGCCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTCCGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTtgt                                                                                                                                                                                                             <  1:1077152/220‑1 (MQ=255)
                                   gTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGGTCGGTATTCCTGCCGGAATGCAGAGCGGGGAGTTTAAATGGTTGCTGCTGGCAACAGCGGGGGTTTTGTggg                                                                                                                                                                                                          >  1:485930/1‑239 (MQ=255)
                                     tCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGCTATTCCTGCCGGAATGGCGAGCGGTGAGTTTAAATGGTTGGTGCTGGGAACAGCGGTGTTTTTGTGGGtt                                                                                                                                                                                                        >  1:194748/1‑239 (MQ=255)
                                     tCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGGGGGTTTAAATGGTGGCGGCTGGCAACACCGGGGTTTTTTTGGGtt                                                                                                                                                                                                        >  1:397860/1‑239 (MQ=255)
                                     tCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCTAAATTTTTTACGCCGAATATTATCGGCTGCCTGCTGGGGGCGCCGTTTTAATTTCTTCCAGGGATTCTTGCGGGGGCCGGGGCTCTTATTCCCTCCGGGAAGCACAGGGGTTGGGTTTAAAGGTTGGCGCTGGGAAAAGCGGGGTTTTTGTggggt                                                                                                                                                                                                        >  2:810075/1‑237 (MQ=255)
                                         cATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCCGAGCGGTGGGTTTAAATGGGTGCTGCTGGCAACCGCGGGGTTTTTGGGGGTtggtg                                                                                                                                                                                                    >  1:989278/1‑239 (MQ=255)
                                          aTTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGGATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGTTATTCCTGCCGGAATGCCGAGCGGTGAGGTTAAATGGTTGTCGCTGGCAACACCGGGCGTTTTGTGGGTtggtgg                                                                                                                                                                                                   >  1:313650/1‑239 (MQ=255)
                                          aTTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATTGGCTGCCTGCTGGGGCCGCCGTTTTACTTCCTGCCCGGGGTTCTGGCGGGCGCGCGGATCGCTATTCCCGCCGGTGTGCGGCACGGGGGGTTTTTATGGTGGTCGCTGGCAACCGGCGGGGTTTTTTGGGGtggggg                                                                                                                                                                                                   >  2:299033/1‑239 (MQ=255)
                                          aTTCTGGTCGGTCGATTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATTCTGGTTCTGTCGGTGGCTAAATTTATTACCCCGGATATTATCGGCTGCCTGCTGTTGCCGCCGTTTTCCTTCCTGCCCGGGCTTTTGGCGGGGCCGGGGAGCCACATTTCTGTCGGGATTCAGATGGGGGTGTTTTAATTGGGTCTTATGGGACAAGCGGTGGTTTTGTGGGTtggggg                                                                                                                                                                                                   >  2:538728/1‑239 (MQ=255)
                                                gtcggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTTTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGGTCGGTATTCCTTCCGGAATGCAGGGCGGGGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGtgtggggggg                                                                                                                                                                                                >  1:965284/1‑227 (MQ=255)
                                                gtcggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCCGCCGGGGGCTAAATTTATTACGCCGGATATTATCGGCTGGCTGCTGTGGCCGCCGGTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGCTCGTTATTCCTGCCGGAATGCATGGCGGGGAGGTTAAATGGTCGCTGCTGGCCAAAGCGGGGTTTTTTGGGGGGGgtgggtgggg                                                                                                                                                                                             >  2:142506/1‑237 (MQ=255)
                                                 tcggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGGTCGGTATTTCTGCCGGGATGCAGGGCGGGGGGGTTAAAAGGTTGCGGGTGGCAACAGGGGGGTTTTTGTgggggggggggtggggg                                                                                                                                                                                            >  2:748423/1‑225 (MQ=255)
                                                     tcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGGTCGGTATTCCTGCCGGAATGCAGAGCGGTGAGGTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTTGGGGGTGGTggctggcttg                                                                                                                                                                                           >  1:196524/1‑234 (MQ=255)
                                                                  ccGGCGCGTCCGCTGGGGCCAATGGTGGCGGAAATGCTGGATCTGCCGGTGGTTAAATTTATTCCGCCTAATATTACCGGCTGCCGGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                                           <  1:1035436/239‑1 (MQ=255)
                                                                  ccGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGGATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCCGGGATTCTGGCGGGGGCGGCGGTCGGGATTCCTGCCGGGAGGCCAGACGGGGGGTTTAAATGGTTTGGGCCGGGCACAGGGGGGTTTTTTTGGGGGGGTGGGGGGGGGTGCTGgggggggtgggg                                                                                                                                                                           >  2:92705/1‑231 (MQ=255)
                                                                   cggcgggtgcGCTGGTGCCAATGGGGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTCCGCCGAATATTAGCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                                           <  1:994438/229‑1 (MQ=255)
                                                                    gccgcgTCCGTTGGTGTCATTGGTGGGGGGAATGCTGGATCTGCCGGTGGCTAAATTTTTTGCGCCGAATTTTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCa                                                                                                                                                                         <  2:551167/237‑1 (MQ=255)
                                                                     aCGCGTCCGCTGGTGCCAATGGTGGCGGGGATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCCCCGTTTTAATTCCTTCCAGGGATTCTGGCGGGCGCGGCGCCCCGCATTCCTCCCCGTATTCACACCGGGGGGTTTAAATGGGTGCTGCTGGGCACCGCGGGGTTTTTGTGGGGGGGGggcttgcttgggtgggggggaggtgggac                                                                                                                                                                        >  2:268551/1‑218 (MQ=255)
                                                                                    ccAGTGGGGGCGGGAATGCGGGTTCTGCCGGTGGCTAATTTTTTTCCGCGAAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGa                                                                                                                                                       <  2:55706/239‑1 (MQ=255)
                                                                                       aTGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGGGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTTTTgcgcgcggcgaaaagagccggccg                                                                                                                                                      >  1:868108/1‑219 (MQ=255)
                                                                                          gtggCGGGATTTGTGGATCTTCTGTTGGCTAAATTTATTCCGCCGATTATTCTCGGCTGCCTGCTGTGGCCGCCGTTTTCCTTCCGGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGttt                                                                                                                                                 <  1:460466/239‑1 (MQ=255)
                                                                                           ggtgggggattgtttgttttGCCGGTGGCTTATTTTATTCCGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTTCTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGttt                                                                                                                                                 <  2:385275/220‑1 (MQ=255)
                                                                                            gggggggatgggtGGATCGGCCGGTGGCTAATTTTTTTCCGCCGAATTTTATCGGTTGCCTGCGGGGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGttt                                                                                                                                                 <  2:989278/226‑1 (MQ=255)
                                                                                            ggggggAATTGCTTGTTCTGCCGGTGGTCTAATTTATTAGGCCGAATATTATCGGTTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCATAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGttt                                                                                                                                                 <  1:810075/235‑1 (MQ=255)
                                                                                            ggCGGGAAGTCGGGATGTGCCGGGGGCTAAATTTATTACGCCGAATATGACCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTg                                                                                                                                                <  1:587021/238‑1 (MQ=255)
                                                                                                             ccgctggTGGCTTATTTTTTTACGCGAAATTTTTCGGGCGGCCTGCGGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCGGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcg                                                                                                                               <  2:581658/234‑1 (MQ=255)
                                                                                                             cTGCCGGTGTCTAAATTTATTCCGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcg                                                                                                                               <  1:748423/239‑1 (MQ=255)
                                                                                                             cTGCCGGTGGCTAAATTTTTTACGCCGAATTTTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcg                                                                                                                               <  1:922965/239‑1 (MQ=255)
                                                                                                               tccgggtgtttatttttttGCGCCAAATTTTTCCGGCTGCCTGGTGTGGCCGCCTTTTTCCTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcgcg                                                                                                                             <  2:485930/226‑1 (MQ=255)
                                                                                                                     gtcTAAATTTATTCCGCCGAATTTTTGCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGt                                                                                                                        <  1:998774/236‑1 (MQ=255)
                                                                                                                               atgaCGCCGAATAGTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCt                                                                                                             <  1:710412/236‑1 (MQ=255)
                                                                                                                                   cGCCGAATAGTCGCGGCTGCCGGGTGGGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACg                                                                                                         <  1:538728/239‑1 (MQ=255)
                                                                                                                                      cgattttcttCGGGTTCCTGCTGTGGCCGCCGTTTTTCTTCCTGGCAGGGATTCTGGCGGGCCGGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCg                                                                                                      <  2:232209/230‑1 (MQ=255)
                                                                                                                                        aaTATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTTTTGGGCA‑‑GCGGGAAAGCGGCTGACAGTTTTAGTGATTATTTTTCCCCCGGGCGGTTTTTTTGGCTGGAGCCgcg                                                                                                    >  1:718086/1‑237 (MQ=255)
                                                                                                                                        aaTATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGGGAGTTTAAATGGTTGCTGCTGGCAACAGCGGGGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGGGGGTGTAGTGCACCGGGGAAAAGCGGTGCCCGTTTGTGCGATATTTTTTCCCCCGGCGGTTTGTtgtgattacggcgcg                                                                                                     >  2:102397/1‑228 (MQ=255)
                                                                                                                                                    tgctggatggggacgccgTTTGACGTCCTGCCAGGGATGCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTTACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCAt                                                                                        <  1:678022/226‑1 (MQ=255)
                                                                                                                                                             gggCCGCCGTTTTTCTGCCTGCGGGGTATTGTGGCGGGTGCGTGTGTGGATATTCCTGCCGGAATGGAGAGCGGTGGGTTTAAATGGTTGCTGCTGGCATCAGCGGTGTTTTTGTGGTTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGAGTTCTGCCATCGGCgtggt                                                                              <  2:693390/238‑1 (MQ=255)
                                                                                                                                                              ggcgctctttttTCTTCTTGCCAGGGGTTGTGGCGGGCGCGGGGGTCGATTTTCCTGCCGGAATGCAGGGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCgtggt                                                                              <  2:313650/230‑1 (MQ=255)
                                                                                                                                                                 cgccgTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgcg                                                                          <  1:859957/239‑1 (MQ=255)
                                                                                                                                                                     tttttACTCCCTGCCGGGTGTTATGGGGGGCGCGGCGTTTGATTTTCCTGCCGGAATGCAGAGCGGTTAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtgg                                                                      <  2:229850/238‑1 (MQ=255)
                                                                                                                                                                            ttCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcgc                                                                                                                                                                          <  1:515103/135‑1 (MQ=255)
                                                                                                                                                                                      aggttCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTGAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCAc                                                     <  1:571812/236‑1 (MQ=255)
                                                                                                                                                                                                gggcgGGGGGTTGGATTTTCCTGCCGGAATGCAGGGCGGGGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGc                                           <  2:359668/239‑1 (MQ=255)
                                                                                                                                                                                                gggcgGGGGGATCGATATTCCTGGCGGAATGCAGGGCGGTGGGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATg                                            <  1:464142/238‑1 (MQ=255)
                                                                                                                                                                                                      cggcgATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGgt                                       <  1:92705/238‑1 (MQ=255)
                                                                                                                                                                                                                 atTCCTGCCGGAATGCAGAGCGGTGAGTTTAGATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGatat                           <  1:614988/239‑1 (MQ=255)
                                                                                                                                                                                                                            aaTGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTaa                  <  1:339043/237‑1 (MQ=255)
                                                                                                                                                                                                                                  caggtcGGTGGGTTTAAATGGTTGCTTCTGGCAACAGCGGTGTTTTTGTTGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGt            <  1:68509/234‑1 (MQ=255)
                                                                                                                                                                                                                                  agagCGGTTAGTTTAAATGGTTGCTGCTGGCAACAGGGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGt            <  1:910487/238‑1 (MQ=255)
                                                                                                                                                                                                                                  agagCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGGGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGt            <  1:350604/238‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTTGTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCTACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTg          <  1:102397/237‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGGGGC‑‑GGGGGGAAGGGACTGACCGGTTGTGTCATTATTTTTTCCGCGGGCGGTTGTTTTGGCTGAGGCCCTTTGTTTTTTCCATCGTGTGGGGGGGGGGGGGGGGGGTTATTTCTCCCCCCGGTTTTCCCGGGGGATATGGAGATTTTGTTTAAAGAgggggggg       >  2:339043/1‑239 (MQ=255)
                                                                                                                                                                                                                                          gttaGTTTAAATGGTTGCTGGTGGCAACAGCGGTGTTTTTGGGGGGTGGTGGCTGGCTGTGCTGGCGGTTATGGCGGA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGt      <  1:204615/235‑1 (MQ=255)
                                                                                                                                                                                                                                             aGTTTAAATGGTTGCTGCTGGCAACAGGGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAg  <  1:458020/239‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
GATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCA‑‑GCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAG  >  NC_000913/71645‑72115

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