Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,604,202 G→A *242* (TAG→TAA
E15K (GAA→AAA) 
yjjQ →
bglJ →
putative transcriptional regulator
bgl operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,2020GA100.0% 222.7 / NA 68*242* (TAG→TAA
E15K (GAA→AAA) 
yjjQ
bglJ
putative transcriptional regulator
bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/39);  total (29/39)

TTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCG  >  NC_000913/4604069‑4604317
                                                                                                                                     |                                                                                                                   
ttCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAaaaaaa                                                                                                                >  1:1282713/1‑139 (MQ=255)
      gCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                          <  1:109483/139‑1 (MQ=255)
      gCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                          <  2:120426/139‑1 (MQ=255)
      gCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                          >  1:945231/1‑139 (MQ=255)
           gCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                     <  2:270729/139‑1 (MQ=255)
           gCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                     <  1:1373969/139‑1 (MQ=255)
                 aTGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                               <  1:238993/139‑1 (MQ=255)
                   gTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCg                                                                                             <  2:522568/139‑1 (MQ=255)
                   gTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCg                                                                                             <  2:1282713/139‑1 (MQ=255)
                      tGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                          <  2:825402/123‑1 (MQ=255)
                      tGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                          >  1:825402/1‑123 (MQ=255)
                         aTGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                       <  2:254030/139‑1 (MQ=255)
                                  tAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                               >  1:1167381/1‑106 (MQ=255)
                                  tAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                               <  2:1167381/106‑1 (MQ=255)
                                     cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                >  1:1103181/1‑118 (MQ=255)
                                     cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                <  2:1103181/118‑1 (MQ=255)
                                     cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCa                                                                           >  1:902405/1‑139 (MQ=255)
                                      gcgcAGGGTGCAAGAGAACTTAACTTTCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                    <  2:972203/113‑1 (MQ=255)
                                      gcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGAAAGACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                    >  1:972203/1‑113 (MQ=255)
                                       cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGTCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                         >  2:152233/1‑139 (MQ=255)
                                       cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                         >  2:498988/1‑139 (MQ=255)
                                       cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                         >  1:311020/1‑139 (MQ=255)
                                       cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGa                                                                         >  1:374502/1‑139 (MQ=255)
                                                   gagaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGt                                                             <  1:498988/139‑1 (MQ=255)
                                                     gaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTaa                                                           <  1:184875/139‑1 (MQ=255)
                                                        cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATc                                                                                              <  1:225057/101‑1 (MQ=255)
                                                        cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATc                                                                                              >  2:225057/1‑101 (MQ=255)
                                                          taacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                               <  1:705537/98‑1 (MQ=255)
                                                          taacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                               >  2:705537/1‑98 (MQ=255)
                                                                gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGc                                                >  2:1419038/1‑139 (MQ=255)
                                                                  cTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                              <  2:311020/139‑1 (MQ=255)
                                                                  cTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                              <  1:11963/139‑1 (MQ=255)
                                                                   ttCTGTTTATGGAGATGTCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                             <  1:152233/139‑1 (MQ=255)
                                                                   ttCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                             <  1:1123181/139‑1 (MQ=255)
                                                                     cTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCAt                                           <  2:945231/139‑1 (MQ=255)
                                                                       gTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATAc                                         <  1:775877/139‑1 (MQ=255)
                                                                       gTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTAAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATAc                                         <  1:724506/139‑1 (MQ=255)
                                                                             ggAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATAc                                   <  1:113621/139‑1 (MQ=255)
                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCgg                                                               >  1:793656/1‑106 (MQ=255)
                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCgg                                                               <  2:793656/106‑1 (MQ=255)
                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTtat                                                    >  1:233168/1‑117 (MQ=255)
                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTtat                                                    >  2:76374/1‑117 (MQ=255)
                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTtat                                                    <  2:233168/117‑1 (MQ=255)
                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTtat                                                    <  1:76374/117‑1 (MQ=255)
                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCACTc                              >  1:976754/1‑139 (MQ=255)
                                                                                    ccGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAg                            >  1:1469452/1‑139 (MQ=255)
                                                                                        aGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGt                        <  1:666016/139‑1 (MQ=255)
                                                                                        aGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGt                        <  1:1419038/139‑1 (MQ=255)
                                                                                          aaTGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGa                          <  2:246864/135‑1 (MQ=255)
                                                                                          aaTGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGa                          >  1:246864/1‑135 (MQ=255)
                                                                                                    aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                          >  1:333510/1‑45 (MQ=255)
                                                                                                    aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                          <  2:333510/45‑1 (MQ=255)
                                                                                                    aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATaa                                                  >  2:14788/1‑101 (MQ=255)
                                                                                                    aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATaa                                                  <  1:14788/101‑1 (MQ=255)
                                                                                                    aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATg            <  2:65169/139‑1 (MQ=255)
                                                                                                      ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGcc          <  2:902405/139‑1 (MQ=255)
                                                                                                         aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGcc          >  1:1433270/1‑136 (MQ=255)
                                                                                                         aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGcc          <  2:1433270/136‑1 (MQ=255)
                                                                                                         aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCCATACCAGTCAGGAGTCATTTCAGGATGCCATg       >  2:247796/1‑139 (MQ=255)
                                                                                                              agaagaGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGa                          >  2:631408/1‑115 (MQ=255)
                                                                                                              agaagaGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGa                          <  1:631408/115‑1 (MQ=255)
                                                                                                              agaagaGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTcgcg  >  2:334669/1‑139 (MQ=255)
                                                                                                              agaagaGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTcgcg  >  2:937550/1‑139 (MQ=255)
                                                                                                                             cACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                             >  2:369396/1‑81 (MQ=255)
                                                                                                                             cACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                             <  1:369396/81‑1 (MQ=255)
                                                                                                                             cACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTcg    >  1:1453006/1‑122 (MQ=255)
                                                                                                                             cACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTcg    <  2:1453006/122‑1 (MQ=255)
                                                                                                                                 cATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                              >  2:359856/1‑76 (MQ=255)
                                                                                                                                 cATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                              <  1:359856/76‑1 (MQ=255)
                                                                                                                                     |                                                                                                                   
TTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCG  >  NC_000913/4604069‑4604317

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