Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,648,458 T→G 8.1% N1277T (AAC→ACC)  yfhM ← bacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,648,4580TG 8.1% 156.4 / 5.6 59N1277T (AAC→ACC) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base G (2/3);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

ATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGA  >  NC_000913/2648158‑2648758
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                             
atCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCACAGAAGACGCTCATCACCCCGGt                                                                                                                                                                                                                                                                                                               >  1:35936/1‑301 (MQ=255)
       catcaCTATTACGCGGCGTTTTCAGCGCCAGCGCACTCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTTACGGGCCAACACCAGCGCAGCGCAAGACTGTACGGCGAATTTACTGGCCTTGAGATTATCCGCGTCCGGGATCGCCATACTGCCCGGATCTTGACAATCACGCAGCAGACACTCCTTACCCCGCTTCATGGc                                                                                                                                                                                                                                                                                                       >  1:247542/1‑301 (MQ=255)
       catcaCTATTACACGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCGGCGGTAAACCAGAAGCGGCATCTGCACGATGCCCCCAGATTTCACGCAGTGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCGGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGATATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTACATAGCGCAGCAGACACTCATTACCCCGGTTAATggg                                                                                                                                                                                                                                                                                                       >  2:512324/1‑300 (MQ=255)
             tattaggggtggTTGTCTGCGCCTGCGCAATCGCTTCGTCACCACGCGTCGCATCACCCAGGGTTTGCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTGTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                                                                                                                                                                                 <  1:112326/290‑2 (MQ=255)
                tACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgctgc                                                                                                                                                                                                                                                                                              <  2:624296/301‑5 (MQ=255)
                    cggcgTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACg                                                                                                                                                                                                                                                                                          <  2:574249/301‑1 (MQ=255)
                    cggcgTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACg                                                                                                                                                                                                                                                                                          <  2:35987/301‑1 (MQ=255)
                        tgtttCGGCGCCAGCGCAATCGCTTCTTCACCCCGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTGTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGt                                                                                                                                                                                                                                                                                      <  1:238671/299‑1 (MQ=255)
                                  ccagcgcAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGACCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGttt                                                                                                                                                                                                                                                                            >  2:65067/1‑300 (MQ=255)
                                   ccgcgcGATCGCGTCTTCACCATGCTGCGCATCACCCAGTGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCg                                                                                                                                                                                                                                                                           <  1:64999/299‑1 (MQ=255)
                                             gCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGa                                                                                                                                                                                                                                                                 >  2:229641/1‑301 (MQ=255)
                                                 tttgacctcGCGTCGCATCACCCTTGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                                                                                                                                                                             <  2:291804/297‑1 (MQ=255)
                                                          gcgTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCAt                                                                                                                                                                                                                                                    <  2:84473/301‑1 (MQ=255)
                                                                     cccATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCg                                                                                                                                                                                                                                         >  2:469546/1‑301 (MQ=255)
                                                                                cAGCGCAACGCCACGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCCCGATGCTCCCAGCTTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGACTTTGAGATTATCCGCGTACGGGATCGAAATCATGCCCGGAACTTGTAAATAGCGCAGCAGACGCCCATTACCCCGGTAAATGGCCTATGTCGGCCCCCTCTAACCCTCATCGCCTACCCCGACCAAGaacacaaacacgcacacccacaaccatcac                                                                                                                                                                                                                              >  1:33517/1‑273 (MQ=255)
                                                                                   cgcAACGCCACGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTCNCGCGCATTTACTGGCTTTCCGATTATCCCCGTACGGGATCGCNNNNATGCCCGGCTCNNGTACATAGCCCCGCANNNCCTCATTACNNNNNTCAATGGCATCTGTCCGCAAGCTGTAAACCCTTTCGCCTGCGCGAACCCGAACACCCAACACCTACACCATCACACAGTActct                                                                                                                                                                                                                           >  1:402200/1‑301 (MQ=255)
                                                                                          ccAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAc                                                                                                                                                                                                                    >  2:197299/1‑301 (MQ=255)
                                                                                                tGCAGCAGCGGTAAACCAGAATCGGCAGTTGCACGATGCTCCCGGATTTCACGCAGCGCACCCAGCGTAGCCTTTTGCTGACTGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAAGTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGtttt                                                                                                                                                                                                              <  1:456754/301‑1 (MQ=255)
                                                                                                      agcGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATccc                                                                                                                                                                                                        >  2:629021/1‑301 (MQ=255)
                                                                                                             aaCCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGAGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGCCCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGCCACCGTTTTTATCCCACAgcgc                                                                                                                                                                                                 >  1:392244/1‑301 (MQ=255)
                                                                                                                        gCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGt                                                                                                                                                                                      <  2:247809/301‑1 (MQ=255)
                                                                                                                                 ggATGCTCCCGGATTTCCCGCAGGGCACCCAGCGGTGCGTTTGGCGGACGGGCCAACACCGTCGTAGCGTAAGAGTGTACGGCGAATTGACTGGCTTTGAGATTATCCGCGTACGGGAACGACATGATGCCCGGATCTGGTAAATAGCGCAGCAGACGATCATGACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGtt                                                                                                                                                                             <  1:504328/300‑1 (MQ=255)
                                                                                                                                    tGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCa                                                                                                                                                                          <  2:143573/301‑1 (MQ=255)
                                                                                                                                      cTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAtt                                                                                                                                                                        >  2:632673/1‑301 (MQ=255)
                                                                                                                                        cccAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTAACACGTTGCAtttt                                                                                                                                                                      >  1:305054/1‑300 (MQ=255)
                                                                                                                                         ccgggttttACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGGGATTATCGGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCACCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc                                                                                                                                                                     <  1:100070/296‑1 (MQ=255)
                                                                                                                                         ccAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc                                                                                                                                                                     >  2:132533/1‑301 (MQ=255)
                                                                                                                                                 cACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAACTCCATCCCGTAAGCCGTCAGCCCGTACTCTTCGTCCCCGATTTTATCCCACAGCGCAAAGCCGCCGTTATCCCGTTGCCTTTGCAGCAAAcc                                                                                                                                                             >  1:379948/1‑300 (MQ=255)
                                                                                                                                                             ccggtggggcgttttGCTGAGGGGCGAACGCCAGCGCAGGGTAAGACTGTACGGGGAATTTACTGGGTTGGAGATTAGCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGat                                                                                                                                                 <  1:584238/290‑1 (MQ=255)
                                                                                                                                                                                   gggcaacctaagcggagggTGAGACTGTACGGGGAATTGTCTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAAGAGCGCAGCAGACGCTCATTACGCCGGTTAATGGCGTCTGTCGGCAAGCTGTAACCCTGTGCGCCTGCGCGGACCAGGAAATCCATCAGGTAAGCCGTCAGCCAGTACTCGTCGTCACCGTTTTTAACCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTtc                                                                                                                           <  1:511652/291‑1 (MQ=255)
                                                                                                                                                                                   ggCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTtc                                                                                                                           <  2:305458/301‑1 (MQ=255)
                                                                                                                                                                                    gCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGATCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACATCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATAACCACCGATGCACGGCGattat                                                                                                                          >  2:216372/1‑299 (MQ=255)
                                                                                                                                                                                    ccccacacCAGTGAGGCGTAAGTCTGTACGGCGAATTTACTGGCTGTGAGATTATCCGCGTACGGGATCGACATCATGCGCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTtct                                                                                                                          <  1:197056/297‑1 (MQ=255)
                                                                                                                                                                                      cacctccggggcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGGACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCtca                                                                                                                        <  1:631796/291‑1 (MQ=255)
                                                                                                                                                                                      cAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGAGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCtca                                                                                                                        <  2:392758/301‑1 (MQ=255)
                                                                                                                                                                                         cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCtcatca                                                                                                                     >  2:432451/1‑301 (MQ=255)
                                                                                                                                                                                          cccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCAc                                                                                                                    <  2:288841/300‑1 (MQ=255)
                                                                                                                                                                                            cggggcggcgtaggcttttACGGCGAATTTACTGGCTTTGTGATTATCCGCGTAGGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTg                                                                                                                  <  1:556081/283‑1 (MQ=255)
                                                                                                                                                                                                                 ggCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTCTCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCAACACTGACCCCTTTGAAGCCCAAcccc                                                                                             >  1:342667/1‑301 (MQ=255)
                                                                                                                                                                                                                      ttttACTGTCTTTGAGATTATCCGCGTACGGGATTGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAg                                                                                        <  2:187646/300‑1 (MQ=255)
                                                                                                                                                                                                                           gttgCTTGGTGATTATCCGGGTAGGGGATCGACAGCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTggg                                                                                   <  1:363632/298‑1 (MQ=255)
                                                                                                                                                                                                                           cTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAActgcaggcgttggg                                                                                   <  2:188371/301‑14 (MQ=255)
                                                                                                                                                                                                                           cTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTggg                                                                                   >  2:546430/1‑301 (MQ=255)
                                                                                                                                                                                                                             ggCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCt                                                                                 >  2:421114/1‑301 (MQ=255)
                                                                                                                                                                                                                                     atttgtCCGCTTACGGGTTGGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAActgcaggcgttgggcatcaaaggc                                                                         <  1:431887/296‑24 (MQ=255)
                                                                                                                                                                                                                                      tttATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGta                                                                        <  2:178802/300‑1 (MQ=255)
                                                                                                                                                                                                                                                tACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAACCGGGAAATGCCGATATCGACCGCTGCCCGGCGTTTCTCATCCCTGTCGCCTTTGCTGCCCAACCCCCGCAGTTCGGCTGCATTCGtcaaacatcacg                                                              >  1:364720/1‑290 (MQ=255)
                                                                                                                                                                                                                                                tACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACg                                                              >  2:458974/1‑301 (MQ=255)
                                                                                                                                                                                                                                                       tCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCCCATCCCTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAg                                                       >  1:18467/1‑301 (MQ=255)
                                                                                                                                                                                                                                                        cGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAga                                                      >  2:141911/1‑300 (MQ=255)
                                                                                                                                                                                                                                                        cGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGACTTTGATACCCAACGCCTGCAGTTGGGCTGCGATGGTATAAAATGAAGGAAACAgg                                                      >  1:122792/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                tGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGc                                              <  2:380444/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        gcggACGCTCATTACCCCGGTTAATGGCTTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACgg                      <  1:545757/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        gcagACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACgg                      >  2:627956/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          agACGCCATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGaa                    >  2:243740/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                            gacgccATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATa                   >  2:210760/6‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                   ttACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGTCTGAGCGGACCATGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTa           <  2:402714/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   ttAACCGGGGTAATGAGCGTCTGTCGGCCCGTTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCCCGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCtttt            <  1:243473/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     aacccGGTTAATGGCGTCTGTCGGCACGCTGTAACCGTGTTGCGCTGGGCGGACCAGGAAATCCGTCGCGTAGGCCGTCAGCCAGTACACTTCGTCACCGTTTTTAGCCCACAGAGAAAAGCCGCCGTTATAACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTGGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAAc         <  1:229398/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                           tttAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTg   <  2:55108/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                           gTTAGTGGCGTCTGTCGGCCCGCTGTAACCCTGTTCGCTTGGGCGGACCAGGAAATCCACCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTGTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGAACGGATACGCTTTTAACTCTTTg   <  1:468945/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                            ttAATGGGGTCGGTCGGCATGCTGTAACCGTGTTCGCCTGCGCGGACCAGGAAATCCATCGCGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTGGCAGCAAACGGGAAATGCCGAGATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGa  <  1:359875/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                             
ATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGA  >  NC_000913/2648158‑2648758

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: