Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,187,586	C→G	6.2%	L31V (CTG→GTG)	PP_4567 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,187,586	0	C	G	6.2%	198.9 / 4.0	65	L31V (CTG→GTG)	PP_4567	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base G (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

ATGGCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAG  >  NC_002947/5187443‑5187730                                                                                                                                                |                                                                                                                                                 aTGGCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGaa                                                                                                                                             <  8:324764/149‑1 (MQ=255)           aaTGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACgag                                                                                                                                    >  7:284321/1‑148 (MQ=255)            aTGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCTGCGCGCAGTTGAAGCGGGGTACTGCATCAACAGCGCGTACAAAAAGGTCCTGTTCCAGGCGGAACAGGAcctttt                                                                                                                                  <  3:338365/149‑6 (MQ=255)               aGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACg                                                                                                                                <  2:205678/148‑1 (MQ=255)                    ttGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTAcc                                                                                                                          <  5:90348/149‑1 (MQ=255)                         gATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCg                                                                                                                     >  2:111087/1‑149 (MQ=255)                              cgcATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTc                                                                                                                >  4:225977/1‑149 (MQ=255)                               gcATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCg                                                                                                               >  6:151769/1‑149 (MQ=255)                               gcATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCg                                                                                                               <  7:58682/149‑1 (MQ=255)                               gcATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCg                                                                                                               >  8:263822/1‑149 (MQ=255)                                 aTCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGt                                                                                                             >  6:267837/1‑149 (MQ=255)                                 aTCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCAGGCGGAACAGGAccttt                                                                                                                                   >  7:285816/1‑123 (MQ=255)                                 aTCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCAGGCGGAACAGGAccttt                                                                                                                                   <  8:285816/127‑5 (MQ=255)                                  tCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTc                                                                                                            >  7:296598/1‑149 (MQ=255)                                   cTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTc                                                                                                            >  8:270314/1‑148 (MQ=255)                                      taCGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAgg                                                                                                        <  1:111087/149‑1 (MQ=255)                                      taCGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAgg                                                                                                        <  7:260740/149‑1 (MQ=255)                                        cggcggCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGc                                                                                                      <  8:284321/149‑1 (MQ=255)                                         ggcggcCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGcc                                                                                                     >  5:191920/1‑149 (MQ=255)                                         ggcggcCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGcc                                                                                                     >  2:231683/1‑149 (MQ=255)                                         ggcggcCTGCTGATTGGCCTGGCCTTGTGACTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGGCGAGAACGGCTACCCACCGGCCTCGGTCGAGGGcc                                                                                                     >  1:24733/1‑149 (MQ=255)                                             cccTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCtgt                                                                                                  <  2:49888/147‑1 (MQ=255)                                                   tgATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGcc                                                                                             >  8:49355/1‑147 (MQ=255)                                                     aTTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGTCGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAg                                                                                          <  2:24733/148‑1 (MQ=255)                                                      ttGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGc                                                                                        >  7:37388/1‑149 (MQ=255)                                                         gcctggcctTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGt                                                                                     <  5:164496/149‑1 (MQ=255)                                                         gcctggcctTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGt                                                                                     <  5:151769/149‑1 (MQ=255)                                                         gcctggcctTGTGCCTGGTGCTGCGCCCCGCGCGCAGATGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTCTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCATCGGTCGAGGGCCTCTGGGCCAAGGCCGt                                                                                     >  5:44833/1‑149 (MQ=255)                                                           ctggcctTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACg                                                                                                                                >  2:261550/1‑104 (MQ=255)                                                           ctggcctTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACg                                                                                                                                <  1:261550/104‑1 (MQ=255)                                                           ctggcctTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTg                                                                                    >  5:289085/1‑148 (MQ=255)                                                            tggcctTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTTCAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTgg                                                                                   >  2:196769/1‑148 (MQ=255)                                                             ggcctTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGaa                                                                                 <  8:296598/149‑1 (MQ=255)                                                                      gCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTAt                                                                        >  7:49185/1‑149 (MQ=255)                                                                       cTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGTTATg                                                                       <  7:49355/149‑1 (MQ=255)                                                                       cTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATg                                                                       <  3:225977/149‑1 (MQ=255)                                                                       cTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATg                                                                       <  5:267837/149‑1 (MQ=255)                                                                          gTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTTCAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTCGGCCAAGGCCGTGGAAGGCGGCTATGcc                                                                     <  1:196769/148‑1 (MQ=255)                                                                           tgctgcGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCAt                                                                   <  7:270314/149‑1 (MQ=255)                                                                               gcgcCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGa                                                                >  3:211955/1‑148 (MQ=255)                                                                                   ccGGCTCGCAGTTGTAGGGGGTCACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGTCTATGCCATCGacaaca                                                           <  6:62063/149‑1 (MQ=255)                                                                                        gcgcAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAgg                                                                                         <  5:265972/114‑1 (MQ=255)                                                                                        gcgcAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAgg                                                                                         >  6:265972/1‑114 (MQ=255)                                                                                           cAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTc                                                   >  1:242717/1‑149 (MQ=255)                                                                                              ttGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCAt                                                >  5:326931/1‑149 (MQ=255)                                                                                              ttGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCAt                                                >  3:43178/1‑149 (MQ=255)                                                                                                        gTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCa                                      <  2:242717/149‑1 (MQ=255)                                                                                                             gCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGc                                                                      >  1:331304/1‑112 (MQ=255)                                                                                                             gCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGc                                                                      <  2:331304/112‑1 (MQ=255)                                                                                                              cATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCTccc                                >  6:257362/1‑149 (MQ=255)                                                                                                               aTGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGcccc                               <  6:289085/149‑1 (MQ=255)                                                                                                                   aCAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGc                           >  5:234518/1‑149 (MQ=255)                                                                                                                         cgTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGAc                        <  5:50146/146‑1 (MQ=255)                                                                                                                         cgTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGAc                        >  6:50146/1‑146 (MQ=255)                                                                                                                            aCAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGAc                        <  3:100263/143‑1 (MQ=255)                                                                                                                            aCAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGAc                        >  4:100263/1‑143 (MQ=255)                                                                                                                             cAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCg                 <  4:211955/149‑1 (MQ=255)                                                                                                                               aaaaGGTCCTGTTCCGCCTGGAACAGGACGAGACCGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGACTATGGCATCGCCCCGGGCGACGTGATCgcc               >  3:21506/1‑147 (MQ=255)                                                                                                                                aaaGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGaca              >  3:208379/1‑149 (MQ=255)                                                                                                                                aaaGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGaca              >  4:316018/1‑149 (MQ=255)                                                                                                                                   ggTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACAccc           <  3:316018/149‑1 (MQ=255)                                                                                                                                     tCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCcgc         >  8:3368/1‑149 (MQ=255)                                                                                                                                      ccTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCcgc         >  3:52446/1‑148 (MQ=255)                                                                                                                                       cTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGa       <  8:49185/149‑1 (MQ=255)                                                                                                                                              gCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAg  >  3:188545/1‑147 (MQ=255)                                                                                                                                                |                                                                                                                                                 ATGGCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAG  >  NC_002947/5187443‑5187730

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