Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,581,441	G→A	7.3%	S165N (AGC→AAC)	PP_3160 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,581,441	0	G	A	7.3%	156.4 / 2.8	55	S165N (AGC→AAC)	PP_3160	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base A (2/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GTGCTGATGCGGTTCAATGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCCCAGCCAGGATAAAA  >  NC_002947/3581293‑3581589                                                                                                                                                     |                                                                                                                                                      gTGCTGATGCGGTTCAATGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCagcagc                                                                                                                                                      >  1:672731/1‑150 (MQ=255) gTGCTGATGCGGTTCAATGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCagcagc                                                                                                                                                      <  1:179127/150‑1 (MQ=255)             ttCAATGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTTACTATCCGCCCAGCAGCTTCAT‑TGCgggc                                                                                                                                         <  1:807339‑M2/150‑5 (MQ=255)              tCAATGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGTTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCgggcg                                                                                                                                        >  1:825302‑M2/1‑145 (MQ=255)                 aTGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGtt                                                                                                                                     <  2:1817783/150‑1 (MQ=255)                  tGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCAAGTTAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGttt                                                                                                                                    <  1:233994/150‑1 (MQ=255)                       cAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGtttattta                                                                                                                               <  1:648345/150‑1 (MQ=255)                        aGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTAc                                                                                                                              >  2:798313/1‑150 (MQ=255)                           cAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCaa                                                                                                                           <  1:658567/150‑1 (MQ=255)                              cTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCg                                                                                                                        >  1:986778/1‑150 (MQ=255)                                     ccGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCcagcag                                                                                                                                                       >  1:1909753/1‑113 (MQ=255)                                     ccGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCcagcag                                                                                                                                                       <  2:1909753/113‑1 (MQ=255)                                     ccGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGa                                                                                                                 <  2:1524825/150‑1 (MQ=255)                                     ccGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGa                                                                                                                 <  1:452757/150‑1 (MQ=255)                                     ccGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCC‑GCAA‑TGAAGCTGCtgggcggatagtcaacgtgaccagttgga                                                                                                                <  2:1536825‑M2/150‑30 (MQ=255)                                             tGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTc                                                                                                         <  2:672731/150‑1 (MQ=255)                                                tGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCC‑GCAATGAAGC‑TGCtgggcggatagtcaacgtgaccagttggaacgtgaccact                                                                                                     >  1:1536825‑M2/1‑110 (MQ=255)                                                    ggggCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGca                                                                                                  <  1:518296/150‑1 (MQ=255)                                                         cAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCata                                                                                             >  1:1310357/1‑150 (MQ=255)                                                         cAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCC‑GCAA‑TGAAGCTGCtgggcggatagtcaacgtgaccagttggaacgt                                                                                                            >  2:1449680‑M2/1‑101 (MQ=255)                                                         cAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCC‑GCAA‑TGAAGCTGCtgggcggatagtcaacgtgaccagttggaacgt                                                                                                            <  1:1449680‑M2/134‑34 (MQ=255)                                                                tATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCAGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTgg                                                                                      >  2:495294/1‑150 (MQ=255)                                                                 aTTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTggg                                                                                     >  2:244136/1‑150 (MQ=255)                                                                    gTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGtt                                                                                                                                     <  1:446433/99‑1 (MQ=255)                                                                    gTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGtt                                                                                                                                     >  2:446433/1‑99 (MQ=255)                                                                    gTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTAACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATAATATCAGTggggcc                                                                                  >  2:619580/1‑148 (MQ=255)                                                                          gCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAg                                                                            <  2:752980/150‑1 (MQ=255)                                                                           cAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCAGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGcat                                                                                                 <  2:96853/128‑1 (MQ=255)                                                                           cAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCAGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGcat                                                                                                 >  1:96853/1‑128 (MQ=255)                                                                                   aTCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGcatca                                                                                               >  1:1307133/1‑122 (MQ=255)                                                                                   aTCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGcatca                                                                                               <  2:1307133/122‑1 (MQ=255)                                                                                    tCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCtt                                                                  >  2:152964/1‑150 (MQ=255)                                                                                        gTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCCTATCAGTGGGGTCAAAAAGCTCCGGCCTTGagc                                                              >  1:81491/1‑148 (MQ=255)                                                                                         tGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTggg                                                                                     >  2:1549821/1‑126 (MQ=255)                                                                                         tGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTggg                                                                                     <  1:1549821/126‑1 (MQ=255)                                                                                         tGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGaaaa                                                             <  2:986778/150‑1 (MQ=255)                                                                                           aaCAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATggg                                                                                                                   >  2:73996/1‑94 (MQ=255)                                                                                           aaCAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATggg                                                                                                                   <  1:73996/94‑1 (MQ=255)                                                                                            aCAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATTGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAAccc                                                          >  2:1928877/1‑150 (MQ=255)                                                                                                 gTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGgt                                                     <  1:244136/150‑1 (MQ=255)                                                                                                      tCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCa                                                >  1:1132911/1‑150 (MQ=255)                                                                                                         gAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGaa                                             >  2:1162723/1‑150 (MQ=255)                                                                                                                  cGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGa                                    <  1:1217802/150‑1 (MQ=255)                                                                                                                   gTTCCAACTGGTCACGTTTACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATAAAAAAAAAAAAAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGAt                                   <  2:807339/150‑1 (MQ=255)                                                                                                                   gTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGAt                                                                                                                      >  2:252218/1‑67 (MQ=255)                                                                                                                   gTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGAt                                                                                                                      <  1:252218/67‑1 (MQ=255)                                                                                                                    ttCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATg                                  >  2:583133/1‑150 (MQ=255)                                                                                                                     tCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTggg                                                                                     <  2:1652539/98‑1 (MQ=255)                                                                                                                     tCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTggg                                                                                     >  1:1652539/1‑98 (MQ=255)                                                                                                                        aaCTGGTCACGTTGACTATCCGCCCATCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGCGCAGAAGGCAGTGGATGAAtt                              >  2:1744985/1‑150 (MQ=255)                                                                                                                                   ttGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCtt                                                                  <  2:100984/103‑1 (MQ=255)                                                                                                                                   ttGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCtt                                                                  >  1:100984/1‑103 (MQ=255)                                                                                                                                   ttGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTAccc                   <  2:1132911/150‑1 (MQ=255)                                                                                                                                     gACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTc                 <  1:798313/150‑1 (MQ=255)                                                                                                                                          tCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCccagc            >  2:1430986/1‑150 (MQ=255)                                                                                                                                            cGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCccagcca          >  1:1403039/1‑150 (MQ=255)                                                                                                                                            cGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCccagcca          <  2:1115135/150‑1 (MQ=255)                                                                                                                                             gCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCccagccag         >  1:1654495/1‑150 (MQ=255)                                                                                                                                             gCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCccagccag         >  1:811191/1‑150 (MQ=255)                                                                                                                                                    agcTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCGTGAAAACCCTCGGTGCGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCCCAGCCAGGATaaaa  <  1:1744985/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GTGCTGATGCGGTTCAATGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCCCAGCCAGGATAAAA  >  NC_002947/3581293‑3581589

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑