Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,524,184 C→T A438A (GCC→GCT)  mrcA → peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,524,1840CT98.6% 241.5 / NA 73A438A (GCC→GCT) mrcApeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (27/44);  minor base G (1/0);  total (28/45)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGG  >  NC_000913/3523939‑3524430
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
gcCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTAt                                                                                                                                                                                                                                                   <  1:102432/251‑1 (MQ=255)
 ccAATCCGCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATg                                                                                                                                                                                                                                                  <  2:378300/251‑1 (MQ=255)
  cAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCC‑TGGTGGCTGGCCCAAGTGCCGGAAGTGAACTCGGCGCTGGGGTCGATCAATCCGCAAAACGGTGCTGTTATgg                                                                                                                                                                                                                                                 >  2:329001/1‑251 (MQ=255)
         cAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGt                                                                                                                                                                                                                                          <  1:106766/251‑1 (MQ=255)
            cAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGGGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATGCTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAAGCCGCAAAACGGTGCTGTTATGGCGCTGGTCgt                                                                                                                                                                                                                                       <  2:218017/251‑2 (MQ=255)
             aaGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAACGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGACCTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGt                                                                                                                                                                                                                                      >  1:354418/1‑251 (MQ=255)
                        aTGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGAttt                                                                                                                                                                                                                           <  1:316894/251‑1 (MQ=255)
                                      gTCGACAGGCGCATTGGGTAAGAAAGGCGTTCTCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAGGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGTCGATGCA‑TGGGGGCTGGAACAAGTGCCGAAAGAGCAATCGGCGCAGGTGTCGATCAATCCGCAAAAAGGTGCTGGTATTGCGCTGGTCGGGGGCGTTGAATTCCATCatggcatgt                                                                                                                                                                                                             >  2:136577/1‑243 (MQ=255)
                                                  aTTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTCCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGGCAGGTTGGCGATGCA‑TGGTGGCTGTCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCGGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCa                                                                                                                                                                                                 <  1:36079/251‑1 (MQ=255)
                                                           ggAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCAc                                                                                                                                                                                        >  1:348206/1‑251 (MQ=255)
                                                            gAAGGCGTTCGCTGGGCGCGTCCTTACAGTTCGGATACTCAGCAAGGACCGACGCCGCGTACAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACt                                                                                                                                                                                       >  2:253761/1‑251 (MQ=255)
                                                             aaGGCGTTCGCGGGGCGCTGCCTTACCGTTCGGATACTCAGCAAGGCCCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGGCAGGTTGGAGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCGGCCCAGGCACTg                                                                                                                                                                                      <  2:112643/251‑1 (MQ=255)
                                                                          ggCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGtt                                                                                                                                                                         >  2:134594/1‑251 (MQ=255)
                                                                               gTCCTTACCGTTCGGATACACAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTGTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAAc                                                                                                                                                                    <  2:43730/251‑1 (MQ=255)
                                                                                    tACCGTTCTGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCACATCTGGGTTCGTCAGATTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGACGTGAACTCGGCGCTGGGGTCGATCAATCCGCAAAACGGTGCGGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCa                                                                                                                                                                >  2:176693/1‑250 (MQ=255)
                                                                                    tACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCACCATCa                                                                                                                                                                >  2:12027/1‑250 (MQ=255)
                                                                                    tACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATc                                                                                                                                                                 <  2:391383/249‑1 (MQ=255)
                                                                                       cGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCAACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAAcc                                                                                                                                                            >  2:240735/1‑251 (MQ=255)
                                                                                                                  ccgcGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCa                                                                                                                                                                                                 >  2:105429/1‑187 (MQ=255)
                                                                                                                  ccgcGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCa                                                                                                                                                                                                 <  1:105429/187‑1 (MQ=255)
                                                                                                                  ccgcGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATaa                                                                                                                                 <  1:283723/251‑1 (MQ=255)
                                                                                                                   cgcgTAAAGTGACCGCTGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATaa                                                                                                                                 >  2:335726/1‑250 (MQ=255)
                                                                                                                   cgcgTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATaa                                                                                                                                 >  2:209429/1‑250 (MQ=255)
                                                                                                                    gcgTAAAGTGACCTATTTTCTGCAACAGGGTCATCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TTGTGGCTGGCACAAGTGCCGTAAGTGAACTCGGCGCTGGTGTCGATCAAGCCTCAAATCGGTGGTGTTATGCCGCTGGTCGGTGGCTTTGATTACAATCAGGGCAAGTTTAACCGCTCCCCCCAGGCTCTTCGTCAGGTGGGTTCCAACATCAAACCGTTCCTTTACACCCCGGCGATGGATaaa                                                                                                                                <  2:226968/250‑1 (MQ=255)
                                                                                                                     cgTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAgg                                                                                                                              <  1:12027/251‑1 (MQ=255)
                                                                                                                      gTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGCTGGa                                                                                                                                    >  2:83787/1‑244 (MQ=255)
                                                                                                                      gTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGa                                                                                                                                    <  1:83787/244‑1 (MQ=255)
                                                                                                                       tAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTc                                                                                                                            >  1:270306/1‑251 (MQ=255)
                                                                                                                             gACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGc                                                                                                                      >  1:239210/1‑251 (MQ=255)
                                                                                                                             gACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAAGCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACg                                                                                                                       <  1:329001/250‑1 (MQ=255)
                                                                                                                                cGATGCTCTGAAAACGGGTCAGGAAATCAGGGCTCGCAAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTAGGCGCTGGTGTCGAACAATCCGCAAAACGGTGATCTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGACACACAAGCAATGCGTCAGGTGGGTTCCAACATCAAACCGTATCTCTACCACGCGGCGATGGACACAGGTCTGAGGCTgg                                                                                                                   >  2:279848/1‑251 (MQ=255)
                                                                                                                                    gTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAg                                                                                                               <  2:313124/251‑1 (MQ=255)
                                                                                                                                      tcggCAAACGGGGCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGGGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGGGCTGTTATGGCGCGGGTCGGTGGCTGTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGGCAGGTGGGTGCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAATGTCTGACGCTGGCAAGt                                                                                                              <  2:201949/247‑1 (MQ=255)
                                                                                                                                       cTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTAt                                                                                                            >  2:31406/1‑251 (MQ=255)
                                                                                                                                          cAAACGGGTCAGCAAATCTGGTTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCTCTGCGTCAGGTGGGTTCCAACATCAAACCGTGCCTCTACACCGCGGCGATGGATAAAGGTCTGTCGCTGGCAAGTATGtt                                                                                                         <  2:403609/251‑1 (MQ=255)
                                                                                                                                           aaaCGGGTAGGCAAATCGGGGTTCTGCAGGTTGGCGATGCAGTGG‑GGCTGGCACGAGTGCGGAAAGTGAACTCGGCGCTGGTGTCGATCAAGCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTGTGAGTTCAATCAGAGCAAGTTTAACCGTGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAACCCGTGCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGtt                                                                                                         <  2:36165/250‑1 (MQ=255)
                                                                                                                                                 gTCAGCAAATCTGGGTACGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCAGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGATGGCAAGTATGTTGAACg                                                                                                    >  2:126080/1‑249 (MQ=255)
                                                                                                                                                    aGCAAATCTGGGTTCGTCAGGTTGGGGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGAtgtg                                                                                               <  2:79661/251‑1 (MQ=255)
                                                                                                                                                       aactCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAAGCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGTCAAGTATGTTGAACGATGTGCCa                                                                                            <  2:24653/248‑1 (MQ=255)
                                                                                                                                                         aTCTGGGTTAGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAAt                                                                                          >  1:325924/1‑251 (MQ=255)
                                                                                                                                                           cTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAAttt                                                                                        >  1:113609/1‑251 (MQ=255)
                                                                                                                                                             gggTTCGCCGGGGTGGCTATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTTGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGCACGATCTGCCAATTtct                                                                                      <  2:335287/251‑1 (MQ=255)
                                                                                                                                                               gTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCg                                                                                    <  1:80670/251‑1 (MQ=255)
                                                                                                                                                                   gTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTGCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTgg                                                                                <  2:224944/251‑1 (MQ=255)
                                                                                                                                                                      aGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGAt                                                                             <  2:348206/251‑1 (MQ=255)
                                                                                                                                                                             cGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTg                                                                      <  1:31406/251‑1 (MQ=255)
                                                                                                                                                                              gATGCA‑TGGTGGCCGGCACAAGTGCCGGACGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGATGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGc                                                                     >  1:128350/1‑251 (MQ=255)
                                                                                                                                                                               aTGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGcc                                                                    <  1:3321/251‑1 (MQ=255)
                                                                                                                                                                                 gCA‑TGTTGGCCGCCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGGTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGGTCCAACATCAAACCTTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCGGGCAAGTATGTTGAACGATGTGCCAATTTCGCTCGGGGATGCAAGTGCCg                                                                   <  2:354418/250‑1 (MQ=255)
                                                                                                                                                                                   a‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGtt                                                                <  1:176693/251‑1 (MQ=255)
                                                                                                                                                                                       gtggCTGGCACAAGTGGCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTg                                                             <  2:77725/251‑1 (MQ=255)
                                                                                                                                                                                          gctggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACt                                                          <  1:63666/251‑1 (MQ=255)
                                                                                                                                                                                             ggcACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTg                                                         <  1:358458/249‑1 (MQ=255)
                                                                                                                                                                                                acaAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCGGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGc                                                    <  2:325924/251‑1 (MQ=255)
                                                                                                                                                                                                    gTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCgaa                                                <  2:321191/251‑1 (MQ=255)
                                                                                                                                                                                                     tGCCTGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTg                                                                                                                          <  1:29541/176‑1 (MQ=255)
                                                                                                                                                                                                     tGCCGGAAGTGAACTCGGCGCTGGTGTAGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTg                                                                                                                          >  2:29541/1‑176 (MQ=255)
                                                                                                                                                                                                      gcctgAAGTGAACTCGGCGCTGGTGTCGATCAAACCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGCTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATGAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTTCCAATTTCTCGCTGGGATGCAATTGCCGGTTCTGACTGTCAGCCgaaga                                              <  1:371279/247‑1 (MQ=255)
                                                                                                                                                                                                       ccGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCACTTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCgaaga                                              >  1:299663/1‑250 (MQ=255)
                                                                                                                                                                                                               gAACTCGGCTCTGGTGTCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACAGCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTTACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTcacca                                      <  1:50500/250‑1 (MQ=255)
                                                                                                                                                                                                                              gtCGATCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTc                      <  1:209429/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCCACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAAtt                 >  1:74271/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tCAATCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAAtt                 >  1:45226/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aaTCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCg               <  1:98488/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aaTCCGCAAAACGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGGGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCg               >  2:101581/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aaaCGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTc       <  1:319345/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aaCGTTGCTGTTATGGTGCTGGTCGGTTGCTTTGATTTCCAGCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTCCACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGGTGTGCCAATTTCGCGCTGGGATTCAAGTCCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGTAGTATGCTGGTCCAATTCGCTTACGTCa      <  2:326896/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aaCGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTc       >  1:15153/1‑250 (MQ=255)
                                                                                                                                                                                                                                              aaCGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCa      <  1:134594/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aCGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTATGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTAACCACCGCAGTATGCTGGTCCAATTCGCTTAAGTc       >  1:398754/1‑249 (MQ=255)
                                                                                                                                                                                                                                                cGGTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAgg    <  1:216818/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAgggg  <  1:240735/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTGCTGTTATGGCGCTGGTCGGTGGCTTTGATTGCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAgggg  <  1:395837/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
GCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCA‑TGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCGTTCCTCTACACCGCGGCGATGGATAAAGGTCTGACGCTGGCAAGTATGTTGAACGATGTGCCAATTTCTCGCTGGGATGCAAGTGCCGGTTCTGACTGGCAGCCGAAGAACTCACCACCGCAGTATGCTGGTCCAATTCGCTTACGTCAGGGG  >  NC_000913/3523939‑3524430

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