Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,425,661	G→A	9.7%	intergenic (+63/‑442)	PP_1249 → / → PP_1250	DUF4223 domain‑containing protein/group II intron‑encoding maturase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,425,661	0	G	A	9.7%	198.6 / 9.4	73	intergenic (+63/‑442)	PP_1249/PP_1250	DUF4223 domain‑containing protein/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base G (34/31);  new base A (4/3);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGT  >  NC_002947/1425519‑1425809                                                                                                                                               |                                                                                                                                                      aCCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggg                                                                                                                                                 >  3:328950/1‑149 (MQ=255) aCCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGACCCCCGGCCATGTGCATGACCgggg                                                                                                                                                 >  5:353727/1‑149 (MQ=255)    tGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggggcc                                                                                                                                              >  8:213576/1‑147 (MQ=255)      ttCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggggcctt                                                                                                                                            <  7:213576/149‑5 (MQ=255)        ccACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                           <  4:328950/148‑1 (MQ=255)         cACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgt                                                                                                                                         >  3:38130/1‑149 (MQ=255)         cACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGACCCCCGGCCATGTGCATGACCgggggccttgga                                                                                                                                         <  6:353727/149‑8 (MQ=255)           cccTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCACCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                                       >  7:20735/1‑149 (MQ=255)             cTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgc                                                                                                                                     <  5:246274/149‑1 (MQ=255)              ttCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgcg                                                                                                                                    >  8:32212/1‑149 (MQ=255)              ttCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgcg                                                                                                                                    >  3:6656/1‑149 (MQ=255)                  gTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAg                                                                                                                                >  3:73621/1‑149 (MQ=255)                            ccAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGcgcg                                                                                                                      >  5:255580/1‑149 (MQ=255)                             cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGcgcg                                                                                                                      <  4:6656/148‑1 (MQ=255)                             cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGt                                                                                                                     >  4:34787/1‑149 (MQ=255)                                gATCATCGGTGGCTGCGGCACGATAGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGACGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTgc                                                                                                                  >  6:189936/1‑149 (MQ=255)                                  tcatcGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTgc‑gc                                                                                                               <  6:11937/149‑1 (MQ=255)                                       gttGGCTGGGGCCCGATCGATAAACTGCCTAAGCAGCAGTAATCTTGGCTGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGTCCCCCGGTCATGCACATGTCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTTC‑GCGTAt                                                                                                           <  6:195439/146‑1 (MQ=39)                                         tGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATgc                                                                                                         >  2:313111/1‑148 (MQ=255)                                             tGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGgcgc                                                                                                                       >  7:304678/1‑131 (MQ=255)                                             tGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGgcgc                                                                                                                       <  8:304678/131‑1 (MQ=255)                                              gCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgcgc                                                                                                                                   >  1:251386/1‑118 (MQ=255)                                              gCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgcgc                                                                                                                                   <  2:251386/118‑1 (MQ=255)                                               cGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAAcc                                                                                                   <  1:264196/148‑1 (MQ=255)                                                    cgatcgatAAACTGCCTCAGCAGAAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTg                                                                                             >  7:240745/1‑149 (MQ=255)                                                       tcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGtt                                                                                          <  3:34787/149‑1 (MQ=255)                                                         gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTtgt                                                                                        <  7:321128/149‑1 (MQ=255)                                                                  ccTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCaa                                                                               <  5:368690/149‑1 (MQ=255)                                                                      agcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                             <  2:9030/147‑1 (MQ=255)                                                                      agcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                             <  7:32212/147‑1 (MQ=255)                                                                        cagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                         <  1:313111/149‑1 (MQ=255)                                                                        cagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGAGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                         >  1:325479/1‑149 (MQ=255)                                                                        aggcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                           <  6:82598/145‑1 (MQ=255)                                                                          gcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑TCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGa                                                                       <  7:31957/149‑1 (MQ=255)                                                                           cagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGAc                                                                      >  7:261988/1‑149 (MQ=255)                                                                           cagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGAc                                                                      <  8:180602/149‑1 (MQ=255)                                                                            agTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACt                                                                     >  4:180069/1‑149 (MQ=255)                                                                               aaTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTg                                                                   >  1:152101/1‑148 (MQ=255)                                                                                   ttGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                           <  4:247738/136‑1 (MQ=255)                                                                                   ttGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                           >  3:247738/1‑136 (MQ=255)                                                                                             ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                                    <  6:255580/149‑1 (MQ=255)                                                                                               tGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTa                                                  >  8:290278/1‑149 (MQ=255)                                                                                               tGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTa                                                  >  7:315737/1‑149 (MQ=255)                                                                                                 aTCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                                <  4:38130/149‑1 (MQ=255)                                                                                                     cgACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                         <  4:354584/120‑1 (MQ=255)                                                                                                     cgACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                         >  3:354584/1‑120 (MQ=255)                                                                                                     cgACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACAcc                                            <  8:315737/149‑1 (MQ=255)                                                                                                     cgACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACAcc                                            >  8:332298/1‑149 (MQ=255)                                                                                                       aCGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTg                                                                                             >  5:352507/1‑98 (MQ=255)                                                                                                       aCGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTg                                                                                             <  6:352507/98‑1 (MQ=255)                                                                                                           ttGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATgc                                                                                                         >  4:36270/1‑82 (MQ=255)                                                                                                           ttGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATgc                                                                                                         <  3:36270/82‑1 (MQ=255)                                                                                                            tGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTg                                                                   >  5:150093/1‑119 (MQ=255)                                                                                                            tGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTg                                                                   <  6:150093/119‑1 (MQ=255)                                                                                                               ggcgTGATCCAAGGCCCCCGGTCATGCACATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgggg                                  <  7:332298/149‑1 (MQ=255)                                                                                                                 cgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGaa                                <  8:261988/149‑1 (MQ=255)                                                                                                                  gTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑TCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAg                               <  8:20735/149‑1 (MQ=255)                                                                                                                     aTCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGCGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                            <  2:325479/149‑1 (MQ=255)                                                                                                                        cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctct                                                   >  4:315799/1‑123 (MQ=255)                                                                                                                        cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctct                                                   <  3:315799/123‑1 (MQ=255)                                                                                                                        cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGtt                           <  6:122580/147‑1 (MQ=255)                                                                                                                        cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGtt                           >  5:122580/1‑147 (MQ=255)                                                                                                                        cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAg                         >  1:15930/1‑149 (MQ=255)                                                                                                                              ccccGGCCATGTGCATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg                   >  1:378393/1‑149 (MQ=255)                                                                                                                               cccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg                   >  7:245779/1‑148 (MQ=255)                                                                                                                                ccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg                  >  2:99538/1‑148 (MQ=255)                                                                                                                                ccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg                  >  4:323798/1‑148 (MQ=255)                                                                                                                                ccGGCCATGTGCATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg                  >  1:283678/1‑148 (MQ=255)                                                                                                                                   gTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAg              >  1:279824/1‑149 (MQ=255)                                                                                                                                    ccatgtgcaTGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg             <  2:378393/142‑1 (MQ=255)                                                                                                                                        gCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgtg          >  1:274338/1‑149 (MQ=255)                                                                                                                                        gCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgtgt         >  6:7031/1‑149 (MQ=40)                                                                                                                                             tGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCgg    <  2:279824/149‑1 (MQ=35)                                                                                                                                               gCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCggg   <  7:290278/148‑1 (MQ=34)                                                                                                                                               gCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGt  <  1:99538/149‑1 (MQ=34)                                                                                                                                               gCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGAGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCggg   <  8:245779/148‑1 (MQ=21)                                                                                                                                               |                                                                                                                                                      ACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGC‑GCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGT  >  NC_002947/1425519‑1425809

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