Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,271,667	A→G	7.9%	intergenic (+51/‑49)	PP_t32 → / → PP_t33	tRNA‑Asp/tRNA‑Asp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,271,667	0	A	G	7.9%	251.3 / 8.3	88	intergenic (+51/‑49)	PP_t32/PP_t33	tRNA‑Asp/tRNA‑Asp
Reads supporting (aligned to +/- strand):  ref base A (38/43);  new base G (4/3);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCC  >  NC_002947/2271520‑2271810                                                                                                                                                      |                                                                                                                                                 gCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACAt                                                                                                                                                 <  5:14873/149‑1 (MQ=255) gCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACAt                                                                                                                                                 <  4:158701/149‑1 (MQ=255)   cAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGc                                                                                                                                               >  3:167168/1‑149 (MQ=255)       gagGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTg                                                                                                                                           <  5:296244/149‑1 (MQ=255)        agGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCATGTGCCTTgt                                                                                                                                          <  6:150551/149‑2 (MQ=255)         ggAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGAGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCATGTGCCTTgtg                                                                                                                                         <  5:173086/149‑3 (MQ=255)          gAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGc                                                                                                                                       >  3:258634/1‑149 (MQ=255)                aCTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTTACAGAT‑GGGTGTCTGGCAC‑CGGTCTGAAGCCACAAGGCATGTGCCTTgtg                                                                                                                                         >  4:50519/1‑140 (MQ=255)                aCTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGTCTGGCAC‑CGGTCTGAAGCCACAAGGCATGTGCCTTgtg                                                                                                                                         <  3:50519/142‑3 (MQ=255)                 cTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATccc                                                                                                                                <  4:258634/149‑1 (MQ=255)                    cgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGt                                                                                                                             <  7:237012/149‑1 (MQ=255)                     gcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGtt                                                                                                                            <  7:69789/149‑1 (MQ=255)                         cGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCaa                                                                                                                        <  2:226409/149‑1 (MQ=255)                               ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGgc                                                                                                                  <  1:15834/149‑1 (MQ=255)                               ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGgc                                                                                                                  <  4:167168/149‑1 (MQ=255)                                tCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATc                                                                                                                      <  8:144647/144‑1 (MQ=255)                                tCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATc                                                                                                                      >  7:144647/1‑144 (MQ=255)                                   gTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAATCCACAAGGCACAAGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCaa                                                                                                              <  6:8611/149‑1 (MQ=255)                                   gTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGTCACAAGGCACATGCCTTG‑GGCTGATCCCa                                                                                                                               <  5:41370/132‑1 (MQ=255)                                   gTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGTCACAAGGCACATGCCTTG‑GGCTGATCCCa                                                                                                                               >  6:41370/1‑132 (MQ=255)                                   gTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGGCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCaa                                                                                                              <  5:310266/149‑1 (MQ=255)                                           gAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                       <  8:257957/148‑1 (MQ=255)                                           gAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                       <  5:183434/148‑1 (MQ=255)                                           gAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                       <  3:18440/148‑1 (MQ=255)                                              tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGa                                                                                                    >  7:239576/1‑148 (MQ=255)                                              tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGAt                                                                                                   <  5:133034/149‑1 (MQ=255)                                                ccGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATAc                                                                                                 >  3:56591/1‑149 (MQ=255)                                                  ggCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGc                                                                                               <  4:56591/149‑1 (MQ=255)                                                     cTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGc                                                                                            <  8:225330/149‑1 (MQ=255)                                                        tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCg                                                                                           <  6:132844/147‑1 (MQ=255)                                                        tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGt                                                                                         >  7:153856/1‑149 (MQ=255)                                                             ccGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTc                                                                                    <  8:153856/149‑1 (MQ=255)                                                                      gcgGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGt                                                                             >  1:238437/1‑148 (MQ=255)                                                                           ttCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATa                                                                      <  8:239576/149‑1 (MQ=255)                                                                              gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGTTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCa                                                                                   >  5:200494/1‑133 (MQ=255)                                                                              gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGTTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCa                                                                                   <  6:200494/133‑1 (MQ=255)                                                                              gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCg                                                                   <  2:238437/149‑1 (MQ=255)                                                                                 tCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCACTCGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGc                                                                 >  1:239332/1‑149 (MQ=255)                                                                                   ccGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTTGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGc                                                                                               <  4:317711/116‑1 (MQ=255)                                                                                   ccGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGc                                                                                               >  3:317711/1‑116 (MQ=255)                                                                                    cGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCgg                                                                                                                    <  5:218136/94‑1 (MQ=255)                                                                                    cGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCgg                                                                                                                    >  6:218136/1‑94 (MQ=255)                                                                                           tGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCg                                                      >  2:207851/1‑149 (MQ=255)                                                                                            gcgcCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCgg                                                     >  2:260074/1‑149 (MQ=255)                                                                                                cATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGt                                                 >  7:190027/1‑149 (MQ=255)                                                                                                cATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGt                                                 >  2:249785/1‑149 (MQ=255)                                                                                                  tttttCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTTAACGAGACGCAGCGGTAGTTAAGTCGGTTAGAATACAGGCCTGTCACGCCGGCGGTcg                                               >  3:70/1‑149 (MQ=255)                                                                                                        aGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAg                                                                                             >  7:164365/1‑97 (MQ=255)                                                                                                        aGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAg                                                                                             <  8:159875/97‑1 (MQ=255)                                                                                                        aGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAg                                                                                             >  7:159875/1‑97 (MQ=255)                                                                                                        aGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAg                                                                                             <  8:164365/97‑1 (MQ=255)                                                                                                        aGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGtt                                         >  3:201578/1‑149 (MQ=255)                                                                                                             cAGAT‑GGGTGCCTGGCAC‑CGGTATGAAGCCACAAGGCACATGCATTG‑GGCTGAACCAAGTTTCAATCGGGCGCAAGAAAGAAAGATACGCACCGGTAGTTAAGTCGGTTAGAAAACCGGACTCACACGCCGGGGGTCGCGGATTCGAga                                    >  7:205789/1‑148 (MQ=255)                                                                                                              aGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCgg                                                     <  1:292395/131‑1 (MQ=255)                                                                                                              aGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTc                                   >  2:84956/1‑149 (MQ=255)                                                                                                              aGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTc                                   <  2:239332/149‑1 (MQ=255)                                                                                                              aGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCACAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCgg                                                     >  2:292395/1‑131 (MQ=255)                                                                                                                aT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTg‑gg                                                                                                                                        >  3:291722/1‑46 (MQ=255)                                                                                                                aT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTg‑gg                                                                                                                                        <  4:291722/46‑1 (MQ=255)                                                                                                                aT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTccc                                 >  1:33769/1‑149 (MQ=255)                                                                                                                   gggTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCt                                                                                                          >  5:58275/1‑74 (MQ=255)                                                                                                                   gggTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCt                                                                                                          <  6:58275/74‑1 (MQ=255)                                                                                                                         cTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTg                         >  3:267925/1‑149 (MQ=255)                                                                                                                         cTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACCCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTg                         >  5:53865/1‑149 (MQ=255)                                                                                                                          tGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGc                           <  3:103976/146‑1 (MQ=255)                                                                                                                          tGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGc                           >  4:103976/1‑146 (MQ=255)                                                                                                                           ggCAC‑CGGTCTGAAGCCTCAAGGCACATGCCTTG‑GGCTGATCTAAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTgc                        <  4:70/148‑1 (MQ=255)                                                                                                                           ggCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTgcg                       >  4:137211/1‑149 (MQ=255)                                                                                                                            gCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTcgc                                              <  3:215271/125‑1 (MQ=255)                                                                                                                            gCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTcgc                                              >  4:215271/1‑125 (MQ=255)                                                                                                                            gCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTc                                        <  4:197081/131‑1 (MQ=255)                                                                                                                            gCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTgc                        <  2:33769/147‑1 (MQ=255)                                                                                                                              aC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCa                    >  1:122448/1‑149 (MQ=255)                                                                                                                                 cGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATTGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTAAGAATACCGGCCTGGCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCTTc                  >  3:249963/1‑149 (MQ=255)                                                                                                                                   gTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCtt                <  8:190027/149‑1 (MQ=255)                                                                                                                                   gTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCtt                >  6:38270/1‑149 (MQ=255)                                                                                                                                     cTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAAAGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCgg                                                                  <  2:259044/97‑1 (MQ=255)                                                                                                                                     cTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAAAGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCgg                                                                  >  1:259044/1‑97 (MQ=255)                                                                                                                                        aaGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGt             <  4:201578/147‑1 (MQ=255)                                                                                                                                           gcccaAGGCATGTGCCTTGTGGCTGATCCCAGTTACAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc          >  8:93315/4‑148 (MQ=14)                                                                                                                                           ccACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCaa        <  1:249785/149‑1 (MQ=255)                                                                                                                                            cccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTg                         >  3:114337/3‑132 (MQ=14)                                                                                                                                            cccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc          >  3:267628/3‑147 (MQ=255)                                                                                                                                             acaAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAgg      >  7:149979/1‑149 (MQ=255)                                                                                                                                             acaAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAgg      >  5:77237/1‑149 (MQ=255)                                                                                                                                             acaAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCACGTCCGCTGCGCCATCTTCGTTTCAAgg      >  1:307439/1‑149 (MQ=255)                                                                                                                                                 ggCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGGTGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc  <  3:137211/149‑1 (MQ=255)                                                                                                                                                 ggCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc  <  1:260074/149‑1 (MQ=255)                                                                                                                                                    acaTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAATGATACGCAGCGGTAGTTCAGTCGGTTAGAATAcc                                                                    <  2:285826/80‑1 (MQ=255)                                                                                                                                                    acaTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATAcc                                                                    >  1:285826/1‑80 (MQ=255)                                                                                                                                                      aTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGt             >  8:202668/1‑133 (MQ=255)                                                                                                                                                      |                                                                                                                                                 GCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGAT‑GGGTGCCTGGCAC‑CGGTCTGAAGCCACAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCC  >  NC_002947/2271520‑2271810

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