Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	702,991	G→T	6.6%	noncoding (105/120 nt)	PP_5SD →	5S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	702,991	0	G	T	6.6%	260.8 / 7.7	91	noncoding (105/120 nt)	PP_5SD	5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (38/47);  new base T (3/3);  total (41/50)
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.94e-01

ACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCT  >  NC_002947/702912‑703138                                                                                |                                                                                                                                                    accaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGa                                                                               <  1:5440231/150‑1 (MQ=11)  ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaa                                                                              >  2:2316030/1‑150 (MQ=31)   caccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaa                                                                             <  2:517560/150‑1 (MQ=32)   caccTGATCCCATCCCGAACTCAATAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaa                                                                             <  2:4285533/150‑1 (MQ=17)    accTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaaa                                                                            >  2:3237736/1‑150 (MQ=32)     ccTGATCCCATCCCGAACTCAGTAGTGAAAAGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTATGGGTTTTGTCTTTGGTGGGaaaaa                                                                           >  1:169332/1‑150 (MQ=17)       tGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaaaaaa                                                                         >  1:31613/1‑150 (MQ=32)       tGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaaaaaa                                                                         >  1:3375543/1‑150 (MQ=32)        gATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAAc                                                                        <  2:504060/150‑1 (MQ=32)         atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACt                                                                       >  2:2615074/1‑150 (MQ=33)         atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACt                                                                       <  2:2556147/150‑1 (MQ=33)         atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACt                                                                       <  2:2466395/150‑1 (MQ=33)          tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTg                                                                      <  1:5343496/150‑1 (MQ=33)          tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTg                                                                      <  1:1052932/150‑1 (MQ=33)            ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGa                                                                    >  1:2030737/1‑150 (MQ=33)            ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGa                                                                    >  2:2324651/1‑150 (MQ=33)            ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGGGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGa                                                                    >  1:12532/1‑150 (MQ=18)             catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGaa                                                                   >  2:1098830/1‑150 (MQ=33)             catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGGGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGaa                                                                   >  1:4338326/1‑150 (MQ=18)              atcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGGGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAAc                                                                  >  1:947141/1‑150 (MQ=18)               tcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACt                                                                 >  1:4532099/1‑150 (MQ=33)                  cGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGTGGTTTCACCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACACCTATCAGCGCATGAAGGTAGGGGTTTTGTATTTGGTGGGAAAAAAACTGGAAATGCt                                                              >  1:5320719/1‑150 (MQ=11)                   gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAGGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTc                                                             <  1:260585/150‑1 (MQ=18)                   gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTc                                                             <  2:2276785/150‑1 (MQ=33)                   gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTc                                                             <  1:2954270/150‑1 (MQ=33)                     aCTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGa                                                           <  1:353626/150‑1 (MQ=33)                     aCTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGa                                                           <  1:982464/150‑1 (MQ=33)                       tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAAt                                                         >  1:225085/1‑150 (MQ=33)                       tCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAAt                                                         >  1:2251855/1‑150 (MQ=33)                         agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCa                                                       <  1:1834159/150‑1 (MQ=33)                           tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACt                                                     <  2:3859376/150‑1 (MQ=33)                             gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGc                                                   <  2:944923/150‑1 (MQ=33)                              tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGCAAAAAACTGGAACTGCTCGAATCACTGCg                                                  <  1:4218846/150‑1 (MQ=14)                              tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCg                                                  <  2:3144628/150‑1 (MQ=33)                              tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCg                                                  >  2:2346423/1‑150 (MQ=33)                                 aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGa                                               <  2:2396260/150‑1 (MQ=33)                                 aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGa                                               <  2:5125233/150‑1 (MQ=33)                                 aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGa                                               <  1:4437579/150‑1 (MQ=33)                                  aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAg                                              >  2:2380226/1‑150 (MQ=33)                                  aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAg                                              >  2:2380224/1‑150 (MQ=33)                                      tGCATCGCCGATGGGAGTGGGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCAGGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGga                                          <  2:1976259/150‑1 (MQ=16)                                         aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaaaaaa                                                                         <  2:4367266/116‑1 (MQ=32)                                         aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGaaaaaaa                                                                         >  1:4367266/1‑116 (MQ=32)                                         aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGAt                                       >  1:2528895/1‑150 (MQ=33)                                         aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGAt                                       <  1:4122702/150‑1 (MQ=33)                                          tCGGCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTGTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATg                                      >  1:2521700/1‑150 (MQ=18)                                          tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATg                                      <  2:672102/150‑1 (MQ=33)                                          tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATg                                      >  1:2547260/1‑150 (MQ=33)                                              cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCAc                                  <  1:4197831/150‑1 (MQ=33)                                                  ggTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGc                              <  1:1517258/150‑1 (MQ=33)                                                    tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                            <  1:1573335/150‑3 (MQ=18)                                                    tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                            >  2:1573335/1‑148 (MQ=18)                                                    tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                            >  2:1573333/1‑148 (MQ=18)                                                    tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                            >  2:1573234/1‑148 (MQ=18)                                                    tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                            <  1:1573333/150‑3 (MQ=18)                                                    tAGTGTGGGGTTTCCCCATGTGAGAGTATGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGcaa                            <  1:1573234/150‑3 (MQ=18)                                                    tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGt                            >  1:911062/1‑150 (MQ=33)                                                      gtgtGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTg                          <  1:2190286/150‑1 (MQ=33)                                                      gtgtGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTg                          >  1:270353/1‑150 (MQ=33)                                                      gtgtGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTg                          <  1:5224317/150‑1 (MQ=33)                                                         tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCg                       >  1:2064142/1‑150 (MQ=33)                                                         tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCg                       >  2:1977031/1‑150 (MQ=33)                                                         tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCg                       >  2:1976950/1‑150 (MQ=33)                                                         tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCg                       >  2:3478408/1‑150 (MQ=33)                                                         tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCg                       >  1:2318648/1‑150 (MQ=33)                                                         tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTCCGGGAGCGGAGATGGCACCGGCGTTGCCg                       >  1:4780890/1‑150 (MQ=18)                                                          ggggTTTCCGCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCTCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTGTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCgg                      >  2:3517098/1‑150 (MQ=16)                                                           gggTTTCCCCATGTGAGAGTAGGTCATCGTCAATATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACt                                                                       <  1:4505774/100‑1 (MQ=25)                                                           gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAATCt                                                                       >  2:4505774/1‑100 (MQ=18)                                                                  cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCTCCGGCGTTGCCGGTGTTTGCg              >  2:2293838/1‑150 (MQ=25)                                                                  cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGACACCGGCAACGCCGgtgccatct              <  1:5177217/150‑7 (MQ=11)                                                                  cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGACACCGGCAACGCCGgtgccatct              <  1:5177505/150‑7 (MQ=11)                                                                   ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTTGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg             <  1:1327179/150‑1 (MQ=25)                                                                   ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg             <  1:1298131/150‑1 (MQ=35)                                                                   ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg             <  2:76859/150‑1 (MQ=35)                                                                   ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg             <  1:3478408/150‑1 (MQ=35)                                                                   ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg             <  1:909790/150‑1 (MQ=35)                                                                   ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCgg             <  1:2199495/150‑1 (MQ=35)                                                                    cATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCggg            <  2:495844/150‑1 (MQ=35)                                                                      tgtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTa          >  2:140103/1‑150 (MQ=35)                                                                      tgtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTa          >  2:140262/1‑150 (MQ=35)                                                                      tgtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGAAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGTGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTa          >  2:140038/1‑150 (MQ=25)                                                                       gtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTaa         >  2:2865183/1‑150 (MQ=35)                                                                         gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGTCACCGGCGTTGCCGGTGTTTGCGGGTAAAc       <  1:4526834/150‑1 (MQ=25)                                                                          agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAAcc      <  2:911062/150‑1 (MQ=35)                                                                          agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAAcc      <  2:2030737/150‑1 (MQ=35)                                                                             gTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGc   <  2:3872234/150‑1 (MQ=35)                                                                             gTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGc   <  2:2547260/150‑1 (MQ=35)                                                                              tAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCt  <  2:4360223/150‑1 (MQ=35)                                                                              tAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCt  <  2:4360225/150‑1 (MQ=35)                                                                              tAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCt  <  2:2592802/150‑1 (MQ=35)                                                                                |                                                                                                                                                    ACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGGTGGGAAAAAAACTGGAACTGCTCGAATCACTGCGGGAGCGGAGATGGCACCGGCGTTGCCGGTGTTTGCGGGTAAACCCGCT  >  NC_002947/702912‑703138

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