Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	644,176	C→T	intergenic (‑209/+21)	rnk ← / ← rna	regulator of nucleoside diphosphate kinase/ribonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	644,176	0	C	T	100.0%	225.2 / NA	59	intergenic (‑209/+21)	rnk/rna	regulator of nucleoside diphosphate kinase/ribonuclease I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/31);  total (28/31)

CAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGT  >  NC_000913/644035‑644322                                                                                                                                              |                                                                                                                                                   cAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTg                                                                                                                                            <  1:553774/150‑1 (MQ=255) cAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTg                                                                                                                                            <  1:660505/150‑1 (MQ=255)   ggcggcGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTaa                                                                                                                                         <  1:24136/151‑1 (MQ=255)   ggcggcGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTa                                                                                                                                          >  1:90075/1‑150 (MQ=255)       gcgAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCa                                                                                                                                     <  1:763386/151‑1 (MQ=255)           gggTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAgatg                                                                                                                                 <  2:822276/151‑1 (MQ=255)             gTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAgatgat                                                                                                                               >  1:695023/1‑151 (MQ=255)               ttAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGAtt                                                                                                                              <  2:192666/150‑1 (MQ=255)                 aaCTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATtaat                                                                                                                           >  2:161019/1‑151 (MQ=255)                 aaCTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGGTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATtaa                                                                                                                            <  2:474026/150‑1 (MQ=255)                   cTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATtaata                                                                                                                          >  2:303097/1‑150 (MQ=255)                     cacaTGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATtaataa                                                                                                                         <  1:164713/149‑1 (MQ=255)                         tgatgaTACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCt                                                                                                                   >  2:635314/1‑151 (MQ=255)                             gatACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTat                                                                                                               <  2:787403/151‑1 (MQ=255)                                actgactgTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTatca                                                                                                             <  2:646275/150‑1 (MQ=255)                                 ctgactgTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTatcaat                                                                                                           >  1:272933/1‑151 (MQ=255)                                   gactgTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTatcaatca                                                                                                         <  1:13141/151‑1 (MQ=255)                                   gactgTTGCGCACTCTTTGACGTGATTTCCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTagcaa                                                                                                            >  2:469350/1‑148 (MQ=255)                                    actgTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATcac                                                                                                        <  2:632044/151‑1 (MQ=255)                                      tgTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATcaca                                                                                                       >  2:311253/1‑150 (MQ=255)                                         tGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTGTAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCGCAAAgg                                                                                                   >  1:377942/1‑151 (MQ=255)                                         tGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAg                                                                                                    <  1:386803/150‑1 (MQ=255)                                         tGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAg                                                                                                    <  2:62571/150‑1 (MQ=255)                                          gCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAgg                                                                                                   <  1:355570/150‑1 (MQ=255)                                               ctctTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGc                                                                                             >  1:190805/1‑151 (MQ=255)                                                   ttGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGcc                                                                                            <  2:637203/148‑1 (MQ=255)                                                     gAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAg                                                                                        >  1:571230/1‑150 (MQ=255)                                                      aaGTGATTTGCGTCACATTCAGGGAATACCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGtt                                                                                      >  1:766951/1‑151 (MQ=255)                                                              tGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGt                                                                              >  2:406549/1‑151 (MQ=255)                                                                   cacaTTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGt                                                                         <  1:811081/151‑1 (MQ=255)                                                                      aTTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTtgt                                                                      >  1:275429/1‑151 (MQ=255)                                                                        tCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTgg                                                                    >  1:635417/1‑151 (MQ=255)                                                                         cAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTgg                                                                    <  1:718308/150‑1 (MQ=255)                                                                          aGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGc                                                                   >  2:767461/1‑150 (MQ=255)                                                                             gAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGa                                                               <  1:793628/151‑1 (MQ=255)                                                                              aaTTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGaa                                                              <  1:27535/151‑1 (MQ=255)                                                                               aTTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAAt                                                             >  2:860897/1‑151 (MQ=255)                                                                                     cAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGt                                                         <  2:714696/149‑1 (MQ=255)                                                                                      aaTGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTg                                                      >  1:206661/1‑151 (MQ=255)                                                                                           aGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGa                                                  >  2:403362/1‑150 (MQ=255)                                                                                              gcATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGc                                              <  2:84829/151‑1 (MQ=255)                                                                                               cATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCg                                             <  1:416339/151‑1 (MQ=255)                                                                                                aTTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCgg                                            <  2:431167/151‑1 (MQ=255)                                                                                                   aTGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGc                                         >  1:410221/1‑151 (MQ=255)                                                                                                      tataAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGtt                                      >  1:165293/1‑151 (MQ=255)                                                                                                          aaTCTTAATCGCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATg                                  <  2:698183/151‑1 (MQ=255)                                                                                                                    gCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCtt                        >  1:105877/1‑151 (MQ=255)                                                                                                                    gCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCtt                        <  1:311279/151‑1 (MQ=255)                                                                                                                    gCCTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCtt                        >  2:695076/1‑151 (MQ=255)                                                                                                                     ccTTGGTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCttt                       <  2:174836/151‑1 (MQ=255)                                                                                                                         ggTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATc                   <  1:821362/151‑1 (MQ=255)                                                                                                                         ggTTTATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATc                   >  1:187773/1‑151 (MQ=255)                                                                                                                            ttATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCga                 <  1:508865/150‑1 (MQ=255)                                                                                                                             tATGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGTGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCg                  >  1:536505/1‑148 (MQ=255)                                                                                                                               tGGAAGACGAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGGTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCg                  >  2:51822/1‑146 (MQ=255)                                                                                                                                       gAATTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTc     >  1:419849/1‑151 (MQ=255)                                                                                                                                         aTTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTc     >  1:261980/1‑149 (MQ=255)                                                                                                                                         aTTAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAg   <  1:436553/151‑1 (MQ=255)                                                                                                                                           tAGTGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGt  <  2:197748/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   CAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGT  >  NC_000913/644035‑644322

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: ‑