Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	945,383	C→T	5.4%	G80D (GGC→GAC)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	945,383	0	C	T	5.4%	282.7 / 3.1	92	G80D (GGC→GAC)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (40/47);  new base T (2/3);  total (42/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 = 5.88e-01

AGCCCGGCAATCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCTGGAGAAACCGGTACCTG  >  NC_002947/945239‑945531                                                                                                                                                 |                                                                                                                                                      aGCCCGGCAATCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCcgac                                                                                                                                                   >  3:239070/1‑149 (MQ=255) aGCCCGGCAATCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCcgac                                                                                                                                                   >  1:80984/1‑149 (MQ=255) aGCCCGGCAATCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCcgac                                                                                                                                                   >  8:156923/1‑149 (MQ=255)  gCCCGGCAATCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCcgacg                                                                                                                                                  <  7:367132/149‑1 (MQ=255)   cccGGCAATCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCcgacga                                                                                                                                                 <  6:12618/149‑1 (MQ=255)   cccGGCAATCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCcgacg                                                                                                                                                  >  4:385617/1‑148 (MQ=255)     cGGCAATCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGTCGgcg                                                                                                                                                  >  6:279509‑M1/1‑143 (MQ=255)        cAATCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGATCGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGAcc                                                                                                                                            <  6:431581/149‑1 (MQ=255)        cAATCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGAcc                                                                                                                                            <  2:425888/149‑1 (MQ=255)         taTCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGAcaa                                                                                                                                           >  8:293430/2‑147 (MQ=255)         aaTCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGTCGGCGAAGAcgg                                                                                                                                           <  5:279509/149‑3 (MQ=255)         aaTCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGTCGGCGAAGAcgg                                                                                                                                           <  7:43842/149‑3 (MQ=255)             gCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGTCGGCGAAGAcggcac                                                                                                                                        <  1:227603/148‑6 (MQ=255)             gCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATc                                                                                                                                      <  2:305110/149‑1 (MQ=255)                cAAGTTGGCATGGCACTGCCTGCAGTACAGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCt                                                                                                                                   >  1:25223/1‑149 (MQ=255)                 aaGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTg                                                                                                                                  <  2:80984/149‑1 (MQ=255)                   gTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGcc                                                                                                                                <  8:326952/149‑1 (MQ=255)                        catggcaCTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGAccc                                                                                                                           <  6:407583/149‑1 (MQ=255)                         atggcaCTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGAcccc                                                                                                                          >  3:398675/1‑149 (MQ=255)                            gcaCTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCgg                                                                                                                        <  6:8000/148‑1 (MQ=255)                                  cctgcAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTc                                                                                                                 <  3:385617/149‑1 (MQ=255)                                  cctgcAGTACCGATGCCTGCGACCAGTTCGTCCAGCCAGGCCACTTCTTCGAACGTGAACCCGGTGCACGCTGCCGCATCCTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTc                                                                                                                 <  4:244934/149‑1 (MQ=255)                                      cAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAc                                                                                                              >  6:405389‑M2/1‑147 (MQ=255)                                       aGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGtggaccaggccggggtcggcagcgatatggtcttcgtcgg                                                                                                             >  1:375765‑M1/1‑109 (MQ=255)                                       aGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGtggaccaggccggggtcggcagcgatatggtcttcgtcgg                                                                                                             >  7:108112‑M1/1‑109 (MQ=255)                                       aGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcgg                                                                                                            >  5:167805‑M2/1‑146 (MQ=255)                                         tACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcggc                                                                                                           <  6:443963‑M2/148‑5 (MQ=255)                                           ccGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcggcga                                                                                                         >  7:109951‑M2/1‑142 (MQ=255)                                            cGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcggcgaag                                                                                                       >  2:25223‑M2/1‑141 (MQ=255)                                              aTGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCgaagaa                                                                                                     <  4:398675/149‑1 (MQ=255)                                               tGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGtggaccaggccggggtcggcagcgatatggtcttcgtcggcgaagacg                                                                                                     <  8:108112‑M1/149‑49 (MQ=255)                                               tGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGtggaccaggccggggtcggcagcgatatggtcttcgtcggcgaagacg                                                                                                     <  2:375765‑M1/149‑49 (MQ=255)                                                 ccTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGtt                                                                                                  <  1:97922/149‑1 (MQ=255)                                                 ccTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGt                                                                                                   <  7:375644/148‑1 (MQ=255)                                                  cTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTa                                                                                                 >  5:429332/1‑149 (MQ=255)                                                  cTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAAGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTa                                                                                                 >  1:77234/1‑149 (MQ=255)                                                     cGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTAc                                                                                                <  7:156923/147‑1 (MQ=255)                                                     cGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTAc                                                                                                <  2:77234/147‑1 (MQ=255)                                                       aCCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGAt                                                                                            >  4:6673/1‑149 (MQ=255)                                                          aGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCATGGCGAAGAAGTTACCGATGGt                                                                                         <  1:260042/149‑1 (MQ=255)                                                           gTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGt                                                                                         >  1:80856/1‑148 (MQ=255)                                                           gTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGt                                                                                         >  6:344868/1‑148 (MQ=255)                                                               gTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATc                                                                                     <  5:315208/148‑1 (MQ=255)                                                                  cAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGtggaccaggccggggtcggcagcgatatggtcttcgtcggcgaaga                                                                                                       >  7:183394‑M1/1‑82 (MQ=255)                                                                  cAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGtggaccaggccggggtcggcagcgatatggtcttcgtcggcgaaga                                                                                                       <  8:183394‑M1/128‑47 (MQ=255)                                                                   aGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTc                                                                                <  8:276462/149‑1 (MQ=255)                                                                    gCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTcc                                                                               >  7:348370/1‑149 (MQ=255)                                                                       aGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTccccc                                                                            <  3:139373/149‑1 (MQ=255)                                                                         gCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGa                                                                          >  4:77761/1‑149 (MQ=255)                                                                          ccACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGAc                                                                         >  4:265293/1‑149 (MQ=255)                                                                          ccACTTCTTCGAAGGTGAAACCGGTGAACGCTGCAGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCTATGGTGATCTTTTCCCACGAc                                                                         >  4:260052/1‑149 (MQ=255)                                                                           cACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACt                                                                        >  8:4095/1‑149 (MQ=255)                                                                             cttcttCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTc                                                                      <  5:344868/149‑1 (MQ=255)                                                                                cttcGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAg                                                                    <  7:4095/148‑1 (MQ=255)                                                                                   cGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAgg                                                                   <  3:6673/146‑1 (MQ=255)                                                                                    gAAGGTGAAACCTGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGCGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCATGGTCCAAGGCGAAGAAGTTACAGATCGTGATCTTTTCCCCCGACTTCAgggag                                                                >  8:122766/1‑146 (MQ=255)                                                                                    gAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAgg                                                                   <  7:122766/145‑1 (MQ=255)                                                                                    gAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTg                                                                <  3:403651/148‑1 (MQ=255)                                                                                    gAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGa                                                               >  5:166960/1‑149 (MQ=255)                                                                                    gAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGa                                                               <  3:260052/149‑1 (MQ=255)                                                                                    gAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGa                                                               <  3:265293/149‑1 (MQ=255)                                                                                         tGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGtggaccaggccggggtcggcagcgatatggtcttcgtcggcgaaga                                                                                                       <  4:341751‑M1/105‑47 (MQ=255)                                                                                         tGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGtggaccaggccggggtcggcagcgatatggtcttcgtcggcgaaga                                                                                                       >  3:341751‑M1/1‑59 (MQ=255)                                                                                                 gTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCt                                                   <  1:370799/148‑1 (MQ=255)                                                                                                 gTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCt                                                   >  2:370799/1‑148 (MQ=255)                                                                                                  tGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGc                                                 <  7:293430/149‑1 (MQ=255)                                                                                                           ccGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTgg                                        <  8:348370/149‑1 (MQ=255)                                                                                                               aTCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGa                                    <  5:443963/149‑1 (MQ=255)                                                                                                                 cGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGAcc                                  >  3:422395/1‑149 (MQ=255)                                                                                                                      tGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGc                             >  5:325418/1‑149 (MQ=255)                                                                                                                         gCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCAc                          >  6:158568/1‑149 (MQ=255)                                                                                                                              ccACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGt                     <  2:80856/149‑1 (MQ=255)                                                                                                                                acaGGGTGCCGTCTTCGCCGACGAAGACCGATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGAATTCAGGGTGACGATCAGGTCATGCCCGCGCTGGGTGACCGTGGCCACTTTGTc                    >  1:33091/1‑149 (MQ=255)                                                                                                                                acaGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCt                   <  5:372920/149‑1 (MQ=255)                                                                                                                                   gggTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCTGga                <  7:125531/149‑1 (MQ=255)                                                                                                                                   gggTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCTGga                <  4:239070/149‑1 (MQ=255)                                                                                                                                       atggTCTTCGTCGGCGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCTGGAGaaa            >  4:33272/4‑149 (MQ=255)                                                                                                                                        ccGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCTGGAGAAAc           <  4:422395/149‑1 (MQ=255)                                                                                                                                         cGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttcaccttcgaagaagtggcctcgctggacgaactggtcgcaggca          <  8:108112‑M2/149‑102 (MQ=255)                                                                                                                                         cGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttcaccttcgaagaagtggcctcgctggacgaactggtcgcaggca          <  2:375765‑M2/149‑102 (MQ=255)                                                                                                                                         cGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCTGGAGAAAcc          >  1:271010/1‑149 (MQ=255)                                                                                                                                         cGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCTGGAGAAAcc          >  1:346456/1‑149 (MQ=255)                                                                                                                                          gTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCTGGAGAAACCg         <  6:325418/149‑1 (MQ=255)                                                                                                                                           tCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttcaccttcgaagaagtggcctcgctg                             >  7:183394‑M2/1‑46 (MQ=255)                                                                                                                                           tCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttcaccttcgaagaagtggcctcgctg                             <  8:183394‑M2/128‑83 (MQ=255)                                                                                                                                           tCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttca                                                    <  4:341751‑M2/105‑60 (MQ=255)                                                                                                                                           tCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttca                                                    >  3:341751‑M2/1‑46 (MQ=255)                                                                                                                                            cTTCGCCGtggaccaggccggggtcggcagcgatatggtcttcgtcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttcaccttcgaagaagtggcctcgctggacgaactggtcgcaggcatcg        >  2:25223‑M1/1‑8 (MQ=255)                                                                                                                                             ttCGCCGACGAAGCCC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCTGGAGAAACCGGTa      >  6:201238/1‑149 (MQ=255)                                                                                                                                              tCGCCGtggaccaggccggggtcggcagcgatatggtcttcgtcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttcaccttcgaagaagtggcctcgctggacgaactggtcgcaggcatcgg       >  7:109951‑M1/1‑6 (MQ=255)                                                                                                                                                gCCGtggaccaggccggggtcggcagcgatatggtcttcgtcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttcaccttcgaagaagtggcctcgctggacgaactggtcgcaggcatcggta     <  6:443963‑M1/148‑145 (MQ=255)                                                                                                                                                gCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCTGGAGAAACCGGTACCt   >  3:377978/1‑149 (MQ=255)                                                                                                                                                 ccGtggaccaggccggggtcggcagcgatatggtcttcgtcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttcaccttcgaagaagtggcctcgctggacgaactggtcgcaggcatcggtact   >  5:167805‑M1/1‑3 (MQ=255)                                                                                                                                                 ccGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttcaccttcgaagaagtggcctcgctggacgaactggtcgcaggcatcggtact  >  1:375765‑M2/1‑40 (MQ=255)                                                                                                                                                 ccGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCAcggcgaagacggcaccctgtggcatgcccagtacgatgcggcagcgttcaccggtttcaccttcgaagaagtggcctcgctggacgaactggtcgcaggcatcggtact  >  7:108112‑M2/1‑40 (MQ=255)                                                                                                                                                 |                                                                                                                                                      AGCCCGGCAATCGCCCAAGTTGGCATGGCACTGCCTGCAGTACCGATGCCTGCGACCAGTTCGTCCAGCGAGGCCACTTCTTCGAAGGTGAAACCGGTGAACGCTGCCGCATCGTACTGGGCATGCCACAGGGTGCCGTCTTCGCCGACGAAGACC‑ATATCGCTGCCGACCCCGGCCTGGTCCACGGCGAAGAAGTTACCGATGGTGATCTTTTCCCCCGACTTCAGGGTGACGATCAGGTCCTGCCCGCGCTGGGTGACCGTGGCCACTTTGTCTGGAGAAACCGGTACCTG  >  NC_002947/945239‑945531

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