Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,148,691	G→A	*357* (TAG→TAA)	plsX →	putative phosphate acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,148,691	0	G	A	100.0%	230.0 / NA	70	*357* (TAG→TAA)	plsX	putative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/43);  total (27/43)

GCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTT  >  NC_000913/1148560‑1148823                                                                                                                                    |                                                                                                                                     gCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACg                                                                                                                               <  2:963363/139‑1 (MQ=255)              aaCAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACt                                                                                                                  <  2:169994/139‑1 (MQ=255)                cAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCg                                                                                                                <  1:1136662/139‑1 (MQ=255)                cAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCg                                                                                                                <  2:879360/139‑1 (MQ=255)                cAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCg                                                                                                                <  2:1215130/139‑1 (MQ=255)                   gCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAg                                                                                                             >  1:651187/1‑139 (MQ=255)                   gCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAg                                                                                                             >  1:832010/1‑139 (MQ=255)                         cAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCa                                                                                                       <  1:739772/139‑1 (MQ=255)                         cAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCa                                                                                                       >  2:859805/1‑139 (MQ=255)                          aGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCaa                                                                                                      <  2:1335256/139‑1 (MQ=255)                             cGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTg                                                                                                   >  2:346517/1‑139 (MQ=255)                               gTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGCACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGAc                                                                                                 >  2:1142225/1‑139 (MQ=255)                                 gCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGttt                                                                                                          >  2:454721/1‑128 (MQ=255)                                 gCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGttt                                                                                                          <  1:454721/128‑1 (MQ=255)                                 gCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAg                                                                                                     <  2:1386951/133‑1 (MQ=255)                                 gCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAg                                                                                                     >  1:1386951/1‑133 (MQ=255)                                             tCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGaa                                                                                   <  1:559775/139‑1 (MQ=255)                                             tCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGaa                                                                                   <  1:1142225/139‑1 (MQ=255)                                               cTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGtata                                                                                           >  2:21577/1‑129 (MQ=255)                                               cTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGtata                                                                                           <  1:21577/129‑1 (MQ=255)                                                          gCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCg                                                                      <  1:346517/139‑1 (MQ=255)                                                          gCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCg                                                                      <  1:141236/139‑1 (MQ=255)                                                          gCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCg                                                                      <  1:563499/139‑1 (MQ=255)                                                              cTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACa                                                                  <  2:1486890/139‑1 (MQ=255)                                                                cGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATg                                                                <  1:1548743/139‑1 (MQ=255)                                                                      gAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACg                                                          <  1:85840/139‑1 (MQ=255)                                                                        aTCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGaa                                                        <  2:832010/139‑1 (MQ=255)                                                                        aTCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCACAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGaa                                                        <  2:981817/139‑1 (MQ=255)                                                                                 cccAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGcc                                                                                                                          <  2:1824778/64‑1 (MQ=255)                                                                                 cccAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGcc                                                                                                                          >  1:1824778/1‑64 (MQ=255)                                                                                 cccAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAg                                                                                                     <  1:117242/85‑1 (MQ=255)                                                                                 cccAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAg                                                                                                     >  2:117242/1‑85 (MQ=255)                                                                                 cccAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACt                                                                           >  1:26488/1‑111 (MQ=255)                                                                                 cccAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACt                                                                           <  2:26488/111‑1 (MQ=255)                                                                                    aGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGTTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTAc                                            <  2:1645337/139‑1 (MQ=255)                                                                                    aGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTg                                                >  1:573993/1‑135 (MQ=255)                                                                                    aGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTg                                                <  2:573993/135‑1 (MQ=255)                                                                                    aGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTAc                                            <  2:68766/139‑1 (MQ=255)                                                                                    aGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTAc                                            >  2:1446720/1‑139 (MQ=255)                                                                                     gCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAg                                                                                                                        >  2:1790779/1‑62 (MQ=255)                                                                                     gCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAg                                                                                                                        <  1:1790779/62‑1 (MQ=255)                                                                                     gCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAGGTGACTGAGCGTa                                                                    <  1:1705053/114‑1 (MQ=255)                                                                                     gCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAGGTGACTGAGCGTa                                                                    >  2:1705053/1‑114 (MQ=255)                                                                                     gCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGAc                                                                            <  2:1347494/106‑1 (MQ=255)                                                                                     gCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACt                                           >  1:31968/1‑139 (MQ=255)                                                                                     gCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTCCAAGTGACGGTATATAACCGAAAAGTGAc                                                                            >  1:1347494/1‑106 (MQ=255)                                                                                      cTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTAc                                            <  2:299774/137‑1 (MQ=255)                                                                                      cTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTAc                                            >  1:299774/1‑137 (MQ=255)                                                                                       tGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCg                                                                                                                <  1:206924/68‑1 (MQ=255)                                                                                       tGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCg                                                                                                                >  2:206924/1‑68 (MQ=255)                                                                                       tGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCa                                                                                                              <  2:768141/70‑1 (MQ=255)                                                                                       tGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCa                                                                                                              >  1:768141/1‑70 (MQ=255)                                                                                       tGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACg                                                                                                >  1:1631288/1‑84 (MQ=255)                                                                                       tGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACg                                                                                                <  2:1631288/84‑1 (MQ=255)                                                                                       tGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCg                                                                      <  2:171926/110‑1 (MQ=255)                                                                                       tGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCg                                                                      >  1:171926/1‑110 (MQ=255)                                                                                        ggTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGc                                        <  1:354204/139‑1 (MQ=255)                                                                                         gTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCa                                       >  1:1708563/1‑139 (MQ=255)                                                                                             tGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTa                                   <  2:31968/139‑1 (MQ=255)                                                                                                   gctgGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAg                                                                                                     >  1:1006376/1‑67 (MQ=255)                                                                                                   gctgGACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAg                                                                                                     <  2:1006376/67‑1 (MQ=255)                                                                                                       gACGGTGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGaa                         <  2:651187/139‑1 (MQ=255)                                                                                                            tGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGa                                                                                                  <  2:1817124/61‑1 (MQ=255)                                                                                                            tGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGa                                                                                                  >  1:1817124/1‑61 (MQ=255)                                                                                                            tGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCg                                                                      <  2:10640/89‑1 (MQ=255)                                                                                                            tGGCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCg                                                                      >  1:10640/1‑89 (MQ=255)                                                                                                              gCAAAAGCGGAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGc                  <  2:1752790/139‑1 (MQ=255)                                                                                                                      ggAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTAc                                                                   >  2:1302818/1‑82 (MQ=255)                                                                                                                      ggAACTCTGCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTAc                                                                   <  1:1302818/82‑1 (MQ=255)                                                                                                                              gCGGTAACAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGAttt  <  1:859805/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                     GCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTT  >  NC_000913/1148560‑1148823

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