Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,605,513	C→G	5.4%	A310G (GCC→GGC)	PP_4077 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,605,513	0	C	G	5.4%	239.8 / 3.6	74	A310G (GCC→GGC)	PP_4077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/38);  new base G (2/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

CCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTCG  >  NC_002947/4605367‑4605653                                                                                                                                                   |                                                                                                                                             cccACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGccc                                                                                                                                            >  4:165767/1‑149 (MQ=255) cccACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGcc                                                                                                                                             >  4:111406/1‑148 (MQ=255)     caGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGt                                                                                                                                        <  8:121313/149‑1 (MQ=255)        gTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGTCGa                                                                                                                                     <  6:93957/149‑1 (MQ=255)         tGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAg                                                                                                                                    <  7:149733/149‑1 (MQ=255)          gggCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGc                                                                                                                                   >  4:182807/1‑149 (MQ=255)            gCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTg                                                                                                                                 <  6:123736/149‑1 (MQ=255)             cATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGa                                                                                                                                <  3:165767/149‑1 (MQ=255)             cATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGa                                                                                                                                <  3:111406/149‑1 (MQ=255)              aTCGAGTTCCATTTCCCGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGaa                                                                                                                               <  7:137390/149‑1 (MQ=255)              aTCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTg                                                                                                                                 >  3:188703/1‑147 (MQ=255)              aTCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGaa                                                                                                                               <  1:159744/149‑1 (MQ=255)               tCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAg                                                                                                                              >  3:99099/1‑149 (MQ=255)                    ttCCATTTCACGCGTCAACTGCCAGACACCCTCCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGcc                                                                                                                            >  4:133555/1‑146 (MQ=255)                      ccATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAg                                                                                                                         >  1:117359/1‑147 (MQ=255)                      ccATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGaa                                                                                                                       >  7:57430/1‑149 (MQ=255)                           tctcgggTCAACTGCCAGACACCCTTCGCGTTGGGAGAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTGAAGCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGa                                                                                                                   <  3:300884/142‑1 (MQ=255)                            cACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAg                                                                                                                  <  1:287326/148‑1 (MQ=255)                               gcgTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgcg                                                                                                               >  7:278639/1‑148 (MQ=255)                                cgTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCa                                                                                                             <  2:117359/149‑1 (MQ=255)                                      cTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGACGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGc                                                                                                       >  1:150068/1‑149 (MQ=255)                                        gCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGccc                                                                                                                           >  2:102963/1‑127 (MQ=255)                                        gCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGccc                                                                                                                           <  1:102963/127‑1 (MQ=255)                                         ccAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTa                                                                                                    <  4:99099/149‑1 (MQ=255)                                         ccAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGAGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTa                                                                                                    >  3:216237/1‑149 (MQ=255)                                          cAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTAc                                                                                                   >  2:91727/1‑149 (MQ=255)                                           aGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTAc                                                                                                   >  7:121578/1‑148 (MQ=255)                                                 ccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCt                                                                                            >  3:310097/1‑149 (MQ=255)                                                    tCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGATGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCtt                                                                                           >  3:189343/1‑147 (MQ=255)                                                     cgcgTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTc                                                                                          <  4:310097/147‑1 (MQ=255)                                                     cgcgTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCa                                                                                         <  1:295882/148‑1 (MQ=255)                                                      gcgTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAgt                                                                                       <  2:150068/149‑1 (MQ=255)                                                      gcgTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAgt                                                                                       <  8:240363/149‑1 (MQ=255)                                                           ggCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCg                                                                                   <  8:158110/148‑1 (MQ=255)                                                               caGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGaa                                                                              <  1:91727/149‑1 (MQ=255)                                                               caGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGaa                                                                              <  3:182807/149‑1 (MQ=255)                                                                agaagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAg                                                                             >  4:279928/1‑149 (MQ=255)                                                                 gaagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGCAATCTTCAGTGTCGATGAAGt                                                                            >  1:321831/1‑149 (MQ=255)                                                                 gaagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAgt                                                                                       >  1:263396/1‑138 (MQ=255)                                                                 gaagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAgt                                                                                       <  2:263396/138‑1 (MQ=255)                                                                    gaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcat                                                                         <  6:187769/149‑1 (MQ=255)                                                                     aTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcat                                                                         >  7:139814/1‑148 (MQ=255)                                                                           gCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcat                                                                         <  2:50991/142‑1 (MQ=255)                                                                           gCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcat                                                                         >  1:50991/1‑142 (MQ=255)                                                                            ccTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAccc                                                                 <  4:188703/149‑1 (MQ=255)                                                                            ccTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAccc                                                                 <  3:279928/149‑1 (MQ=255)                                                                             cTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGaa                                                                                                >  1:144636/1‑117 (MQ=255)                                                                             cTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGaa                                                                                                <  2:144636/117‑1 (MQ=255)                                                                              tGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGCAGCCCAGAAGCGAGGCGCAGTCCGCGTAGGAACTCTTCAGTGTCGATGAAGTCATCAGCACCCg                                                                >  5:84849/1‑148 (MQ=255)                                                                                ttCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGa                                                               <  8:278639/147‑1 (MQ=255)                                                                                ttCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGa                                                               <  4:216237/147‑1 (MQ=255)                                                                                     cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaacaac                                                          <  8:121578/147‑1 (MQ=255)                                                                                      aCCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCAc                                                       >  2:15629/1‑149 (MQ=255)                                                                                       ccccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACt                                                      >  3:285911/1‑149 (MQ=255)                                                                                        cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGTCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaaca                                                            <  5:58536/142‑1 (MQ=255)                                                                                        cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaaca                                                            >  6:58536/1‑142 (MQ=255)                                                                                                tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCGCCCGAACCACCACTGACGGAAGt                                             >  1:37343/1‑149 (MQ=255)                                                                                                  aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTAc                                           >  8:135593/1‑149 (MQ=255)                                                                                                    ccTGTTCAAGCACAGTGCCGAGCCCATCGACGTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATGAGCACCCGAACAACCACTGACGGAAGTACCg                                         <  3:41705/149‑1 (MQ=255)                                                                                                    ccTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCg                                         <  3:133555/149‑1 (MQ=255)                                                                                                        ttCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACGACTGACGGAAGTACCGGt                                       <  6:314816/147‑1 (MQ=255)                                                                                                         tCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGa                                     <  7:135593/148‑1 (MQ=255)                                                                                                           aaGCACAGTGCCGAGCCCATCGACCTGGCCAGCGGGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTc                                                                                          >  4:163211/1‑93 (MQ=255)                                                                                                           aaGCACAGTGCCGAGCCCATCGACCTGGCCAGCGCGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTc                                                                                          <  3:163211/93‑1 (MQ=255)                                                                                                              cacaGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCt                               >  2:94754/1‑149 (MQ=255)                                                                                                                 aGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCg                            >  2:153263/1‑149 (MQ=255)                                                                                                                     ccGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGCAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGGCATCAGCCCCCGAACCACCACTGACGGAAGTACCGGTGAGCACCTGCGGACt                        >  4:49707/1‑149 (MQ=255)                                                                                                                          cccATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCg                   >  6:221378/1‑149 (MQ=255)                                                                                                                           ccATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGc                  <  1:126304/149‑1 (MQ=255)                                                                                                                            gctCGACCTGGGCAGCCTCGCTGGCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACgg                                                 >  7:198327/3‑117 (MQ=38)                                                                                                                            gctCGACCTGGGCAGCCTCGCTGGCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACgg                                                 <  8:198327/115‑1 (MQ=38)                                                                                                                               cGACCTGGGCAGCCTCGCTGGCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCtt              >  2:200738/1‑149 (MQ=255)                                                                                                                                gACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCttt             <  5:221378/149‑1 (MQ=255)                                                                                                                                          gCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTg                              <  8:101359/122‑1 (MQ=255)                                                                                                                                          gCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTg                              >  7:101359/1‑122 (MQ=255)                                                                                                                                           cGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTCg  <  4:189343/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             CCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTCG  >  NC_002947/4605367‑4605653

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