Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	944,826	A→T	6.0%	S266T (TCG→ACG)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	944,826	0	A	T	6.0%	235.0 / 5.2	83	S266T (TCG→ACG)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (41/36);  new base T (2/3);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

GATATTGCCAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAG  >  NC_002947/944680‑944970                                                                                                                                                   |                                                                                                                                                  gaTATTGCCAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑c                                                                                                                                                >  6:337150/1‑149 (MQ=255)     ttGCCAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGTG‑CCgag                                                                                                                                            >  5:49070/1‑147 (MQ=255)     ttGCCAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCggg                                                                                                                                            >  4:195626/1‑149 (MQ=255)     ttGCCAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCggg                                                                                                                                            >  2:356695/1‑149 (MQ=255)     ttGCCAGCCGAGTCTGCGGCGCTGATATCGAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGTGCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCTCACGGAAATCGTCGAG‑CCgg                                                                                                                                             <  8:15633/148‑1 (MQ=255)        ccAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGTG‑ccaagccc                                                                                                                                         >  5:253133/1‑143 (MQ=255)         cAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCg                                                                                                                                        >  6:222497/1‑149 (MQ=255)         cAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCg                                                                                                                                        <  1:86335/149‑1 (MQ=255)          aGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCgg                                                                                                                                       <  1:356695/149‑1 (MQ=255)               cgTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGTG‑CCAAGCCCGGCTCGa                                                                                                                                  <  6:253133/149‑1 (MQ=255)                gTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGTG‑CCGAGCCCGGCTCGAc                                                                                                                                 <  6:49070/149‑1 (MQ=255)                   tGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGGAACCGTCGAG‑CCGGGCTCGGCACGACCgg                                                                                                                              >  2:96848/1‑149 (MQ=255)                   tGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCgg                                                                                                                              >  3:58350/1‑149 (MQ=255)                     cggcgCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGTTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGTG‑CCGAGCCCGGCTCGAcggt                                                                                                                              <  6:105427/147‑4 (MQ=255)                      ggcgCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTaa                                                                                                                           <  5:80253/149‑1 (MQ=255)                       gcgcTGATCTCCAGGGCCTGCCCGTCGCTTTGTGGCGGTGGCAGGGTAAGGCTCAACTGGCCGGCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGCAACCCGGCCAG‑CCGGGCGCGGCACGACCGGTAAg                                                                                                                          <  3:195626/149‑1 (MQ=255)                         gcTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCg                                                                                                                               >  2:321121/1‑142 (MQ=255)                         gcTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCg                                                                                                                               <  1:321121/142‑1 (MQ=255)                           tGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGgtg                                                                                                                      >  1:50848/1‑149 (MQ=255)                           tGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGgtg                                                                                                                      >  2:162775/1‑149 (MQ=255)                              tctcCCGGGCCTGGCCGTCGCTTTGTGGCGGGTGCAGGGTAATGCTGAGCTGGCCGTCCGCACCGGCCCCCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGCG‑CCGGGCTCGGCACGAGCGGTAAGGgggt                                                                                                                     >  8:204435/1‑147 (MQ=255)                                 ccAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTAccc                                                                                                                >  7:29458/1‑149 (MQ=255)                                  cAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGATACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTAccc                                                                                                                <  1:323432/148‑1 (MQ=255)                                  cAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTa                                                                                                                   <  2:53140/145‑1 (MQ=255)                                  cAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCt                                                                                                               >  3:335586/1‑149 (MQ=255)                                   aGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTg                                                                                                              >  1:344865/1‑149 (MQ=255)                                    gggCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGc                                                                                                             >  1:315065/1‑149 (MQ=255)                                       ccTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCa                                                                                                          <  6:53085/149‑1 (MQ=255)                                        cTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCAt                                                                                                         >  3:152950/1‑149 (MQ=255)                                        cTGGCCGTCGCTTTGTGGCGGTTGAAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGCGCTCGGCACGACCGGTAAGGGTGATACCCTGCTCAt                                                                                                         >  6:18825/1‑149 (MQ=255)                                                 gCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGcc                                                                                                 <  5:337150/148‑1 (MQ=255)                                                  cTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAg                                                                                               <  5:222497/149‑1 (MQ=255)                                                   tttGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAgg                                                                                              >  2:292282/1‑149 (MQ=255)                                                     tgtgGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAggt                                                                                             <  4:335586/148‑1 (MQ=255)                                                       tgGCGGTTGCATGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGGAACCGTCGAGACCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAggtcg                                                                                           <  1:96848/149‑1 (MQ=255)                                                          cGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGgtcggtcg                                                                                       >  2:149795/1‑149 (MQ=255)                                                          cGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGgtcggtcg                                                                                       <  1:162775/149‑1 (MQ=255)                                                                cAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCa                                                                                 <  4:152950/149‑1 (MQ=255)                                                                 aGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGgcgc                                                                                   >  2:191388/1‑146 (MQ=255)                                                                 aGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCa                                                                                 <  1:191388/148‑1 (MQ=255)                                                                 aGGGTAATGCTGAACTGGCCGTCAGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAg                                                                                >  6:293462/1‑149 (MQ=255)                                                                   ggTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGTGCCAgtg                                                                              <  1:49308/149‑1 (MQ=255)                                                                    gTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAgtgg                                                                             <  8:29458/149‑1 (MQ=255)                                                                      aaTGCTGAACTGGTCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAgtggtg                                                                           >  5:185268/1‑149 (MQ=255)                                                                      aaTGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAgtggtg                                                                           >  5:228485/1‑149 (MQ=255)                                                                       aTGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGcc                                                                                                 <  2:152885/126‑1 (MQ=255)                                                                       aTGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGcc                                                                                                 >  1:152885/1‑126 (MQ=255)                                                                         gCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAg                                                                        >  6:72914/1‑149 (MQ=255)                                                                          cTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTG‑‑CATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGc                                                                          >  2:90028/1‑144 (MQ=255)                                                                          cTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTG‑‑CATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGc                                                                          <  1:90028/144‑1 (MQ=255)                                                                           tGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGtt                                                                      >  1:155123/1‑149 (MQ=255)                                                                                  gCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGTAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTcc                                                                <  6:26980/148‑1 (MQ=255)                                                                                    cGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGt                                                             <  2:50848/149‑1 (MQ=255)                                                                                       ccGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCa                                                          <  2:344865/149‑1 (MQ=255)                                                                                        cGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAg                                                         >  3:304148/1‑149 (MQ=255)                                                                                        cGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAg                                                         <  3:256402/149‑1 (MQ=255)                                                                                          cACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGt                                                        >  4:52812/1‑148 (MQ=255)                                                                                           aCCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGgtgt                                                                                                                     <  2:8641/86‑1 (MQ=255)                                                                                           aCCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGgtgt                                                                                                                     >  1:8641/1‑86 (MQ=255)                                                                                           aCCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGa                                                      >  4:14435/1‑149 (MQ=255)                                                                                             cGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACa                                                    <  4:304148/149‑1 (MQ=255)                                                                                                    cGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGGCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAg                                             >  1:302261/1‑149 (MQ=255)                                                                                                        gTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGgcgc                                         >  1:301707/1‑149 (MQ=255)                                                                                                           cccAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGcc                                        >  6:17389/1‑147 (MQ=255)                                                                                                            ccAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCgag                                     >  8:102856/1‑149 (MQ=255)                                                                                                             cAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCgaga                                    <  2:301707/149‑1 (MQ=255)                                                                                                              aGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGAc                                   <  1:292282/149‑1 (MQ=255)                                                                                                               gCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACa                                  <  6:292066/149‑1 (MQ=255)                                                                                                                 aaTACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACAt                                 >  1:6903/1‑148 (MQ=255)                                                                                                                  aTACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACGGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCTGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTa                               <  5:293462/149‑1 (MQ=255)                                                                                                                      gCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACcgg                           >  1:73710/1‑149 (MQ=255)                                                                                                                        cACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACcggcg                         >  2:243738/1‑149 (MQ=255)                                                                                                                           cTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGgcggcg                      <  2:315065/149‑1 (MQ=255)                                                                                                                           cTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGgcggcg                      >  5:364779/1‑149 (MQ=255)                                                                                                                             gcgcCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCgt                     >  2:65552/1‑148 (MQ=255)                                                                                                                              cgcCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCgt                     <  2:302261/147‑1 (MQ=255)                                                                                                                               gcCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCgtcag                  <  2:155123/149‑1 (MQ=255)                                                                                                                                 cTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCgtcagt                 >  4:58334/1‑148 (MQ=255)                                                                                                                                    gCACGGAAACCGTCGAG‑CCGGGCTGGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAgg             <  5:18825/149‑1 (MQ=255)                                                                                                                                          aaaCCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACa       >  7:77673/1‑149 (MQ=255)                                                                                                                                              cGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCaa   <  4:58350/149‑1 (MQ=255)                                                                                                                                               gTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAg  <  7:102856/149‑1 (MQ=255)                                                                                                                                               gTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAg  >  6:22219/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GATATTGCCAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAG‑CCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAG  >  NC_002947/944680‑944970

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