Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,976,233	A→C	I8S (ATT→AGT)	motB ←	motility protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,976,233	0	A	C	100.0%	216.4 / NA	63	I8S (ATT→AGT)	motB	motility protein B
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (34/29);  total (34/29)

TTGGGCTGGAGATGGAGATCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATAATCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCAA  >  NC_000913/1976099‑1976367                                                                                                                                       |                                                                                                                                       ttGGGCTGGAGATGGAGATCAGCCACATCACCAG‑AAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGAt                                                                                                                                  >  1:1323503/1‑140 (MQ=255)  tGGGCTGGAGATGGAGATCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGAt                                                                                                                                  <  2:26917/140‑1 (MQ=255)  tGGGCTGGAGATGGAGATCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGAt                                                                                                                                  <  1:472929/140‑1 (MQ=255)  tGGGCTGGAGATGGAGATCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACAACAATACTCGGAt                                                                                                                                  <  1:778933/140‑1 (MQ=255)      cTGGAGATGGAGATCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATgcgc                                                                                                                              >  2:1184214/1‑140 (MQ=255)         gagatggagatCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATgcgc                                                                                                                              <  2:1158814/137‑1 (MQ=255)         gagatggagatCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATgcgc                                                                                                                              >  1:1158814/1‑137 (MQ=255)         gagatggagatCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCtt                                                                                                                            <  1:175149/139‑1 (MQ=255)         gagatggagatCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCtt                                                                                                                            >  2:175149/1‑139 (MQ=255)         gagatggagatCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTg                                                                                                                           >  2:440219/1‑140 (MQ=255)         gagatggagatCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTg                                                                                                                           >  2:910338/1‑140 (MQ=255)         gagatggagatCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTg                                                                                                                           >  1:722325/1‑140 (MQ=255)          agatggagatCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGa                                                                                                                          <  1:1393196/140‑1 (MQ=255)                    cAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGc                                                                                                                <  1:892127/140‑1 (MQ=255)                     aGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCa                                                                                                                   <  1:301465/136‑1 (MQ=255)                     aGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCa                                                                                                                   >  2:301465/1‑136 (MQ=255)                         acaTCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCt                                                                                                           >  2:1325556/1‑140 (MQ=255)                         acaTCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCt                                                                                                           <  2:535622/140‑1 (MQ=255)                         acaTCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCt                                                                                                           <  1:161190/140‑1 (MQ=255)                                 aGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCt                                                                                                           >  1:1029070/1‑132 (MQ=255)                                 aGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCt                                                                                                           <  2:1029070/132‑1 (MQ=255)                                     aaaaGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCt                                                                                               >  1:1063214/1‑140 (MQ=255)                                                       tCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCa                                                                             <  2:1323503/140‑1 (MQ=255)                                                            aaGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACcgca                                                                        <  1:1432073/140‑1 (MQ=255)                                                               tCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACcgcacgc                                                                     >  2:396733/1‑140 (MQ=255)                                                                 ggCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGCTGCGGATTTTTCACCgcacgcac                                                                   >  1:344515/1‑140 (MQ=255)                                                                          aaTCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTc                                                          <  1:396733/140‑1 (MQ=255)                                                                             cTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGAt                                                                                                                                  <  2:980052/65‑1 (MQ=255)                                                                             cTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGAt                                                                                                                                  >  1:980052/1‑65 (MQ=255)                                                                             cTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATg                                                                <  1:1433458/131‑1 (MQ=255)                                                                             cTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATg                                                                >  2:1433458/1‑131 (MQ=255)                                                                               tCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGtt                                                     <  1:910338/140‑1 (MQ=255)                                                                               tCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGtt                                                     <  1:440219/140‑1 (MQ=255)                                                                                 cACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTc                                                          <  1:877091/133‑1 (MQ=255)                                                                                 cACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTc                                                          >  2:877091/1‑133 (MQ=255)                                                                                  aCGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCaa                                                  <  2:1427160/140‑1 (MQ=255)                                                                                   cGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAAt                                                 >  1:1132393/1‑140 (MQ=255)                                                                                      tCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATaaa                                              <  2:679474/140‑1 (MQ=255)                                                                                        cATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAacg                                            >  2:869223/1‑140 (MQ=255)                                                                                               tgcCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTGGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACgacggacg                                     <  2:1359439/140‑1 (MQ=255)                                                                                               tgcCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACgacggacg                                     >  1:1116394/1‑140 (MQ=255)                                                                                                gcCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGt                                    >  2:117172/1‑140 (MQ=255)                                                                                                gcCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGt                                    >  2:777783/1‑140 (MQ=255)                                                                                                gcCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGt                                    >  1:1038227/1‑140 (MQ=255)                                                                                                gcCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGt                                    >  1:1037917/1‑140 (MQ=255)                                                                                                 ccccGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAacgacg                                         >  2:754679/1‑134 (MQ=255)                                                                                                 ccccGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAacgacg                                         <  1:754679/134‑1 (MQ=255)                                                                                                 ccccGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGtt                                   >  2:613542/1‑140 (MQ=255)                                                                                                 ccccGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGtt                                   >  2:989502/1‑140 (MQ=255)                                                                                                 ccccGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGtt                                   >  1:1280705/1‑140 (MQ=255)                                                                                                   ccGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCg                                 <  2:1132393/140‑1 (MQ=255)                                                                                                      tggcttttggctttGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTg                              >  2:1376398/1‑140 (MQ=255)                                                                                                      tggcttttggctttGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTg                              <  2:1063214/140‑1 (MQ=255)                                                                                                      tggcttttggctttGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTg                              >  1:769070/1‑140 (MQ=255)                                                                                                        gcttttggctttGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGa                            <  2:1037917/140‑1 (MQ=255)                                                                                                            ttggctttGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGtt                                   <  2:1106142/129‑1 (MQ=255)                                                                                                            ttggctttGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGtt                                   >  1:1106142/1‑129 (MQ=255)                                                                                                             tggctttGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATaga                       >  2:246555/1‑140 (MQ=255)                                                                                                                   tGCGTCGTTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGttt                 >  1:1151133/1‑140 (MQ=255)                                                                                                                     cgtcgtTTGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTg               <  2:344515/140‑1 (MQ=255)                                                                                                                          tttGACGACAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTGCGAcc          <  2:769070/140‑1 (MQ=255)                                                                                                                                  cAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCaa  >  2:1255855/1‑140 (MQ=255)                                                                                                                                  cAATACTCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCaa  >  1:1186330/1‑140 (MQ=255)                                                                                                                                    atactCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTg               <  2:774135/125‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       TTGGGCTGGAGATGGAGATCAGCCACATCACCAGAAAAAAGGCCATCATCGCAGTCATAAAGTCGGCATAAGCAATCTTCCACGATCCATGTGCTGCCCCGTGGCTTTTGGCTTTGCGTCGTTTGACGACAATAATCGGATGCGCTTGATTCTTCATGCTTCCTCGGTTGTCGTCTGTTGTTGCGGATTTTTCACCGCACGCACATGCTCTTCCAGTTCAATAAACGACGGACGTTCGCTGGAATAGAGCGTTTTGCGACCAAACTCAA  >  NC_000913/1976099‑1976367

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