Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,947,067 (T)8→7 intergenic (+35/+61) trpT → / ← hdfR tRNA‑Trp/flhDC operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,947,0600T.100.0% 259.5 / NA 54intergenic (+28/+68)trpT/hdfRtRNA‑Trp/flhDC operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (24/30);  total (24/30)

CAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCAGCCATATG  >  NC_000913/3946912‑3947206
                                                                                                                                                    |                                                                                                                                                  
cAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAttt                                                                                                                                                  <  1:628811/151‑1 (MQ=255)
   ggggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttttt                                                                                                                                               <  2:632366/151‑1 (MQ=255)
        cgcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttttttatct                                                                                                                                          >  1:512380/1‑147 (MQ=255)
         gcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTg                                                                                                                                        >  2:462599/1‑151 (MQ=255)
              ttCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAAt                                                                                                                                    >  1:547388/1‑150 (MQ=255)
                cGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAAc                                                                                                                                 <  2:373196/151‑1 (MQ=255)
                  aGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAAccc                                                                                                                               <  2:425893/151‑1 (MQ=255)
                      ccgtccgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCtctc                                                                                                                           <  1:608008/151‑1 (MQ=255)
                        gtccgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCg                                                                                                                         <  2:409577/151‑1 (MQ=255)
                        gtccgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCg                                                                                                                         <  2:128629/151‑1 (MQ=255)
                        gtccgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCg                                                                                                                         <  1:398652/151‑1 (MQ=255)
                           cgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAg                                                                                                                      <  2:547236/151‑1 (MQ=255)
                             ttCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGt                                                                                                                     >  2:496202/1‑150 (MQ=255)
                              tCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTa                                                                                                                    >  2:207659/1‑150 (MQ=255)
                               ccGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTaa                                                                                                                   <  2:536328/150‑1 (MQ=255)
                                    aCCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCtt                                                                                                             <  2:540099/151‑1 (MQ=255)
                                     cccTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCtt                                                                                                             <  2:312540/150‑1 (MQ=255)
                                      ccTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCt                                                                                                           <  2:47935/151‑1 (MQ=255)
                                         aaTTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTAcc                                                                                                        <  1:477395/151‑1 (MQ=255)
                                               gggCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATc                                                                                                  >  2:584672/1‑151 (MQ=255)
                                                        tCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagc                                                                                         >  1:178007/1‑151 (MQ=255)
                                                         cAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagcc                                                                                        <  1:509059/151‑1 (MQ=255)
                                                         cAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagcc                                                                                        >  2:371634/1‑151 (MQ=255)
                                                             tGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTc                                                                                    >  1:300319/1‑151 (MQ=255)
                                                              ggTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCa                                                                                   <  2:228232/151‑1 (MQ=255)
                                                               gTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCat                                                                                  <  2:178002/151‑1 (MQ=255)
                                                                    gCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTaa                                                                             <  2:373760/151‑1 (MQ=255)
                                                                             tcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTc                                                                     >  1:296200/1‑150 (MQ=255)
                                                                                     cGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCAc                                                            <  1:423273/151‑1 (MQ=255)
                                                                                     cGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCAc                                                            >  2:539487/1‑151 (MQ=255)
                                                                                      gggTGTTGGGAGTTAGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCGTGCCAGCCATTCATATTAATACACTTCATCCAGCAc                                                            <  2:39025/150‑1 (MQ=255)
                                                                                          gtTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTa                                                        >  2:618604/1‑150 (MQ=255)
                                                                                                 gTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaa                                                <  2:459177/151‑1 (MQ=255)
                                                                                                   tCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaa                                               <  2:193701/150‑1 (MQ=255)
                                                                                                     gAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaa                                               >  1:500087/1‑148 (MQ=255)
                                                                                                            tcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGa                                      >  2:302733/1‑150 (MQ=255)
                                                                                                             ccGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAAt                                    <  1:264609/151‑1 (MQ=255)
                                                                                                                 cccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAAc                                <  1:637856/151‑1 (MQ=255)
                                                                                                                  ccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACg                               <  2:335109/151‑1 (MQ=255)
                                                                                                                  ccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACg                               >  1:484118/1‑151 (MQ=255)
                                                                                                                     gCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAt                            >  2:300997/1‑151 (MQ=255)
                                                                                                                       cAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAttt                          >  1:327139/1‑151 (MQ=255)
                                                                                                                       cAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAttt                          <  1:646069/151‑1 (MQ=255)
                                                                                                                        aGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAtttt                         >  1:258314/1‑151 (MQ=255)
                                                                                                                          aaaTCATCCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAttttt                        >  1:254485/1‑150 (MQ=255)
                                                                                                                            atcatcCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATc                     <  1:520583/151‑1 (MQ=255)
                                                                                                                            atcatcCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATc                     <  1:257592/151‑1 (MQ=255)
                                                                                                                              catcCTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGc                   >  1:89395/1‑151 (MQ=255)
                                                                                                                                  cTTAGCGAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTAtt               >  1:329276/1‑151 (MQ=255)
                                                                                                                                   ttAGCGAAAGCTAAGGAATTTTTTATCTGAAATAACCCTCTCCGCAGTAAATACTCCTACCGGCATCCTTGGCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAACGATCGCGAATCAACGCATTTTTATCGCTAttt              >  2:649789/1‑151 (MQ=255)
                                                                                                                                       cgcaaGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCATTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGc           >  2:539925/4‑150 (MQ=255)
                                                                                                                                        gAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCAg         <  1:507767/151‑1 (MQ=255)
                                                                                                                                        gAAAGCTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCAg         <  1:216150/151‑1 (MQ=255)
                                                                                                                                             cTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCAGCCat     >  1:178006/1‑150 (MQ=255)
                                                                                                                                             cTAAGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCAGCCat     >  1:401920/1‑150 (MQ=255)
                                                                                                                                               aaGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCAGCCata    <  1:84881/149‑1 (MQ=255)
                                                                                                                                               aaGGATTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCAGCCATATg  <  1:279541/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCAGCCATATG  >  NC_000913/3946912‑3947206

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