Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,766,151 G→T 10.2% R161L (CGA→CTA) ‡ PP_5599 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,766,1510GT10.2% 179.6 / 7.3 58R161L (CGA→CTA) ‡PP_5599LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (28/24);  new base T (2/4);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

AAGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGAATGACAATGACCTG  >  NC_002947/3766005‑3766296
                                                                                                                                                  |                                                                                                                                                  
aaGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCga                                                                                                                                                   <  5:326742/148‑1 (MQ=255)
aaGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCga                                                                                                                                                   <  2:304837/148‑1 (MQ=255)
        ccccGCAACTGAGGCTACGCATCAGGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCTCGAAGAAcct                                                                                                                                          <  7:220593/149‑2 (MQ=255)
             cAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCTCGAAGAAcctgtata                                                                                                                                     <  8:1072/149‑7 (MQ=255)
                cTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCt                                                                                                                                  >  2:256935/1‑149 (MQ=255)
                   aGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGagaa                                                                                                                                                >  5:284198/1‑132 (MQ=255)
                   aGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGagaa                                                                                                                                                <  6:284198/132‑1 (MQ=255)
                   aGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCtttt                                                                                                                               <  6:232411/149‑1 (MQ=255)
                     gCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGa                                                                                                                             >  2:48703/1‑149 (MQ=255)
                      cTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGAc                                                                                                                            <  1:193603/149‑1 (MQ=255)
                          gCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCa                                                                                                                          >  8:7211/1‑147 (MQ=255)
                              cACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGAcc                                                                                                                    >  7:306184/1‑149 (MQ=255)
                                  ccAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCTCGAAGAAcct                                                                                                                                          >  2:233424/1‑122 (MQ=255)
                                  ccAGGGGACTCGCGCAGTATTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCTCGAAGAAcct                                                                                                                                          <  1:233424/123‑2 (MQ=255)
                                   cAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGtt                                                                                                               <  1:256935/149‑1 (MQ=255)
                                          cTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAgc                                                                                                       <  7:102745/149‑1 (MQ=255)
                                                  gTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATc                                                                                               <  2:314223/149‑1 (MQ=255)
                                                  gTCTTTACGAACAAGTAGCCTGCGGAGAGCTAGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATc                                                                                               >  1:33158/1‑149 (MQ=255)
                                                        aCGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGATCAAAGACCCCGTTGACAAGCGCGACATCCGCTc                                                                                          <  1:231554/148‑1 (MQ=255)
                                                        aCGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTc                                                                                          <  3:155932/148‑1 (MQ=255)
                                                          gAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCTCGAAGAAcctgt                                                                                                                                        >  8:341938/1‑98 (MQ=255)
                                                          gAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCTCGAAGAAcctgt                                                                                                                                        <  7:341938/101‑4 (MQ=255)
                                                          gAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTcgac                                                                                       >  2:257216/1‑149 (MQ=255)
                                                          gAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGACCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTcgac                                                                                       >  1:86032/1‑149 (MQ=255)
                                                                gTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCt                                                                                 <  7:7211/149‑1 (MQ=255)
                                                                    ccTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGt                                                                             <  8:306184/149‑1 (MQ=255)
                                                                     cTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCaaa                                                                                                                        >  5:79391/1‑106 (MQ=255)
                                                                     cTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCaaa                                                                                                                        <  6:79391/106‑1 (MQ=255)
                                                                         ggAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGaaa                                                                        >  8:90357/1‑149 (MQ=255)
                                                                           agagCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAAcc                                                                      <  2:33158/149‑1 (MQ=255)
                                                                                   aTGCGGTAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCg                                                               <  1:257216/149‑1 (MQ=255)
                                                                                      cGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGata                                                            <  2:40957/148‑1 (MQ=255)
                                                                                            tCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTca                                                     >  3:310489/1‑149 (MQ=255)
                                                                                               tGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGaaa                                                                        >  1:165982/1‑127 (MQ=255)
                                                                                               tGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGaaa                                                                        <  2:165982/127‑1 (MQ=255)
                                                                                                  tATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGcc                                               >  2:44462/1‑149 (MQ=255)
                                                                                                      tCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTggg                                           <  1:44462/149‑1 (MQ=255)
                                                                                                        cGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTgggg                                          >  1:145096/1‑148 (MQ=255)
                                                                                                             ttCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTggg                                      >  5:81068/1‑147 (MQ=255)
                                                                                                                      cccAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGccc                            <  2:86032/148‑1 (MQ=255)
                                                                                                                       ccAAAAGCATTCGTATACCGGTTCTTCGAGAAGAACCACTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGAATGcc                             >  6:313408/1‑146 (MQ=255)
                                                                                                                       ccAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAACCCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCcag                          >  8:352241/1‑149 (MQ=255)
                                                                                                                       ccAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCcag                          >  5:250363/1‑149 (MQ=255)
                                                                                                                       ccAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCcag                          <  1:48703/149‑1 (MQ=255)
                                                                                                                       ccAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCcag                          >  3:128009/1‑149 (MQ=255)
                                                                                                                       ccAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCcag                          >  7:206099/1‑149 (MQ=255)
                                                                                                                       ccAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGAATCGATATGACTCACATGCCTGGGGTGGGCAGATTGcc                             >  3:353369/1‑146 (MQ=255)
                                                                                                                         aaaaGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGaaa                                                                        <  3:278580/101‑1 (MQ=255)
                                                                                                                         aaaaGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGaaa                                                                        >  4:278580/1‑101 (MQ=255)
                                                                                                                               aTTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACtt                  >  3:64221/1‑149 (MQ=255)
                                                                                                                               aTTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACtt                  <  4:310489/149‑1 (MQ=255)
                                                                                                                               aTTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGAATCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACtt                  <  4:353369/149‑1 (MQ=255)
                                                                                                                                   gTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaa               <  6:81068/148‑1 (MQ=255)
                                                                                                                                     ataCAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaatga            >  2:62740/1‑149 (MQ=255)
                                                                                                                                       aCAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaatgac           >  4:266079/1‑148 (MQ=255)
                                                                                                                                       aCAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGAGACATCCGCTAGACGATTATTGGTCGAAACCTGATTCGATATGACTCAAATGCCTGGGGTGCGCAGATTGCACAGCAGTACTTGaattac           >  8:265773/1‑148 (MQ=255)
                                                                                                                                         aGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaa               >  3:17585/1‑142 (MQ=255)
                                                                                                                                         aGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaa               <  4:17585/142‑1 (MQ=255)
                                                                                                                                           gTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaatgacaatga      >  3:50012/1‑149 (MQ=255)
                                                                                                                                                 cGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTCACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACCTGAATGACAATGACCTg  <  3:332967/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGAATGACAATGACCTG  >  NC_002947/3766005‑3766296

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