Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	734,326	C→T	5.8%	A331T (GCG→ACG)	ndh ←	NADH dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	734,326	0	C	T	5.8%	211.9 / 2.2	69	A331T (GCG→ACG)	ndh	NADH dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base T (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

CATCAGGTTGCCCACCGCCGAGAAGCGCGACAGCGACACCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTGGTTGATGCGGTTGGTTTCCAGGCCATCGATG  >  NC_002947/734189‑734474                                                                                                                                           |                                                                                                                                                     cATCAGGTTGCCCACCGCCGAGAAGCGCGACAGCGACACCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTATCCAG‑GCGCGCCTTCAGGCTTGGCGCCCGCATCGAAGCCTGCTGGTGCGCCGCCTGTGcgcgcg                                                                                                                                           >  6:94814/1‑149 (MQ=255)  aTCAGGTTGCCCACCGCCGAGAAGCGCGACAGCGACACCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCgcgg                                                                                                                                          >  1:115720/1‑149 (MQ=255)       gTTGCCCACCGCCGAGAAGCGCGACAGCGACACCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCa                                                                                                                                     >  5:22471/1‑149 (MQ=255)           cccACCGCCGAGAAGCGCGACAGCGACACCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGt                                                                                                                                  >  5:342149/1‑148 (MQ=255)           cccACCGCCGAGAAGCGCGACAGCGACACCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGt                                                                                                                                  >  3:357210/1‑148 (MQ=255)                cgccgAGAAGCGCGACAGCGACACCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGAt                                                                                                                            >  8:280606/1‑149 (MQ=255)                             gacagcgacaCCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATc                                                                                                                           <  7:326318/137‑1 (MQ=255)                             gacagcgacaCCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATc                                                                                                                           >  8:326318/1‑137 (MQ=255)                              acagcgacaCCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGc                                                                                                              <  3:224128/149‑1 (MQ=255)                               cagcgacaCCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCTGCGGCACGTTACGATCGCTGCCGGGTTGCg                                                                                                             <  6:350699/149‑1 (MQ=255)                               cagcgacaCCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCg                                                                                                             >  7:15471/1‑149 (MQ=255)                               cagcgacaCCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCg                                                                                                             <  4:357210/149‑1 (MQ=255)                                agcgacaCCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGTCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCgg                                                                                                            >  7:261808/1‑149 (MQ=255)                                agcgacaCAAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGTCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCgg                                                                                                            >  7:256716/1‑149 (MQ=255)                                     cacCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAgg                                                                                                       >  4:313951/1‑149 (MQ=255)                                        cAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAGAGCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAggcg                                                                                                     >  1:97223/1‑149 (MQ=255)                                          gCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCCGGTTTAGGGCAGgcggcg                                                                                                  >  6:269431/1‑149 (MQ=255)                                                cATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCa                                                                                            >  8:63449/1‑149 (MQ=255)                                                 aTAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCAc                                                                                           >  7:309269/1‑149 (MQ=255)                                                 aTAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCAc                                                                                           >  7:267598/1‑149 (MQ=255)                                                 aTAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCAc                                                                                           <  7:280993/149‑1 (MQ=255)                                                  tAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCAcc                                                                                          <  2:97223/149‑1 (MQ=255)                                                     tCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGaa                                                                                       <  8:162206/149‑1 (MQ=255)                                                      ccTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAg                                                                                      >  3:271877/1‑149 (MQ=255)                                                       cTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAgg                                                                                     >  7:287820/1‑149 (MQ=255)                                                       cTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAg                                                                                      <  8:15471/148‑1 (MQ=255)                                                           tACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGaa                                                                                 <  2:115720/149‑1 (MQ=255)                                                           tACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGa                                                                                  <  5:46486/148‑1 (MQ=255)                                                           tACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCAGCCTGTGCTAGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGaa                                                                                 <  5:269431/149‑1 (MQ=255)                                                             cTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAg                                                                                <  3:353009/148‑1 (MQ=255)                                                              tCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGAt                                                                              >  6:258776/1‑149 (MQ=255)                                                               cGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGAtg                                                                             >  8:9443/1‑149 (MQ=255)                                                                gTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGTCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGAtg                                                                             <  8:261808/148‑1 (MQ=255)                                                                gTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGTCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGAtg                                                                             <  8:256716/148‑1 (MQ=255)                                                                gTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGAtgt                                                                            <  8:309269/149‑1 (MQ=255)                                                                gTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGAtgt                                                                            >  6:333767/1‑149 (MQ=255)                                                                gTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGAtg                                                                             <  6:342149/148‑1 (MQ=255)                                                                 tACGTCGGCAGGGCCTTGTGCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGAtgtt                                                                           <  7:65941/149‑1 (MQ=255)                                                                     tCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTgtcg                                                                       >  2:253173/1‑149 (MQ=255)                                                                       ggCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTgtcgt                                                                      <  7:9443/148‑1 (MQ=255)                                                                         cAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTcgc                                                                   >  6:75308/1‑149 (MQ=255)                                                                           gggCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCgg                                                                 <  6:22471/149‑1 (MQ=255)                                                                            ggCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCggg                                                                                                           <  3:190358/106‑1 (MQ=255)                                                                            ggCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCggg                                                                                                           >  4:190358/1‑106 (MQ=255)                                                                               cTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGggtgg                                                             <  6:9171/149‑1 (MQ=255)                                                                                   ttCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTg                                                         <  8:267598/149‑1 (MQ=255)                                                                                      tcCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAg                                                      >  8:42841/1‑149 (MQ=255)                                                                                      tcCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAg                                                      >  3:127607/1‑149 (MQ=255)                                                                                        cAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAggcggc                                                                                                   <  2:51364/102‑1 (MQ=255)                                                                                        cAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAggcggc                                                                                                   >  1:51364/1‑102 (MQ=255)                                                                                        cAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCa                                                       <  4:271877/146‑1 (MQ=255)                                                                                        cAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCa                                                       <  7:2334/146‑1 (MQ=255)                                                                                        cAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAgg                                                     <  2:159835/148‑1 (MQ=255)                                                                                         aG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGt                                                   <  6:170692/149‑1 (MQ=255)                                                                                             cgcgcCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAgg                                                <  7:280606/149‑1 (MQ=255)                                                                                              gcgcCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAggg                                               <  4:127607/149‑1 (MQ=255)                                                                                                 ccTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGgcgc                                            <  8:287820/149‑1 (MQ=255)                                                                                                  cTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCAGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCa                                           >  5:363170/1‑149 (MQ=255)                                                                                                   ttCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGaa                                                                                 <  3:79169/110‑1 (MQ=255)                                                                                                   ttCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGaa                                                                                 >  4:79169/1‑110 (MQ=255)                                                                                                       ggCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGAc                                       >  6:19747/1‑148 (MQ=255)                                                                                                           ttGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCt                                  >  6:77827/1‑149 (MQ=255)                                                                                                              cgcCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTGGt                               <  6:112231/149‑1 (MQ=255)                                                                                                              cgcCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTGGt                               >  5:112231/1‑149 (MQ=255)                                                                                                                    cATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTgg                                <  2:25673/142‑1 (MQ=255)                                                                                                                    cATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTgg                                >  1:25673/1‑142 (MQ=255)                                                                                                                           gCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGTTCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTGGTTGATGCggttggt                  <  5:333767/149‑1 (MQ=255)                                                                                                                                  ggTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTGGTTGATGCGGTTGGTTTCCAgg           >  5:330619/1‑149 (MQ=255)                                                                                                                                   gTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTGGTTGATGCGGTTGGTTTCCAGGc          >  6:296130/1‑149 (MQ=255)                                                                                                                                      cgccgcCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTGGTTGATGCGGTTGGTTTCCAGGCCAt       <  4:378170/149‑1 (MQ=255)                                                                                                                                        ccgccTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTGGTTGATGCGGTTGGTTTCCAGGCCATc      >  3:369151/1‑148 (MQ=255)                                                                                                                                           ccTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTGGTTGATGCGGTTGGTTTCCAGGCCATCGATg  <  5:77827/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                     CATCAGGTTGCCCACCGCCGAGAAGCGCGACAGCGACACCAGCGAGCCATAGTCCTTGTACTCGTACGTCGGCAGGGCCTTGTTCTCCAG‑GCGCGCCTTCAGGCTTTGCGCCAGCATCGAAGCCTGCTGGTGCGCCGCCTGTGCGCGCGGCGGCACGTTACGATCGCTGCCGGGTTGCGGGCAGGCGGCGCAGTCACCGAAGGCGAAGATGTTGTCGTCGCGGGTGGTTTGCAGGGTAGGGCGCACGACCAGCTGGTTGATGCGGTTGGTTTCCAGGCCATCGATG  >  NC_002947/734189‑734474

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: ‑