Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,490,902	T→A	9.1%	S178C (AGC→TGC)	clpV ←	protein ClpV1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,490,902	0	T	A	9.1%	264.7 / 14.6	90	S178C (AGC→TGC)	clpV	protein ClpV1
Reads supporting (aligned to +/- strand):  ref base T (40/41);  new base A (3/5);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

GTCGACCATCTGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGA  >  NC_002947/3490753‑3491044                                                                                                                                                       |                                                                                                                                               gTCGACCATCTGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTGTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCgctgggct                                                                                                                                                >  2:1741155/1‑150 (MQ=255)  tCGACCATCTGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCgctgggctg                                                                                                                                               <  2:8335/150‑1 (MQ=255)     aCCATCTGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCTATCCCGCCGTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGcc                                                                                                                                            >  1:1758810/1‑150 (MQ=255)     aCCATCTGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGcc                                                                                                                                            >  1:811415/1‑150 (MQ=255)     aCCATCTGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGcc                                                                                                                                            >  1:1919567/1‑150 (MQ=255)        aTCTGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGc                                                                                                                                         <  1:1278073/150‑1 (MQ=255)           tGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTCCCGCGGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGCCAGCCCAGCggc                                                                                                                                      <  1:559407/150‑2 (MQ=255)            gCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGc                                                                                                                                     >  2:236968/1‑150 (MQ=255)              cGTACTTCGAACTCACGGCCCAGCACCTGTTTTATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                   <  1:1020270/150‑1 (MQ=255)              cGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                   <  2:1328989/150‑1 (MQ=255)              cGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGAGGCGTCCGCGAGTTTACTCATGGACTGGTCCAGCGCCGGGGTCTTGCCGCTGGACTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                   <  1:285358/150‑1 (MQ=255)                tACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGATCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGgc                                                                                                                                 <  1:844737/150‑1 (MQ=255)                tACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGgc                                                                                                                                 <  1:1576971/150‑1 (MQ=255)                   ttCGAACTCACGGCCCAGCACCGGGTCTATGCGGCCTTCGCGGGCGCTTTGGGTAGGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGGG‑CCGGTGGGCTGGCCAGCGGACTGgcctg                                                                                                                              >  2:2364486/1‑150 (MQ=255)                   ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGgcctg                                                                                                                              >  2:1586983/1‑150 (MQ=255)                   ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGgcctg                                                                                                                              >  1:1110425/1‑150 (MQ=255)                   ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCAGGCCAGCCCGCTGgccag                                                                                                                              <  2:167267/150‑3 (MQ=255)                   ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCAGGCCAGCCCGCTGgccag                                                                                                                              >  1:167267/1‑148 (MQ=255)                   ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCAGGCCAGCCCGCTGgccag                                                                                                                              <  2:1380948/150‑3 (MQ=255)                   ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTACTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGgcctg                                                                                                                              >  1:1512611/1‑150 (MQ=255)                   ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCGCGACGTGCTCACTGGTG‑CCGCTGGGCTGGCCCGCGGGCTGgcctc                                                                                                                              >  2:1708558/1‑149 (MQ=255)                      gAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCCGCTGGCCAGCCCAGCggcgc                                                                                                                                    <  1:1914811/142‑4 (MQ=255)                      gAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCCGCTGGCCAGCCCAGCggcgc                                                                                                                                    >  2:1914811/1‑139 (MQ=255)                       aaCTCACCGCCCAGCACCGGGTCGATGCGTCCTTCGCGGGCGCTTTGGGTAAGGTTGACCTTGTACTGGTCCAGCTCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGcctgcctg                                                                                                                          <  1:1968945/150‑1 (MQ=255)                          tCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGgcctgcccag                                                                                                                         <  2:1458363/148‑4 (MQ=255)                          tCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGgcctgcccag                                                                                                                         >  1:1458363/1‑145 (MQ=255)                          tCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCAGGCCAGCCCGCTGgccagccc                                                                                                                           >  2:250361/1‑145 (MQ=255)                          tCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCAGGCCAGCCCGCTGgccagccc                                                                                                                           <  1:250361/146‑2 (MQ=255)                            aCGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTc                                                                                                                     <  1:1428585/150‑1 (MQ=255)                              ggCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCgg                                                                                                                   <  1:662898/150‑1 (MQ=255)                                 ccAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcgc                                                                                                                <  1:1586983/150‑1 (MQ=255)                                 ccAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcgc                                                                                                                >  2:1324444/1‑150 (MQ=255)                                  cAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcg                                                                                                                 <  2:1281092/148‑1 (MQ=255)                                  cAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcg                                                                                                                 >  1:1281092/1‑148 (MQ=255)                                  cAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTCTGGGTAAGGCTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGCCGTGCTCCCCGGCG‑CCGCTGCGCCGGCCCGCGGGCTGGCCTGCCCGGTCGCGgcgcc                                                                                                               >  1:1018058/1‑149 (MQ=255)                                   aGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTa                                                                                                              <  2:1152510/150‑1 (MQ=255)                                   aGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTa                                                                                                              <  1:836655/150‑1 (MQ=255)                                   aGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACATGCTCACTCGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTa                                                                                                              >  2:2244913/1‑150 (MQ=255)                                                  cGGCCTTCGCCCGCTCTTTGGCTAAGCTTGACCTTGTACTTGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTCCTCCACGTGCTCTTTCTCG‑CCGCTGGGCTGGCCAGCGGCCTGGCCTGCCTTTACTCGGCGCTACCCTCCACCAGcgcg                                                                                               <  2:815483/150‑1 (MQ=255)                                                      cTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCaa                                                                                           >  1:1460015/1‑150 (MQ=255)                                                        tCGCGGGCGCTTTGGGTCTCGGTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAgg                                                                                         <  1:1328555/150‑1 (MQ=255)                                                         cgcgGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCTCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGt                                                                                        <  2:346625/150‑1 (MQ=255)                                                            gggCGCTTTGGGTACGGTTGACAGGGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCa                                                                                     <  2:1758810/150‑1 (MQ=255)                                                            gggCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCa                                                                                     <  2:232607/150‑1 (MQ=255)                                                            gggCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCa                                                                                     >  1:2313370/1‑150 (MQ=255)                                                                  tttGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCCGGCGCa                                                                               >  2:538640/1‑150 (MQ=255)                                                                  tttGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCa                                                                               >  1:1478079/1‑150 (MQ=255)                                                                  tttGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCa                                                                               <  1:1426240/150‑1 (MQ=255)                                                                    tGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGACGATGGAATTGCGCGAGGTGCGCACTGGCA‑CACGTGTGCTGGCCAGCGGGCTTGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAgg                                                                             <  1:1741155/150‑1 (MQ=255)                                                                        tAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcgc                                                                                                                >  1:1238335/1‑111 (MQ=255)                                                                        tAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcgc                                                                                                                <  2:1238335/111‑1 (MQ=255)                                                                        tAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCAt                                                                         <  2:2025524/150‑1 (MQ=255)                                                                         aaGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATc                                                                        <  1:951967/150‑1 (MQ=255)                                                                            gTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAg                                                                                    <  2:1367953/135‑1 (MQ=255)                                                                            gTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAg                                                                                    >  1:1367953/1‑135 (MQ=255)                                                                               gACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGtt                                                                  >  2:1859590/1‑150 (MQ=255)                                                                                 ccGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGa                                                                <  1:333477/150‑1 (MQ=255)                                                                                  cGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAt                                                               >  1:339286/1‑150 (MQ=255)                                                                                   gtgtACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAtt                                                              >  2:1402228/1‑150 (MQ=255)                                                                                        cTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTc                                                         <  1:1423559/150‑1 (MQ=255)                                                                                        cTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTc                                                         <  1:236968/150‑1 (MQ=255)                                                                                           gTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAg                                                      >  1:1821614/1‑150 (MQ=255)                                                                                                gcgcCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCg                                                 <  2:1302446/150‑1 (MQ=255)                                                                                                 cgcCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCgg                                                >  1:1935511/1‑150 (MQ=255)                                                                                                  gcCGGGGTCTTGCCGCTGGCCTTGCTCGACTTGCTCACTGGCG‑CCGCTGGGCTGGCCAGAGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCTCGGCAAGGTTCAGGCGCAGTTCATCGGCGTTGATTTTTTCCAGCTCCGGt                                               <  2:811415/150‑1 (MQ=255)                                                                                                  gcCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGt                                               >  1:72728/1‑150 (MQ=255)                                                                                                    cGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGc                                             >  1:1695976/1‑150 (MQ=255)                                                                                                    cGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGc                                             >  2:66925/1‑150 (MQ=255)                                                                                                      gggTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACGAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCg                                           >  2:1277774/1‑150 (MQ=255)                                                                                                        gTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAg                                         >  2:585499/1‑150 (MQ=255)                                                                                                        gTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAg                                         <  2:1821614/150‑1 (MQ=255)                                                                                                           ttGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCg                                      >  1:2379933/1‑150 (MQ=255)                                                                                                             gCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGAt                                    >  2:2270789/1‑150 (MQ=255)                                                                                                              ccGCTGCCCGTGGTCGGCGGGCTCCCTGGCG‑CCGCGGGGCTGGCCAGCGGGCTGCCCTGCCTGGACCCGGCGCTACCCCCCCCCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTGGATTTTTTCCAGCGCCGGTGCCGAGGCGATc                                   <  2:1018058/150‑1 (MQ=255)                                                                                                               cGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                   >  1:1016607/1‑53 (MQ=255)                                                                                                               cGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                   <  2:1016607/53‑1 (MQ=255)                                                                                                                       ttGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGa                          >  2:875363/1‑150 (MQ=255)                                                                                                                           tCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCa                      <  1:585499/150‑1 (MQ=255)                                                                                                                             gACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGt                    >  2:186323/1‑150 (MQ=255)                                                                                                                               cGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAATGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCAcgc                            <  2:362542/140‑1 (MQ=255)                                                                                                                               cGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAATGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCAcgc                            >  1:362542/1‑140 (MQ=255)                                                                                                                                gTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCa                 <  2:2313370/150‑1 (MQ=255)                                                                                                                                   cTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCa              >  2:1925996/1‑150 (MQ=255)                                                                                                                                      aCTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAttttt                                                           <  1:1707088/102‑1 (MQ=255)                                                                                                                                      aCTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAttttt                                                           >  2:1707088/1‑102 (MQ=255)                                                                                                                                      aCTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCg           >  1:1578232/1‑150 (MQ=255)                                                                                                                                      aCTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGGTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCg           >  2:607337/1‑150 (MQ=255)                                                                                                                                       cTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACAATACACCAGCGAGGCAAGGTTACGGCCAAGGTCATCCGCGTTGATTTTTTCTAGCTCCGGTTCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGt          <  2:1479427/150‑1 (MQ=255)                                                                                                                                          gcg‑cCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGGGCGTTACGCGACGCAGTTCAGCATCGTCaa       <  1:2244913/150‑1 (MQ=255)                                                                                                                                          gcg‑cCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCaa       >  2:1300961/1‑150 (MQ=255)                                                                                                                                           cg‑ccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAg      <  1:1324444/150‑1 (MQ=255)                                                                                                                                           cg‑ccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAg      <  2:1578232/150‑1 (MQ=255)                                                                                                                                                gctgggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCaga  <  1:1088096/150‑1 (MQ=255)                                                                                                                                                  tgggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGc                                                                     >  2:348229/1‑81 (MQ=255)                                                                                                                                                  tgggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGc                                                                     <  1:348229/81‑1 (MQ=255)                                                                                                                                                    ggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTc         <  1:1945199/139‑1 (MQ=255)                                                                                                                                                    ggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTc         >  2:1945199/1‑139 (MQ=255)                                                                                                                                                       |                                                                                                                                               GTCGACCATCTGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCG‑CCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGA  >  NC_002947/3490753‑3491044

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