Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,426,491	A→T	6.0%	D147V (GAC→GTC) ‡	PP_3042 →	terminase large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,426,491	0	A	T	6.0%	201.8 / 2.7	66	D147V (GAC→GTC) ‡	PP_3042	terminase large subunit
Reads supporting (aligned to +/- strand):  ref base A (30/32);  new base T (2/2);  total (32/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.69e-01

CTCCGGCCAACATCTTGACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGACGTTAACCAAGAAGGTGATC  >  NC_002947/3426347‑3426628                                                                                                                                                  |                                                                                                                                          cTCCGGCCAACATCTTGACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGacac                                                                                                                                        >  8:91947/1‑148 (MQ=255)         aaCATCTTGACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAAtt                                                                                                                               >  5:52184/1‑149 (MQ=255)          aCATCTTGACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTc                                                                                                                              <  7:310370/149‑1 (MQ=255)               ttGACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGTGTCCAAGGTAttacg                                                                                                                            >  5:105591/1‑143 (MQ=255)               ttGACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGTGTCCAAGGTAttacg                                                                                                                            <  6:105591/146‑4 (MQ=255)                 gACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGAtggtg                                                                                                                       >  2:344166/1‑149 (MQ=255)                 gACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGAtg                                                                                                                          <  8:129487/146‑1 (MQ=255)                 gACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGAtg                                                                                                                          <  6:52184/146‑1 (MQ=255)                 gACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGAtg                                                                                                                          <  6:80221/146‑1 (MQ=255)                 gACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGAt                                                                                                                           <  7:201109/145‑1 (MQ=255)                    gCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAgg                                                                                                                                   <  4:172880/134‑1 (MQ=255)                    gCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAgg                                                                                                                                   >  3:172880/1‑134 (MQ=255)                                cGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCa                                                                                                        <  7:300496/149‑1 (MQ=255)                                       ttGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGc                                                                                                 >  2:200024/1‑149 (MQ=255)                                         ggCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCg                                                                                               >  7:125317/1‑149 (MQ=255)                                               gCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGc                                                                                                 <  3:51139/141‑1 (MQ=255)                                               gCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCGTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGc                                                                                                 >  4:51139/1‑141 (MQ=255)                                                cGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGc                                                                                        <  1:344166/149‑1 (MQ=255)                                                   gtgtCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGc                                                                                     >  5:292243/1‑149 (MQ=255)                                                   gtgtCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGCTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGAAACCAAGGAATTCGATGCAGGTGCATTCTATGCCACTACAGCCGGTTCGGCAGc                                                                                     >  4:240545/1‑149 (MQ=255)                                                        ttCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAAcc                                                                                >  3:156993/1‑149 (MQ=255)                                                        ttCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAAcc                                                                                >  1:249589/1‑149 (MQ=255)                                                         tCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGc                                                                                                          >  3:212992/1‑122 (MQ=255)                                                         tCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGc                                                                                                          <  4:212992/122‑1 (MQ=255)                                                         tCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAAcc                                                                                >  3:48854/1‑148 (MQ=255)                                                          cACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAAc                                                                                 <  7:91947/146‑1 (MQ=255)                                                          cACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCt                                                                               <  4:48854/148‑1 (MQ=255)                                                          cACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATCGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCAAACCt                                                                               >  5:310531/1‑148 (MQ=255)                                                                   gCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCt                                                                           >  3:275139/1‑143 (MQ=255)                                                                   gCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCt                                                                           <  4:275139/143‑1 (MQ=255)                                                                     aaGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTc                                                                   >  2:77106/1‑149 (MQ=255)                                                                     aaGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTc                                                                   >  5:213914/1‑149 (MQ=255)                                                                          cATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTc                                                              >  1:52748/1‑149 (MQ=255)                                                                          cATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTc                                                              <  1:193110/149‑1 (MQ=255)                                                                          cATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTc                                                              >  7:94719/1‑149 (MQ=255)                                                                           aTTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCg                                                             <  1:172340/149‑1 (MQ=255)                                                                              gATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGt                                                           <  8:125317/148‑1 (MQ=255)                                                                                     tACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTAt                                                   >  8:178149/1‑149 (MQ=255)                                                                                        cAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTAGATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATgg                                                 >  1:81599/1‑149 (MQ=255)                                                                                        cAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGc                                                >  4:309307/1‑149 (MQ=255)                                                                                             cTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGa                                           <  1:200024/149‑1 (MQ=255)                                                                                              tGAAAGCGCGTGTTGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAg                                          <  3:240545/149‑1 (MQ=255)                                                                                              tGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAg                                          >  8:70068/1‑149 (MQ=255)                                                                                                        tgtgGCGTCGGTTCGTTCCAGAGTCTCGCGTA‑ATAGCTTGGAGACCAAGGAATTCGATGGTGGTGCTTTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCt                                <  6:310531/149‑1 (MQ=255)                                                                                                          tgGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTgg                              <  6:292243/149‑1 (MQ=255)                                                                                                          tgGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTgg                              <  3:309307/149‑1 (MQ=255)                                                                                                             cGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAg                           <  1:77106/149‑1 (MQ=255)                                                                                                                 ggTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTg                        <  2:52748/148‑1 (MQ=255)                                                                                                                  gttcgttcCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTg                               <  8:113818/140‑1 (MQ=255)                                                                                                                  gttcgttcCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTg                               >  7:113818/1‑140 (MQ=255)                                                                                                                  gttcgttcCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGa                      <  4:156993/149‑1 (MQ=255)                                                                                                                   ttcgttcCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTg                        <  6:213914/146‑1 (MQ=255)                                                                                                                   ttcgttcCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGAc                     >  3:274237/1‑149 (MQ=255)                                                                                                                      gttcCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGACGt                   >  1:195127/1‑148 (MQ=255)                                                                                                                      gtgccAGATATTCCCGTA‑ATAACGTCGACACCAATGAATTCGATGGGTTTTAATTGTTTGCCACTCCAGCCGGTTCGGCAGCCAACCTATCTGAGCTTGAGGATCGTTATGTTTATGGCGATGATGTTGATCTATGGGAGGTGGACGt                   <  8:286014/145‑1 (MQ=255)                                                                                                                          cAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGACg                    >  5:284470/1‑143 (MQ=255)                                                                                                                          cAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGACg                    <  6:284470/143‑1 (MQ=255)                                                                                                                                 tCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGATGGGAGGTGGACGTTAACCaagaagt       >  5:9896/1‑148 (MQ=255)                                                                                                                                   ccGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGa                            >  8:195197/1‑126 (MQ=255)                                                                                                                                   ccGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGa                            <  7:195197/126‑1 (MQ=255)                                                                                                                                   atcga‑aTTCCTTGGTGTCCAAGGTATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGACGTTa                 <  4:308352/133‑1 (MQ=255)                                                                                                                                   atcga‑aTTCCTTGGTGTCCAAGGTATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGACGTTa                 >  3:308352/5‑137 (MQ=255)                                                                                                                                   acGTA‑ATACCTTGGACAGCAATGAATTCGATGGTGGTGCATTGTTTGCCTCTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGACGTTAACCAAGAAgg       <  8:94719/146‑1 (MQ=255)                                                                                                                                    cGTA‑ATAGCTTGGACTCCATGGAATTCGATGGTGGTGAATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGTCTATGAGGTTGATCGCTGGGAGGTTGACGTAAACCAAGAAGGTGa    <  2:249589/149‑1 (MQ=255)                                                                                                                                      tA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGACGTTAACCAAGAAGGTGATc  >  5:235289/1‑149 (MQ=255)                                                                                                                                              ttGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGACGTTaa                <  2:163423/128‑1 (MQ=255)                                                                                                                                              ttGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGACGTTaa                >  1:163423/1‑128 (MQ=255)                                                                                                                                                  |                                                                                                                                          CTCCGGCCAACATCTTGACGCTACTGCCCACCGGGGGCTTGGCGAAGCGGGTGTCTTCACGGATTGGCAAGACCATTGATTCAGTACCAGAACTGAAAGCGCGTGTGGCGTCGGTTCGTTCCAGAGACTCCCGTA‑ATACCTTGGACACCAAGGAATTCGATGGTGGTGCATTGTTTGCCACTACAGCCGGTTCGGCAGCCAACCTATCTGAGCTGTCGGCTCGGTATGTGTATGGCGATGAGGTTGATCGCTGGGAGGTGGACGTTAACCAAGAAGGTGATC  >  NC_002947/3426347‑3426628

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