Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,521 T→C L204L (TTG→CTG)  radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,5210TC98.7% 227.2 / ‑4.7 75L204L (TTG→CTG) radADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (27/46);  total (29/46)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

GCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCA  >  NC_000913/4626274‑4626752
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
gcaaaTGAAAACGCTGTATGTCCCCGGCGAAGAGTCGCTGCAACAGGTGGCAAGGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGa                                                                                                                                                                                                                                         <  2:275106/247‑1 (MQ=255)
   gATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGacgc                                                                                                                                                                                                                                      >  1:406871/1‑249 (MQ=255)
                 attccACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAAcg                                                                                                                                                                                                                       <  2:253632/247‑1 (MQ=255)
                 aTGTCACCGTCGACGAGTCGCTGCAACAGGTAGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGAATGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGCCTCGATCCAGGTGATGCATATGGCGGCTGTAACGTCATCGCCTGGACGCGTGGCGCAGGTGCGTGAAACGGGGGCTTATCTGACGGCGTTCGCCAAAAcc                                                                                                                                                                                                                       >  2:431334/1‑250 (MQ=255)
                 aTGTCACCGGAGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGCCAATCTCAATATGTTGTCGGAAACCAGCCTCGAACAGCTATGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGAGATGCATATAGCGGATGGACAATCTACACACGAGAGCGTGGCGCAGGTGCGTGAAACAGCGGGTTATTTGAACCGCTTTGCCaacaac                                                                                                                                                                                                                       >  2:607631/1‑247 (MQ=255)
                    tCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGAGATGCATATGGCGGAAGTACAGGCATCGCCAGGCAGCATGGCGCCGGGGCGTGAAACGGCGGCTTATCTGAAGCGCTTAGCCCAAAagcgc                                                                                                                                                                                                                    >  2:560386/1‑251 (MQ=255)
                     cACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAAcgcgc                                                                                                                                                                                                                    <  2:104884/250‑1 (MQ=255)
                        cGGCGACGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGgtg                                                                                                                                                                                                                <  2:149632/251‑1 (MQ=255)
                        cGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGACTGCCGACTGACAATCTCAATATGTTGTCTGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGgtg                                                                                                                                                                                                                <  1:396613/251‑1 (MQ=255)
                           cGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTgg                                                                                                                                                                                                             <  1:147855/251‑1 (MQ=255)
                             aaGAGTCGCTGCAACAGGTTGCTATTCGCGCTCATCGCTTTGGCCTGCCGACTGACAATCTCAAGATGTTGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGGCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCg                                                                                                                                                                                                           <  2:406088/251‑1 (MQ=255)
                              agagTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGa                                                                                                                                                                                                          <  1:565282/251‑1 (MQ=255)
                                     cTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCCGGGGCGTGAAACGGCGGCTTATCTGCCGCGCTTCGCCAAAACGCGCGGTGTGGCGGTTGTCAt                                                                                                                                                                                                  >  1:418580/1‑251 (MQ=255)
                                      tGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGGAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCGTCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGGCGCGCTTCGCCACAACTCGCGGTGTGGCGATTGGCAtg                                                                                                                                                                                                 >  1:117198/1‑251 (MQ=255)
                                         aaCAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGGCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggtt                                                                                                                                                                                              <  2:663535/251‑2 (MQ=255)
                                          aCAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggttg                                                                                                                                                                                             <  1:196819/251‑3 (MQ=255)
                                             gggggCAATGCGCGCTAATCGCCGTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACGCGATCCAGGTGATGCATTTGGCGGATGTACAGTCATCGACTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTAGCTGACGCGCTTCGCCAAAACGCGCGGTGTGTCGATTGTCAtgtttggac                                                                                                                                                                                          <  2:357696/248‑8 (MQ=255)
                                             ggTGGCAAAGCGCGCTCAACTCCTTGGCCTGCCGACGGCAAATCGCAATAGGTTGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggttggac                                                                                                                                                                                          <  2:184053/251‑6 (MQ=255)
                                              gTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGac                                                                                                                                                                                                                                        <  1:565705/206‑1 (MQ=255)
                                                 gCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGt                                                                                                                                                                                      <  1:557203/251‑1 (MQ=255)
                                                   aaTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTaa                                                                                                                                                                                    <  1:363720/251‑1 (MQ=255)
                                                   aaTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTa                                                                                                                                                                                     >  1:198512/1‑250 (MQ=255)
                                                        gcgcTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGAGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGGGGCGGTTGTCAtggttggaaccgtaaccaa                                                                                                                                                                                >  1:351375/1‑235 (MQ=255)
                                                              aTCGCCTTGGCCTGCAGACTGACAATCTTAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCTTGATTGCCGAAGAAGAGCAACTGAAGCTGATGGTAATTGACTCGATCCAGGTGATTCATATGGCTGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGTTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATggc                                                                                                                                                                         <  1:606464/251‑2 (MQ=255)
                                                               tCGCCATGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATAAGGCGGATGTTCAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGGTTATATGCCGCGCATCGCCAAAACGCCCGGGGCGGGGATTGTCAtgggtgggaacgaaacaaaaaaagggc                                                                                                                                                                        >  2:580576/1‑227 (MQ=255)
                                                                       gCCTTCCGACTGATAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCGGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGCTCCAGGTGATGCATATGGCGGATTTACAGTCATCGCCTGGCAGCGTGGCGGAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAATATGGCTCgctggca                                                                                                                                                                <  2:308959/251‑2 (MQ=255)
                                                                        ccTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCgctggcgg                                                                                                                                                               <  1:360040/251‑3 (MQ=255)
                                                                          tgcagACTAACAAGCTCAATATGTTGTCGGATACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGACGAGCAAACGAAGATGATGGAAATTGACTCGATCCAGGTGATGCATATGGGGGATGTACAGTCATCGCCTGGCAGCGTGGCGCTGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGATGGGGCGATTGTCCTAGTTGGACCCGTAACCAACGATGGGtggctggcgggg                                                                                                                                                             >  2:355585/5‑247 (MQ=255)
                                                                            ccGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGccc                                                                                                                                                          >  1:558923/1‑251 (MQ=255)
                                                                                  gACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAAACAAAGATGGCTCGCTGGCGGGCCCTAAAGt                                                                                                                                                    >  1:476342/1‑251 (MQ=255)
                                                                                          cAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTg                                                                                                                                                <  1:492655/247‑1 (MQ=255)
                                                                                              atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACt                                                                                                                                         <  1:431334/250‑1 (MQ=255)
                                                                                              atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACCCTg                                                                                                                                        >  1:578348/1‑251 (MQ=255)
                                                                                                gttgtCTGAAACCAGCATCGAGCAGATCTGCCGGATGGCCGAAGAAGAGCATCCGAAGCGGATGTGAATTGACTCGAGCCAGGTGATGCATATGGCGGATGTACAGGTATCGCCTGGCATCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGATGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTa                                                                                                                                      <  2:328378/251‑1 (MQ=255)
                                                                                                      ggAAACCAGCATCGAGCAGATATGCATGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGGGGCGCAAGGGCGTGAAACAGCGGCATATCTGACGCGCTTCGCCAAAACGCGCGGGGGGGCGATTGTCAAGGGTGGGCACGTAAACAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGAc                                                                                                                                 >  1:365655/1‑250 (MQ=255)
                                                                                                              gCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTg                                                                                                                        <  1:105339/251‑1 (MQ=255)
                                                                                                                   gagcAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCAGATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCtttt                                                                                                                   <  1:26883/248‑1 (MQ=255)
                                                                                                                             gCCTGATTGCCGAAGAAGAGCAACCGATGCTGATAGTATTTGTCTAGATTCTGTTGCTGCATATGGCGTATGGTCAGTCTGCGCCTGGCAGCGTTGCGCCGGTGCGTGTAACGGCGTCTTATCTGACTCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTATCCAAAGATGGTTCGCTGGCGGGGTATGAAGTGCTGGAACACTGTATCTACTGTTCGGTGCttttttctga                                                                                                             <  2:95665/247‑7 (MQ=255)
                                                                                                                             gCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGGGATGCATATGGCGGATGTACAGGCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGc                                                                                                            <  2:371001/248‑1 (MQ=255)
                                                                                                                                tGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGAGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCCGGAACACTGTATCGCCTGTTCCGTGCTTTTGGATGGCGATGcc                                                                                                      >  1:285309/1‑251 (MQ=255)
                                                                                                                                        gaagaaGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCg                                                                                              >  1:286672/1‑251 (MQ=255)
                                                                                                                                          agaagaTCAGCCGAAGTTGGTGGTAATTGACTCGTGGCAGGTGATGCATATGGCGTATGTACTGGCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAACCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGta                                                                                            <  2:223582/251‑2 (MQ=255)
                                                                                                                                          agaagaGCAGCCGAAGCTGTTGGTAATTGACGCGTTCCAGGTTAGGCATATGGCGGAGGTCTAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGTCGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCGCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGGTCGGTGCTTTTGGATGGCGATGCCGTCTCCCGtt                                                                                            <  2:6351/251‑1 (MQ=255)
                                                                                                                                             agagTAACCGAAGCGGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATTTTCAGTCAGCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGTCGCGCGTCGCCAAAACGCGCGGTGTGGCTATTGGCATGGTTGGACACGTAACCAAATATGGCTCGCTGGCGGGCACTAAAGTGCGGGAACACTGTACCGACTGTGCGGTGCTGTTGGATGGCGACGCTGACTCCAGTTTTc                                                                                         <  2:388471/251‑1 (MQ=255)
                                                                                                                                                  aaCCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCAct                                                                                    <  1:214730/251‑2 (MQ=255)
                                                                                                                                                     cGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAATATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCActtta                                                                                 <  1:394882/251‑5 (MQ=255)
                                                                                                                                                      gAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCActttac                                                                                <  1:633276/251‑6 (MQ=255)
                                                                                                                                                            gATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGc                                                                           <  1:505656/250‑1 (MQ=255)
                                                                                                                                                              gggtaatggacgtgATCCAGGTGATGCATATGGCGGAGGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGGGTGAAACGGTGGCTTATCTGACGCGCTACGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCGGTCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCAt                                                                        <  2:155283/238‑1 (MQ=255)
                                                                                                                                                                  aatggatgcgagccAGGTGATGCAAATGGCGGATGGACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCGGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGTCTCCCGTTTTCGCACGTTAGGCAGCCATaaaa                                                                    <  2:32440/239‑1 (MQ=255)
                                                                                                                                                                   aTTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATaaaaa                                                                   <  1:509558/251‑1 (MQ=255)
                                                                                                                                                                   aTTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAAAACTGTATCGACTGTTCGGGGCTTTTGGATGGCGGTGCCGACTCCCGTTTTCGCACTTTACGCAGCCATaaaa                                                                    >  1:508403/1‑250 (MQ=255)
                                                                                                                                                                     tGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAAc                                                                  <  1:277807/250‑1 (MQ=255)
                                                                                                                                                                         ccGATCCGGGTGATGCAAATGGCTGATTTACAGTCATCGCCGGGCGGCGTGGCGCAGGTGCGTGAAACGGCGGTTTATCTGACGCGTGTCGTCAAAACGCGCGGTGTGGCGATTTGCATGGTTGGACACGTAACCAAAGCTGGTCGGCGGGCGTGGTCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCtt                                                             <  2:376233/250‑1 (MQ=255)
                                                                                                                                                                               cAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACCCTGTATCGACTGTTCGGGGCTTTTGGATGGCGGTGCCGACTCCCGTTTTCGCACTTTACGCCGCCCTAAAAACCGCTTCGgcgg                                                       >  1:10935/1‑250 (MQ=255)
                                                                                                                                                                                  gggaggCATATGGTGTATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAACCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGGCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCGGGAACACTGTATCGACTGTTCGGGGCTTTTGGATGGCGATGCCGACTCCCGTGTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGt                                                    <  2:654811/246‑1 (MQ=255)
                                                                                                                                                                                   tGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTTTCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCCTGGTTGGACACGTAACCCAAGATGGCTCGCTTGCGGGCCCTAAAGTGCTGGAACACAGTATTCACTGTTCGGGGCTTTTGGATGGCGATGCCGCCTCCCGTTTTCGCACTTTTCGCCGCAATTAAAACCGCTTctgcacgggg                                                   >  2:99834/1‑242 (MQ=255)
                                                                                                                                                                                    gATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTg                                                   <  1:202487/250‑1 (MQ=255)
                                                                                                                                                                                    gATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGa                                                  <  1:6086/251‑1 (MQ=255)
                                                                                                                                                                                         atatGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGAAACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGTCTGTTCGGTGCTTTTGGATGGCGATGCCGAATCCCGTTTTTGCACTTTAAGCAGCCCTAAAAAACGCTTCGGCGAGGTGAATg                                               >  1:300109/1‑249 (MQ=255)
                                                                                                                                                                                              gCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGCTGGCTCGCTGGCGGGCCCTAAAGTGCTGGCAAACTGTATCGGCTGTTCGGTGCTTTTGGATGGCGCTGCCGACTCCCGTTTTCGCACTTTACGCAGGCATAAAAAACGCTTCGGCGCGGTGAATGAGCTggg                                        >  1:410471/1‑251 (MQ=255)
                                                                                                                                                                                               cGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACTGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATTGCTCgctggc                                                                                                                                                                 <  2:446261/130‑1 (MQ=255)
                                                                                                                                                                                               cGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCCAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCgctggc                                                                                                                                                                 >  1:446261/1‑130 (MQ=255)
                                                                                                                                                                                                  aTGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTc                                    >  1:72568/1‑251 (MQ=255)
                                                                                                                                                                                                  aTGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTc                                    >  1:527580/1‑251 (MQ=255)
                                                                                                                                                                                                     tACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTc                                 >  1:492594/1‑251 (MQ=255)
                                                                                                                                                                                                      aCAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTc                                 <  1:119060/250‑1 (MQ=255)
                                                                                                                                                                                                              cGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGCTGCCGACTCCCGTTTTCGCACTTTACGCCGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGAcc                        >  1:221082/1‑251 (MQ=255)
                                                                                                                                                                                                               gCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTTGACACGTAACCAAAGATAGCTCGCTGGCGGGACCTAAAGTGCTGGAACACTGTATCGGATGTTCGGTGCTTTTGGAAGGCGCGGCCGCCTCCCGTTTTCGCCCTTTAAGCCGCCATAAAAACCGCTTCGGCGCGGGGCAAGGGCTGGGCGTCTTcactagtaccc                       >  2:458050/1‑241 (MQ=255)
                                                                                                                                                                                                               gCCTGGCAGCGGGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGACATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGGACACTGGATCGACTGTTCGGTGGTTTTGGATGGCGATGCCGACCCCCGTTTTCGCCCTTTACCCCGCAAAAAAAACACCTTCGGCGCGGTGAATGAGCAGGGCGTCTTccgcatggccc                       >  2:348741/1‑241 (MQ=255)
                                                                                                                                                                                                                ccTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGGACACTGTATCGACTGGTCGGGGCTTTTGGATGGCAATGCCGACTCCCGTTTTCGCACTTTACGCAACCATAAAAACCGCTTCGGCGCGGGGGATGAGGTGGGCGTCTTCGCGAttgccg                       >  2:545849/1‑245 (MQ=255)
                                                                                                                                                                                                                   ggcagggggtggagggggCGTGAAACGGCGGCTTATCTGACGAGATTCCGCAAAACGCGTGGTGTGGGGATGTGAATGGTTGGACACGTAACCAAAGATGTCTCGCTGGTGGGACCTAAAGTGCTGGAAAAATGTATCGACTGTTCTGTGCTTTTTGAAGGCGATGGCGACTCCCGTTTGCGTACTTTACGCAGCAATAAAACCCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGCCCGaaaa                   <  2:265528/234‑3 (MQ=255)
                                                                                                                                                                                                                    gcaGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAg                  <  1:120830/249‑1 (MQ=255)
                                                                                                                                                                                                                                  gTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTGGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGt    <  1:296828/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCa  <  1:545849/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
GCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCA  >  NC_000913/4626274‑4626752

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