Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	201,361	A→G	7.2%	K2289K (AAA→AAG)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	201,361	0	A	G	7.2%	190.7 / 4.9	69	K2289K (AAA→AAG)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (37/27);  new base G (2/3);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

CAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTG  >  NC_002947/201213‑201509                                                                                                                                                     |                                                                                                                                                      cAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCaaa                                                                                                                                                       <  7:41532/149‑1 (MQ=32)  aGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAAc                                                                                                                                                      <  1:49602/149‑1 (MQ=33)             aGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACt                                                                                                                                           <  5:164854/149‑1 (MQ=38)                  gAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCa                                                                                                                                      >  3:243247/1‑149 (MQ=39)                     aacctgaacctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGt                                                                                                                                   <  1:218976/149‑1 (MQ=39)                      acctgaacctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGt                                                                                                                                   >  2:46161/1‑148 (MQ=39)                           aacctgGACAAGAACACGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGTCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAAt                                                                                                                               <  7:143219/147‑1 (MQ=255)                             cctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGc                                                                                                                           >  5:106236/1‑149 (MQ=255)                              ctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGcc                                                                                                                          <  6:216636/149‑1 (MQ=255)                                 gACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACtgtg                                                                                                                                        >  2:8570/1‑132 (MQ=40)                                 gACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACtgtg                                                                                                                                        <  1:8570/132‑1 (MQ=40)                                   cAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                     >  7:202846/1‑149 (MQ=255)                                   cAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                     <  7:340918/149‑1 (MQ=255)                                    aaGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                     <  3:93368/148‑1 (MQ=255)                                    aaGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                     >  4:93368/1‑148 (MQ=255)                                      gACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGcc                                                                                                                   >  2:172719/1‑148 (MQ=255)                                         cccGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGc                                                                                                                        <  2:328624/140‑1 (MQ=255)                                         cccGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGc                                                                                                                        >  1:328624/1‑140 (MQ=255)                                           cGGGCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGAC                                                                                                             <  8:276955‑M1/149‑1 (MQ=255)                                              tCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCT                                                                                                           >  2:69318‑M1/1‑148 (MQ=255)                                                 gCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATc                                                                                                       <  3:295038/149‑1 (MQ=255)                                                   aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                     <  7:288733/133‑1 (MQ=255)                                                   aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                     >  8:288733/1‑133 (MQ=255)                                                   aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGt                                                                                                     >  6:118323/1‑149 (MQ=255)                                                    cTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGCTCGt                                                                                                     >  6:254275/1‑148 (MQ=255)                                                       ccGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGAcc                                                                                                 >  2:81838/1‑149 (MQ=255)                                                                     ccAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGc                                                                                   >  5:264085/1‑149 (MQ=255)                                                                     ccAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGc                                                                                   >  8:207591/1‑149 (MQ=255)                                                                       aGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCt                                                                                  >  2:1205/1‑148 (MQ=255)                                                                           aCTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAgg                                                                              >  4:165357/1‑148 (MQ=255)                                                                             tCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCa                                                                           >  7:132124/1‑149 (MQ=255)                                                                               atctcaacgAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGggcgttgttgctcagggtcacgaccaggtcgtg                                                                         <  7:296745‑M1/145‑34 (MQ=255)                                                                                gggcAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGggcgttgttgctcagggtcacgaccaggtcgtgg                                                                        <  1:326165‑M1/148‑35 (MQ=255)                                                                                 ggCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCAt                                                                      <  4:243247/149‑1 (MQ=255)                                                                                    aaCGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGggcgttgttgctcagggtcacgaccaggtcgtg                                                                         >  8:296745‑M1/1‑111 (MQ=255)                                                                                    aaCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATc                                                                     >  5:263275/1‑147 (MQ=255)                                                                                    aaCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCaa                                                                   <  6:106236/149‑1 (MQ=255)                                                                                    aaCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCaa                                                                   <  5:118323/149‑1 (MQ=255)                                                                                       gAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGc                                                                >  1:345611/1‑149 (MQ=255)                                                                                               cATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGaaa                                                        <  1:69318/149‑1 (MQ=255)                                                                                                    tCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACACCGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgc                                                   >  6:238995/1‑149 (MQ=255)                                                                                                    tCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑c                                                                         >  3:233601/1‑127 (MQ=255)                                                                                                    tCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑c                                                                         <  4:233601/127‑1 (MQ=255)                                                                                                    tCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgc                                                   <  1:1205/149‑1 (MQ=255)                                                                                                     aaaGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTTAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgcg                                                  >  2:223307/2‑149 (MQ=255)                                                                                                      aaGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgcgc                                                 >  8:39934/1‑149 (MQ=255)                                                                                                       aGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgcgcg                                                >  2:232605/1‑149 (MQ=255)                                                                                                        ggTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgcgcgc                                               >  7:226537/1‑149 (MQ=255)                                                                                                         gTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACAGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGcc                                              >  4:318379/1‑149 (MQ=255)                                                                                                          tCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACACCGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCg                                             >  6:348969/1‑149 (MQ=255)                                                                                                           cACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCACGTCA‑CCATCAAGGCCCGTTAAACCAGCGCGCCg                                             >  2:90534/1‑148 (MQ=255)                                                                                                            acaaTCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACAGACCTTCAATGTGCACGTCAATCAGCCGCTGGACCACGACCTGATCTTGACCCTGAGCAACAACGATCAGGTAA‑CCATCAAGGCCGGTGAAACCAGCGCGCAGTa                                           >  6:270736/4‑149 (MQ=14)                                                                                                               aTCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGggcgttgttgctcagggtcacgaccaggtcgtgg                                                                        <  8:319322‑M1/118‑35 (MQ=255)                                                                                                               aTCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGggcgttgttgctcagggtcacgaccaggtcgtgg                                                                        >  7:319322‑M1/1‑84 (MQ=255)                                                                                                                  aCGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTTAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACAccc                                      <  1:223307/148‑1 (MQ=255)                                                                                                                  aCGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCa                                     <  1:172719/149‑1 (MQ=255)                                                                                                                         aCCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGAcgcg                              <  8:132124/149‑1 (MQ=255)                                                                                                                          ccAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGAcgcg                              >  4:296137/1‑148 (MQ=255)                                                                                                                             gACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCaa                           <  1:81838/149‑1 (MQ=255)                                                                                                                             gACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAg                          >  6:209362/1‑149 (MQ=255)                                                                                                                              aCCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAgg                         <  5:209362/149‑1 (MQ=255)                                                                                                                              aCCTCGGTAGACGAGAACGTCAAACCGACCTTCACTGTTCACTTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCTCGCCGTACACCCACGACGCGCAAgg                         <  5:270736/149‑1 (MQ=255)                                                                                                                                cTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGc                        >  3:321652/1‑148 (MQ=255)                                                                                                                                   ggTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACCCCCACGACGCGCAAGGCGATg                    >  8:221258/1‑149 (MQ=255)                                                                                                                                       gCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGt                <  2:345611/149‑1 (MQ=39)                                                                                                                                        ccGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACACTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTc               >  7:258567/1‑149 (MQ=255)                                                                                                                                         cGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCt              <  8:226537/149‑1 (MQ=39)                                                                                                                                            gaACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATc           <  1:46161/149‑1 (MQ=40)                                                                                                                                                gTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCCTCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGa       >  3:98365/1‑149 (MQ=255)                                                                                                                                                gTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACTACGCGCAACGCCATCACGTCTATCAGGa       >  4:54486/1‑149 (MQ=255)                                                                                                                                                     aCCGACCTTCACTCTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTg  <  6:311777/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTG  >  NC_002947/201213‑201509

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