Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,597,128	A→T	6.5%	V894D (GTC→GAC)	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,597,128	0	A	T	6.5%	235.4 / 5.0	77	V894D (GTC→GAC)	PP_4071	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (36/36);  new base T (2/3);  total (38/39)
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.51e-01

TGATTGATTGGCTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGG  >  NC_002947/4596981‑4597273                                                                                                                                                     |                                                                                                                                                  tgattgatTGGCTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGAc                                                                                                                                                  >  3:127566/1‑149 (MQ=255)  gattgattGGCTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACg                                                                                                                                                 >  8:350375/1‑149 (MQ=255)      gattgGCTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGc                                                                                                                                               <  6:393646/147‑1 (MQ=255)      gattgGCTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGc                                                                                                                                               >  5:393646/1‑147 (MQ=255)      gattgGCTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCg                                                                                                                                              <  4:378918/148‑1 (MQ=255)       attgGCTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCACGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTAACGCCTTGGCATAGTCGTGGACGGCGa                                                                                                                                             >  1:171582/1‑148 (MQ=255)          ggCTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTAt                                                                                                                                         >  5:109607/1‑149 (MQ=255)           gCTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATg                                                                                                                                        <  8:325063/149‑1 (MQ=255)            cTCCGGCAAGCCAGTGGCAACGCGTTGCTAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATgg                                                                                                                                       >  1:217957/1‑149 (MQ=255)            cTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATgg                                                                                                                                       <  4:61725/149‑1 (MQ=255)            cTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATgg                                                                                                                                       <  2:171582/149‑1 (MQ=255)               cGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTTCCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCt                                                                                                                                     <  5:307178/148‑1 (MQ=255)               cGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGCCGTCCACGACTATgcca                                                                                                                                     >  8:101818/1‑145 (MQ=255)               cGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGCCGTCCACGACTATgcca                                                                                                                                     <  7:101818/148‑4 (MQ=255)                 gCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCg                                                                                                                                  >  7:271149/1‑149 (MQ=255)                  cAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGCCATAGTCGCCGTCCACGACTATgccaaggc                                                                                                                                 <  7:132600/149‑8 (MQ=255)                      ccAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGt                                                                                                                             <  6:109607/149‑1 (MQ=255)                            gCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCAt                                                                                                                        <  6:58700/148‑1 (MQ=255)                            gCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATc                                                                                                                       <  8:271149/149‑1 (MQ=255)                            gCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATc                                                                                                                       >  7:115903/1‑149 (MQ=255)                              aaCGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAATAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTGTCATCg                                                                                                                      >  7:360422/1‑148 (MQ=255)                               aCGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCgta                                                                                                                    <  2:242424/149‑1 (MQ=255)                                 gcgTTGCAAGGTGCCGATGGGGTCAGTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGGTGCGTGGGTTTCATc                                                                                                                       >  7:240917/1‑145 (MQ=255)                                 gcgTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATc                                                                                                                       <  8:240917/144‑1 (MQ=255)                                 gcgTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCgtagt                                                                                                                  <  2:217957/149‑1 (MQ=255)                                     tGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCt                                                                                                              >  4:106912/1‑149 (MQ=255)                                       cAAGGTGCCGATGGGGTG‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCGTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTc                                                                                                             <  6:241084/148‑1 (MQ=255)                                              ccGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGTCGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTc                                                                                                     >  2:345791/1‑149 (MQ=255)                                                   ggggTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTg                                                                                                >  2:133050/1‑149 (MQ=255)                                                     ggTC‑GTGCCCGGTCATTGAAAAGCGGTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCa                                                                                              <  1:283587/149‑1 (MQ=255)                                                            gCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTg                                                                                         >  6:126960/1‑148 (MQ=255)                                                                ggTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCgtagta                                                                                                                 <  1:48871/120‑1 (MQ=255)                                                                ggTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCgtagta                                                                                                                 >  2:48871/1‑120 (MQ=255)                                                                    aTTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGgc                                                                                >  7:218964/1‑149 (MQ=255)                                                                    aTTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGgc                                                                                <  1:44143/149‑1 (MQ=255)                                                                        aaaaGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCt                                                                             >  5:1907/1‑148 (MQ=255)                                                                             cGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGcc                                                                       >  6:304004/1‑149 (MQ=255)                                                                              gTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCt                                                                      >  3:374167/1‑149 (MQ=255)                                                                              aaggctGCGCTGCACGCAGTGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCt                                                                      >  7:354947/6‑149 (MQ=38)                                                                                     cACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCa                                                                   >  4:179598/1‑145 (MQ=255)                                                                                     cACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCa                                                                   <  3:179598/145‑1 (MQ=255)                                                                                      aCTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGc                                                              >  3:349324/1‑149 (MQ=255)                                                                                             gcagcgcagcCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTg                                                        >  2:148238/1‑148 (MQ=255)                                                                                             gcagcgcagcCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTg                                                        >  5:200061/1‑148 (MQ=255)                                                                                               agcgcagcCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAAt                                                     <  7:350375/149‑1 (MQ=255)                                                                                               agcgcagcCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAAt                                                     <  1:189206/149‑1 (MQ=255)                                                                                                      ccTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGTCGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAgc                                              <  1:345791/149‑1 (MQ=255)                                                                                                      ccTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAgc                                              >  1:45104/1‑149 (MQ=255)                                                                                                            cTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCtgt                                        >  7:212775/1‑149 (MQ=255)                                                                                                            cTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCtgt                                        <  8:354947/149‑1 (MQ=255)                                                                                                               ggTCTTGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTGATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCtgttg                                      <  4:127566/148‑1 (MQ=255)                                                                                                                gTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTgttgtt                                    <  6:1907/149‑1 (MQ=255)                                                                                                                 tCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGttgttg                                   <  6:200061/149‑1 (MQ=255)                                                                                                                     ggCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCt                                 <  5:126960/147‑1 (MQ=255)                                                                                                                     ggCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCa                               <  3:88723/149‑1 (MQ=255)                                                                                                                     ggCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTACGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCa                               <  8:360422/149‑1 (MQ=255)                                                                                                                        tACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGc                            >  2:285955/1‑149 (MQ=255)                                                                                                                          cccTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCtt                          >  4:118311/1‑149 (MQ=255)                                                                                                                          cccTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACACCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCtt                          >  3:52661/1‑149 (MQ=255)                                                                                                                          cccTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACACCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCtt                          >  3:47631/1‑149 (MQ=255)                                                                                                                           ccTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTg                         >  2:49149/1‑149 (MQ=255)                                                                                                                               tACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCt                       <  3:106912/147‑1 (MQ=255)                                                                                                                               tACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGt                     <  2:45104/149‑1 (MQ=255)                                                                                                                               tACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGt                     <  8:212775/149‑1 (MQ=255)                                                                                                                                 ccGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGtt                    <  7:247346/148‑1 (MQ=255)                                                                                                                                 ccGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTg                   <  4:349324/149‑1 (MQ=255)                                                                                                                                   gCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGAt                 <  3:118311/149‑1 (MQ=255)                                                                                                                                    ccTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATc                >  3:175631/1‑149 (MQ=255)                                                                                                                                     cTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGGTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATc                >  4:297841/1‑148 (MQ=255)                                                                                                                                         gCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGaaa            <  1:148238/148‑1 (MQ=255)                                                                                                                                           aTAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACgg         >  7:119680/1‑149 (MQ=255)                                                                                                                                           aTAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACgg         >  6:204171/1‑149 (MQ=255)                                                                                                                                              gTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCgt                                                                                                                     >  3:86470/1‑38 (MQ=255)                                                                                                                                              gTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCgt                                                                                                                     <  4:86470/38‑1 (MQ=255)                                                                                                                                               tCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCg      >  6:364552/1‑148 (MQ=255)                                                                                                                                                cGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGaa    <  8:115903/149‑1 (MQ=255)                                                                                                                                                  tGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAgg  <  5:304004/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TGATTGATTGGCTCCGGCAAGCCAGTGGCAACGCGTTGCAAGGTGCCGATGGGGTC‑GTGCCCGGTCATTGAAAAGCGTTGATGCACTGCGTGCAGCGCAGCCTTGCCTCGGTCTGGGCTACCCTGTACCGCCTTGGCATAGTCGTGGACGGCGACTATGGCTGCGTGGGTTTCATCGTAGTAGCTCGGTTTTTCGCCTGCAGCGTGCAACATGGCGCTCAAGCCTGCAAACGCCCGCTGAATGACCAGCGCCTGTTGTTGCTCAGGCTTGCTGTTGATCGAAACGGGCGAAGG  >  NC_002947/4596981‑4597273

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