Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,827 C→T D44D (GAC→GAT creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,636,8270CT100.0% 104.2 / NA 36D44D (GAC→GATcreCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/22);  total (14/22)

CGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGC  >  NC_000913/4636591‑4637068
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
tgATACCCACATTAAAACGTTCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTg                                                                                                                                                                                                                                    >  1:74948/2‑251 (MQ=255)
            tAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGCTGTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTgg                                                                                                                                                                                                                          <  2:67123/250‑1 (MQ=255)
               aaCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGcgcg                                                                                                                                                                                                                      >  1:245171/1‑251 (MQ=255)
                      cGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTTGTACTGGCTATTTTTGTCAAAGAAGTTTAAACGGGCGTGCGAAGGGCAACGGAGGGGGCGGTAAATGATACCGCCACGTTGCTGGGGGGGCTGGggcgggccgct                                                                                                                                                                                                               >  2:364210/1‑244 (MQ=255)
                                 ccgcgcCATCAACCCCGATCTTTCACCGTTTAATACTCATCGCGGCATGGGATACAGCCTGAGGGGCCTGTAATGCGTAGCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTGCGTACTGGCTATTGTTGTCAAAGAAGTTAAACCGGTCATGCTAAGAGCAACGGATGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGctctct                                                                                                                                                                                                      <  1:216340/248‑1 (MQ=255)
                                          aaaacccGATCTGTCACCTAGTAATTCTATTCGCGGCATGGGATATAGCCTGAGGGGCCTTGAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTGTTACTGGTTGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGGCAAAGAAGGTAAACCGGGCGGGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAAc                                                                                                                                                                                           <  2:302516/247‑1 (MQ=255)
                                                        cACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCTTGTAATGCGTATCGGCATGCGGTTGCTGCTTGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTg                                                                                                                                                                              <  1:379396/250‑1 (MQ=255)
                                                                   gttgaTCGCGGTTTGGGATATAGTCTGGGGTGGCTGTAATGCGTATCGGCAGTCGGTTGCTGCCGGGATAGTTTTGACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGGCCAAGAAGTTAAACCGGGCGGGCGAAGAGCAACTGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCGGGCGTGTCCCGAGTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGAGCTGGCGttt                                                                                                                                                                  <  2:207666/247‑1 (MQ=255)
                                                                         cgcgGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTTCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATc                                                                                                                                                              >  1:373166/1‑249 (MQ=255)
                                                                            ggCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTa                                                                                                                                                         <  1:364210/251‑1 (MQ=255)
                                                                             gCATGGGATATATCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGTCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTac                                                                                                                                                        <  1:194777/251‑1 (MQ=255)
                                                                             gCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAAATGGCGCAGGCGTTTAATCAGcaac                                                                                                                                                        >  1:109948/1‑248 (MQ=255)
                                                                                         gCCTTAGGTGCCTGTTATGCGTGTCGGCATGCGGTTGTTGCTGGGCTATTTTTGACTGGTGGCGGTGGCGGCCTGGGTTGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGATCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGTGCTGGCGCGTCCCGATTTGCTCTCTGGGTACCCACCGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGt                                                                                                                                            <  2:9772/251‑1 (MQ=255)
                                                                                               ggggCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCTGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTcgcg                                                                                                                                      <  1:209500/251‑1 (MQ=255)
                                                                                               ggggCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCCACGCAAGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTcgcc                                                                                                                                      >  1:16459/1‑250 (MQ=255)
                                                                                               ggggCCTGTAAGGCGTACCGGCATGCGGTTGCTGCTTGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCGGGCGGAGCTGGCGCGGCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTcgcg                                                                                                                                      <  1:52470/251‑1 (MQ=255)
                                                                                                     tGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACt                                                                                                                                                                               <  1:223895/204‑1 (MQ=255)
                                                                                                        aaTGCGTATCGGCATGCGGTTGCTGCTGGGCTAGTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCg                                                                                                                             <  1:323345/251‑1 (MQ=255)
                                                                                                                       catgttgctGCTTGGCTATTATTTACGGGAGGCGGTGGCGGCCTGGTTCGTACTGGCTATTGTTGTCAAAGAAGTTAAACCGGTCGTGCGAAGAGCAACGGAGGGGACGTTAAGCGATACCTCAACGTTGCTGGCGGAGCGGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTGCGCGCCAATATCGGTGGCATTAACaaa                                                                                                               <  2:282371/247‑1 (MQ=255)
                                                                                                                         ggTTGCTGCTGGGCTCTTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGGACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATAGGTGGCATTAACaaa                                                                                                               >  1:282371/1‑248 (MQ=255)
                                                                                                                                     gCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAatat                                                                                                <  1:118215/251‑1 (MQ=255)
                                                                                                                                        aTTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAt                                                                                             <  1:114126/251‑1 (MQ=255)
                                                                                                                                                                 cTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCATTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAgggt                                                                     >  1:3618/1‑249 (MQ=255)
                                                                                                                                                                    gTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTacaaca                                                                                                                                                    <  1:177844/168‑1 (MQ=255)
                                                                                                                                                                          aCTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGt                                                           <  1:253636/251‑1 (MQ=255)
                                                                                                                                                                           ctggctctttttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGTCGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATTTCGGTGGCCTTACAAAAGTACGCAATGAGAATACGGTCTATATGACCGATTCACAGGGTAAAGTATTGtt                                                          >  2:179178/8‑251 (MQ=255)
                                                                                                                                                                             ggctaTTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAGCAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGtt                                                          <  1:367582/245‑1 (MQ=255)
                                                                                                                                                                                 aTTTTTGTCAAAGACGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCACCGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGAATTAACAAAGGGAGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGAATTGTTCGAttt                                                    >  2:47743/1‑250 (MQ=255)
                                                                                                                                                                                                   aaaCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTgg                                  >  1:75929/1‑251 (MQ=255)
                                                                                                                                                                                                        gggCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAgg                             >  1:58162/1‑251 (MQ=255)
                                                                                                                                                                                                                      cAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTgg               >  1:52661/1‑251 (MQ=255)
                                                                                                                                                                                                                        aCGGAGGGGACGTTAATTGATACCGCAACGGTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGGCCTAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATTGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGaa             <  2:279658/251‑1 (MQ=255)
                                                                                                                                                                                                                              gggACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCTGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGt       <  1:121789/251‑1 (MQ=255)
                                                                                                                                                                                                                               ggACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCg                                                                                <  1:264696/177‑1 (MQ=255)
                                                                                                                                                                                                                               ggACGTTAATCGATACCGCAACGTTGATGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCCTTAACCAAGTGCGCAATGAATATCATGTCTATATGACCg                                                                                >  2:264696/1‑177 (MQ=255)
                                                                                                                                                                                                                                   gttaaTCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGc  <  1:176522/247‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
CGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGC  >  NC_000913/4636591‑4637068

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