Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,648,455 C→G 8.3% R1278P (CGG→CCG)  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,4550CG 8.3% 182.1 / 7.5 60R1278P (CGG→CCG) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base G (2/3);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

ATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAAC  >  NC_000913/2648156‑2648751
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                         
atatCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTAACCCg                                                                                                                                                                                                                                                                                                          >  2:556795/1‑301 (MQ=255)
 tatCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATCGCGCAGCAGACGCTCATTACCCCgg                                                                                                                                                                                                                                                                                                         >  2:505000/1‑301 (MQ=255)
 tatCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGAGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAACAGCGCAGCAAACGCTCATTACCCCgg                                                                                                                                                                                                                                                                                                         >  1:188129/1‑301 (MQ=255)
  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:35987/301‑1 (MQ=255)
                      cggcgTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACg                                                                                                                                                                                                                                                                                   <  2:574249/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)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                         
ATATCCGCTCATCACTATTACGCGGCGTTTTCAGCGCCAGCGCAATCGCTTCTTCACCACGCGTCGCATCACCCATGGTTTTCAGCGCAACGCCAAGTTGCAGCAGCGGTAAACCAGAAGCGGCATCTGCACGATGCTCCCAGATTTCACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAAC  >  NC_000913/2648156‑2648751

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: