Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	422,059	T→C	5.9%	Q630R (CAG→CGG)	PP_0348 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	422,059	0	T	C	5.9%	246.4 / 4.2	84	Q630R (CAG→CGG)	PP_0348	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/42);  new base C (3/2);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

AGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTAC  >  NC_002947/421912‑422207                                                                                                                                                    |                                                                                                                                                      aGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCTGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTg                                                                                                                                                      <  1:131384/149‑1 (MQ=255)    atatCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGc                                                                                                                                                   >  4:124779/1‑149 (MQ=255)       tCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGa                                                                                                                                                 <  4:276218/148‑1 (MQ=255)         gATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑t                                                                                                                                               >  3:262004/1‑147 (MQ=255)         gATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑Tc                                                                                                                                              <  4:195187/148‑1 (MQ=255)         gATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCa                                                                                                                                             <  7:73990/149‑1 (MQ=255)           ttGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAgc                                                                                                                                           <  3:389078/149‑1 (MQ=255)               cATGCGGCCAATCGACAGCACCAGCAAAGCGCGGATGGGCAGAAGCAGCCTGATGAAATCCAGCAAACGCCAGCTCAGAGCCATGACCGACTCCGCAGTCACTTCCATGGACGTAATGCTGTCGGCGTCGGACTTGGCCA‑TCAGAGccc                                                                                                                                       >  6:226080/1‑149 (MQ=255)                 tGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCgg                                                                                                                                     <  5:99910/149‑1 (MQ=255)                   cGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGc                                                                                                                                   >  8:371890/1‑149 (MQ=255)                   cGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCCGGGCGC‑TGATCGCCCagg                                                                                                                                    <  5:363951/148‑1 (MQ=255)                    ggCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGcc                                                                                                                                  >  6:338896/1‑149 (MQ=255)                      ccAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTc                                                                                                                                >  1:261171/1‑149 (MQ=255)                       cAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCg                                                                                                                               >  2:330802/1‑149 (MQ=255)                        actCGACAGCACCCGCAACGCGCGGCTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCgg                                                                                                                              <  8:300397/147‑1 (MQ=255)                        aaTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCgg                                                                                                                              >  2:76212/1‑149 (MQ=255)                        aaTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGGCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCgg                                                                                                                              <  5:48294/149‑1 (MQ=255)                               aGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgta                                                                                                                       <  1:213276/149‑1 (MQ=255)                               aGCACCAGCAACGCGCGGGTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGGTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgta                                                                                                                       <  6:209178/149‑1 (MQ=255)                                 cACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtag                                                                                                                      >  5:77799/1‑148 (MQ=255)                                 cACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtag                                                                                                                      >  6:170506/1‑148 (MQ=255)                                      gCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCg                                                                                                                <  2:261171/149‑1 (MQ=255)                                            cgcgGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGaaa                                                                                                          <  7:371890/149‑1 (MQ=255)                                             gcgGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCgaagaa                                                                                                           >  6:357459/1‑147 (MQ=255)                                               ggTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                       <  5:338896/149‑1 (MQ=255)                                               ggTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                       >  2:279675/1‑149 (MQ=255)                                               ggTTGGGCAGCAGCAGGATGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCTCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCTA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                       <  5:259167/149‑1 (MQ=255)                                                     gcagcagcaGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCt                                                                                                 >  8:256229/1‑149 (MQ=255)                                                       agcagcagGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCttt                                                                                               <  6:39341/149‑1 (MQ=255)                                                                cTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCggg                                                                                      <  5:170506/149‑1 (MQ=255)                                                                  gatgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                     <  4:359345/132‑1 (MQ=255)                                                                  gatgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                     >  3:359345/1‑132 (MQ=255)                                                                   atgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGt                                                                                   <  5:357459/149‑1 (MQ=255)                                                                    tgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTg                                                                                  >  7:29160/1‑149 (MQ=255)                                                                      aaaTCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTgg                                                                                 <  7:528/148‑1 (MQ=255)                                                                       aaTCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGt                                                                                >  8:365945/1‑148 (MQ=255)                                                                        aTCCAGCAAACGCCTGCTCAGCGCCATGACCGACTCCGCATTCACTTCCGTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCTAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTc                                                                               <  5:226080/148‑1 (MQ=255)                                                                        aTCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTc                                                                               >  5:378088/1‑148 (MQ=255)                                                                          ccAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGaaaa                                                                                                         <  7:30933/120‑1 (MQ=255)                                                                          ccAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGaaaa                                                                                                         >  8:30933/1‑120 (MQ=255)                                                                           cAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAAt                                                                           >  4:329098/1‑149 (MQ=255)                                                                           cAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAAt                                                                           >  7:313245/1‑149 (MQ=255)                                                                           cAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaa                                                                            <  1:76212/148‑1 (MQ=255)                                                                            aGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaa                                                                            <  5:390894/147‑1 (MQ=255)                                                                            aGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaa                                                                            >  6:390894/1‑147 (MQ=255)                                                                             gCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGc                                                                         <  5:219926/149‑1 (MQ=255)                                                                                  cGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGc                                                                     >  1:274771/1‑148 (MQ=255)                                                                                    ccAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCa                                                                  <  2:274771/149‑1 (MQ=255)                                                                                         tCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGgcgc                                                             >  7:88021/1‑149 (MQ=255)                                                                                           aGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGgcgc                                                             <  1:330802/147‑1 (MQ=255)                                                                                            tcgcCATGACCGACTCCGCATTCACTTCGCTGGTCGTCATGCTGTCGGCGGCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGcccc                                                          <  8:346828/148‑1 (MQ=255)                                                                                               ccATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTg                                                       <  7:256229/149‑1 (MQ=255)                                                                                                cATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACTTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGt                                                      >  3:127739/1‑149 (MQ=255)                                                                                                 aTGACCGACTCCGCATTCACTTCCCTGGTCGACATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTc                                                     >  7:298382/1‑149 (MQ=255)                                                                                                  tGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTcc                                                    <  7:365945/149‑1 (MQ=255)                                                                                                  tGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTcc                                                    <  8:313245/149‑1 (MQ=255)                                                                                                     ccGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtagt                                                                                                                     <  6:147911/81‑1 (MQ=255)                                                                                                     ccGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtagt                                                                                                                     >  5:147911/1‑81 (MQ=255)                                                                                                      cGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGc                                                 >  6:377328/1‑148 (MQ=255)                                                                                                          tCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCg                                            >  1:318349/1‑149 (MQ=255)                                                                                                            cGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGc                                                                                     <  7:204825/106‑1 (MQ=255)                                                                                                            cGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGc                                                                                     >  8:204825/1‑106 (MQ=255)                                                                                                                ttCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCAc                                      <  3:124779/149‑1 (MQ=255)                                                                                                                   aCTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGCCCTCCGCCATCTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGCCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCTGCAACCGCTGCACCGt                                   >  5:78585/1‑149 (MQ=255)                                                                                                                    cTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTGAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTg                                  <  1:279675/149‑1 (MQ=255)                                                                                                                    cTTCCCTGGTCGTCATGCTGTCGGCGGCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTg                                  <  8:29160/149‑1 (MQ=255)                                                                                                                     ttCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTgg                                 >  3:46803/1‑149 (MQ=255)                                                                                                                       cccTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGGCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTGGGGCGGGGGCAAATGCGTGCTCAGGCGCCCCTTTTCCAGCAACCGCTGCACCGTGGCg                               >  2:188409/1‑149 (MQ=255)                                                                                                                       cccTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGc                                <  6:164764/148‑1 (MQ=255)                                                                                                                       cccTGGGCGTCATGCTGACGGCGGCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGc                                <  8:273413/148‑1 (MQ=255)                                                                                                                           ggTCGTCATGCTGTCGGCGTCGGCCTGGGCGAGTCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGa                              >  1:47358/1‑147 (MQ=255)                                                                                                                            gtcgtcATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGGCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAg                           <  6:77799/148‑1 (MQ=255)                                                                                                                                   tGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAAt                   >  5:8515/1‑149 (MQ=255)                                                                                                                                     cTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATgg                 <  6:8515/149‑1 (MQ=255)                                                                                                                                      tGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGc                <  3:235763/149‑1 (MQ=255)                                                                                                                                         cggcgtcggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACGTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAgg             >  6:167553/1‑149 (MQ=255)                                                                                                                                            cgtcggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGgcgc          >  2:329589/1‑149 (MQ=255)                                                                                                                                            ccgaggcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTgg                                 <  5:299022/122‑1 (MQ=255)                                                                                                                                            ccgaggcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTgg                                 >  6:299022/5‑126 (MQ=255)                                                                                                                                             gtcggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCt         <  2:318349/149‑1 (MQ=255)                                                                                                                                             cgaggcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGAGCTCAGGCGCCCCTTGTCCCGCAACAGCTGCACCGTGGCGACAGGCAGCATTGGCAGGCGCt         >  4:238413/4‑149 (MQ=255)                                                                                                                                                ggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTgg       <  4:46803/148‑1 (MQ=255)                                                                                                                                                 gcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGc     >  7:219290/1‑149 (MQ=255)                                                                                                                                                   cTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCta   >  6:205116/1‑149 (MQ=255)                                                                                                                                                    tGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCtac  >  3:377633/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      AGGATATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTAC  >  NC_002947/421912‑422207

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