Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	905,474	A→T	5.3%	E578V (GAA→GTA) ‡	PP_0788 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	905,474	0	A	T	5.3%	226.5 / 2.2	75	E578V (GAA→GTA) ‡	PP_0788	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base T (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

ACATGCTGCTGGCCAGCATCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGCCCCTGCACACTGCCGACGGGCAG  >  NC_002947/905330‑905622                                                                                                                                                 |                                                                                                                                                     aCATGCTGCTGGCCAGCATCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAg                                                                                                                                                  >  8:95267/1‑149 (MQ=255)   aTGCTGCTGGCCAGCATCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGtt                                                                                                                                                <  6:189449/149‑1 (MQ=255)      ctgctgGCCAGCATCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGTTCCAGTTGCa                                                                                                                                             <  6:247361/149‑1 (MQ=255)            gCCAGCATCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTgc                                                                                                                                        <  7:95267/148‑1 (MQ=255)            gCCAGCATCACCCATCCACATGATCCGCAAGCCTTGCTGCAGCAGTTGCCAGGCCTGCTCAAGGCGCTTCGCGATGTCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTgc                                                                                                                                        <  2:205384/148‑1 (MQ=255)              cAGCATCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGTTCCAGTTGCa                                                                                                                                             >  6:309455/1‑141 (MQ=255)              cAGCATCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGTTCCAGTTGCa                                                                                                                                             <  5:309455/141‑1 (MQ=255)              cAGCATCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTgcgcg                                                                                                                                     >  8:29858/1‑149 (MQ=255)               aGCATCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTgcgcgc                                                                                                                                    >  1:102609/1‑149 (MQ=255)                  aTCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTg                                                                                                                                 >  6:320510/1‑149 (MQ=255)                   tCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTgc                                                                                                                                >  7:219202/1‑149 (MQ=255)                       cccGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAg                                                                                                                            <  7:29858/149‑1 (MQ=255)                        ccGCCACATGAGCCGCAAGCCTTGCAGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCATGCGCAgg                                                                                                                           >  6:297430/1‑149 (MQ=255)                           ccACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCGCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCAGGCGCGCAGGtgc                                                                                                                        <  7:118275/149‑3 (MQ=255)                                    ccGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCt                                                                                                               >  2:183768/1‑149 (MQ=255)                                      gCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGc                                                                                                             >  3:23201/1‑149 (MQ=255)                                            cTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCg                                                                                                       >  4:143420/1‑149 (MQ=255)                                            cTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCg                                                                                                       >  2:102592/1‑149 (MQ=255)                                            cTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCg                                                                                                       >  4:2228/1‑149 (MQ=255)                                                 tgcAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAg                                                                                                                                                  >  6:200771/1‑101 (MQ=255)                                                 tgcAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAg                                                                                                                                                  <  5:200771/101‑1 (MQ=255)                                                   cagcagGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAg                                                                                                  <  8:44074/147‑1 (MQ=255)                                                    agcagGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCTACCCGCGAGTTCTTCCTGCAACTTGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGc                                                                                               <  4:262510/149‑1 (MQ=255)                                                      cagGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTTCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGcc                                                                                              >  1:288564/1‑148 (MQ=255)                                                          tGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATcc                                                                                         >  2:320089/1‑149 (MQ=255)                                                            ccAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCtt                                                                                       <  1:102592/149‑1 (MQ=255)                                                            ccAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCt                                                                                        <  8:177842/148‑1 (MQ=255)                                                             cAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTg                                                                                      <  7:149781/149‑1 (MQ=255)                                                               ggCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTggg                                                                                    <  5:112990/149‑1 (MQ=255)                                                                 ccTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGAGCGCCTGCGCAGGTAAAGAGGCCATGACGAGCGGGGAGGGCCAAACCTTGGGc                                                                                   >  3:4751/1‑148 (MQ=255)                                                                     cTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGAtgt                                                                              >  3:44251/1‑149 (MQ=255)                                                                     cTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGAtgt                                                                              >  3:251591/1‑149 (MQ=255)                                                                     cTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGAt                                                                                <  3:104925/147‑1 (MQ=255)                                                                     cTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGAt                                                                                <  7:62013/147‑1 (MQ=255)                                                                     cTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGAt                                                                                >  4:104925/1‑147 (MQ=255)                                                                        aaGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCt                                                                                                                                          >  3:182154/1‑86 (MQ=255)                                                                        aaGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCt                                                                                                                                          <  4:182154/86‑1 (MQ=255)                                                                         aGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTg                                                                          <  4:44251/149‑1 (MQ=255)                                                                           gcgcTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGt                                                                        <  2:328028/149‑1 (MQ=255)                                                                           gcgcTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGt                                                                        <  5:320510/149‑1 (MQ=255)                                                                             gcTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTgg                                                                         >  6:199406/1‑146 (MQ=255)                                                                                 cgcgATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGaa                                                                    <  3:261798/147‑1 (MQ=255)                                                                                   cgATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTTCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAAc                                                                   >  8:270810/1‑146 (MQ=255)                                                                                   cgATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAag                                                                <  1:320089/149‑1 (MQ=255)                                                                                   cgATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAag                                                                <  5:199406/149‑1 (MQ=255)                                                                                    gATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTTCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAaga                                                               <  7:270810/149‑1 (MQ=255)                                                                                    gATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACaa                                                                 >  3:18815/1‑147 (MQ=255)                                                                                      tGGCCTGGCAATCGTAGCTCTGGACGCGGTGGCGACCCTCGAGTTCTTCCTTCAACTGGAACAGTTGCACCTGCGCGCCTGCGAAGTTAAAGAGTCGATTCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGATCAAgaga                                                             <  5:297430/149‑1 (MQ=255)                                                                                          cTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTa                                                                       <  4:67190/135‑1 (MQ=255)                                                                                          cTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTa                                                                       >  3:67190/1‑135 (MQ=255)                                                                                          cTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATtgt                                                         >  6:29326/1‑149 (MQ=255)                                                                                             gCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATtgtg                                                        <  6:247678/147‑1 (MQ=255)                                                                                             gCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCt                                                      >  7:1797/1‑149 (MQ=255)                                                                                               aaGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCACTCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGa                                                    >  8:245183/1‑149 (MQ=255)                                                                                               aaGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGa                                                    <  3:2228/149‑1 (MQ=255)                                                                                                aGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTGTTCCTGCAACTGGAACAGTTTCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGAc                                                   <  4:267811/149‑1 (MQ=255)                                                                                                 gCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAAt                                                                                           <  1:111049/108‑1 (MQ=255)                                                                                                 gCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAAt                                                                                           >  2:111049/1‑108 (MQ=255)                                                                                                    tAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGAtt                                               <  2:92413/149‑1 (MQ=255)                                                                                                     aGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTg                                              >  1:146518/1‑149 (MQ=255)                                                                                                        tctGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGcc                                            >  7:67517/1‑148 (MQ=255)                                                                                                            gACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCgaa                                         <  8:219202/147‑1 (MQ=255)                                                                                                                     gCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTg                              >  4:24662/1‑149 (MQ=255)                                                                                                                       gACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTggc                            >  8:261159/1‑147 (MQ=255)                                                                                                                       gACCCGCGAGTTCTTACTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCg                            >  2:174100/1‑149 (MQ=255)                                                                                                                           cgcgAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGcc                        >  2:44985/1‑149 (MQ=255)                                                                                                                            gcgAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATtg                                                          >  1:113449/1‑114 (MQ=255)                                                                                                                            gcgAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATtg                                                          <  2:113449/114‑1 (MQ=255)                                                                                                                                gTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGCCCCTGc                   <  4:23201/149‑1 (MQ=255)                                                                                                                                   cttcCTGCAAGTGTAACAGTTGCACCTGCGCGCCTGCGAAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGCCCCTGcac                 <  4:4751/148‑1 (MQ=255)                                                                                                                                   cttcCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGCCCCTGcac                 <  3:24662/148‑1 (MQ=255)                                                                                                                                        tGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGCCCATGCACACTGcc           >  5:252620/1‑149 (MQ=255)                                                                                                                                         gCAACTGTTCCAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGCCCCTGCACACTGcc           >  1:153595/1‑148 (MQ=255)                                                                                                                                          cAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGCCCCTGCACACTGCCGa         >  4:243443/1‑149 (MQ=255)                                                                                                                                                gAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGCCCCTGCACACTGCCGACggg     >  6:145133/1‑147 (MQ=255)                                                                                                                                                 aaCAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGCCCCTGCACACTGCCGACGGGCAg  >  3:301374/1‑149 (MQ=255)                                                                                                                                                 aaCAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGCCCCTGCACACTGCCGACGGGCAg  >  8:45196/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     ACATGCTGCTGGCCAGCATCACCCCGCCACATGAGCCGCAAGCCTTGCTGCAGCAGGTGCCAGGCCTGCTCAAGGCGCTGCGCGATGGCCTGGCAAGCGTAGCTCTGGACGCGGTGGCGACCCGCGAGTTCTTCCTGCAACTGGAACAGTTGCACCTGCGCGCCTGCGCAGGTAAAGAGGCGCTGCCGAGCGGGGAGGGCCAATCCTTGGGCGATGTGCTGGTAGAACAAGAGATTGTGCTGACGATTGCCGAAGAACCCGCCTGCGTGCCCCTGCACACTGCCGACGGGCAG  >  NC_002947/905330‑905622

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