Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,337,784	T→A	8.6%	S330T (TCG→ACG) ‡	PP_2934 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,337,784	0	T	A	8.6%	186.3 / 6.9	70	S330T (TCG→ACG) ‡	PP_2934	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base A (2/4);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

GGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCC  >  NC_002947/3337637‑3337930                                                                                                                                                    |                                                                                                                                                     gggCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTc                                                                                                                                                     >  5:220535/1‑149 (MQ=255) gggCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTc                                                                                                                                                     <  2:285349/149‑1 (MQ=255)          gTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAgc                                                                                                                                               >  4:167344/1‑146 (MQ=255)          gTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAgc                                                                                                                                               <  3:167344/146‑1 (MQ=255)             ttCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑gg                                                                                                                                         <  3:212478/148‑1 (MQ=255)             ttCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGc                                                                                                                                        >  7:336933/1‑149 (MQ=255)             ttCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGc                                                                                                                                        >  4:329836/1‑149 (MQ=255)               cATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCCGGCGCTGGCGACCCCAGCGcctgg                                                                                                                                         <  4:76208/147‑4 (MQ=255)               cATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCCGGCGCTGGCGACCCCAGCGcctgg                                                                                                                                         >  3:76208/1‑144 (MQ=255)                   gCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCg                                                                                                                                  <  2:355851/149‑1 (MQ=255)                    cccACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCCGGCGCTGGCGACCCCAGCGCCTGGCTGaaga                                                                                                                                  <  1:157835/149‑4 (MQ=255)                     ccACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGTGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGtt                                                                                                                                >  8:150792/1‑149 (MQ=255)                     ccACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGtt                                                                                                                                >  4:104688/1‑149 (MQ=255)                       aCTGCGATGCCGCGAATTAGGCGCGAAGGCGGGGGCGTTGCAAGGACGCCTCCGGCGGCGACCGCCGGTATATCCATCATTCGCTGCAGCGATAAACGGCAGCGTTCTTCAGCCAGGCGCTGGCTACCCCAGCGcctgg                                                                                                                                         <  2:191488/139‑4 (MQ=255)                       aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGCGACCCCAGCGcctgg                                                                                                                                         >  1:191488/1‑136 (MQ=255)                            gATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGc                                                                                                                         >  8:10403/1‑149 (MQ=255)                              tGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGa                                                                                                                       >  8:353877/1‑149 (MQ=255)                                 cGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACgg                                                                                                                    >  8:178167/1‑149 (MQ=255)                                  gAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACggg                                                                                                                   >  8:87847/1‑149 (MQ=255)                                      aGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGgcgcg                                                                                                               >  8:171029/1‑149 (MQ=255)                                       gTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGgcgcgc                                                                                                              >  5:31215/1‑149 (MQ=255)                                        tACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCACCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGgcgcgc                                                                                                              >  1:136136/1‑148 (MQ=255)                                           gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCTTCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGGTCAGCTGCGACGGGcgcgccgcg                                                                                                          >  1:104700/1‑149 (MQ=255)                                           gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGtt                                                                                                                                >  1:146972/1‑127 (MQ=255)                                           gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGtt                                                                                                                                <  2:146972/127‑1 (MQ=255)                                              cgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgcca                                                                                                       <  7:10403/149‑1 (MQ=255)                                                aaTGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagc                                                                                                     <  8:336933/149‑1 (MQ=255)                                                     cGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGTAGCGTTCTTCAGCCATGCGTTGGCGACTCTAGCGcctgg                                                                                                                                         <  8:56325/109‑4 (MQ=255)                                                        cgcTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGTGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGt                                                                                             <  7:150792/149‑1 (MQ=255)                                                             tgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAg                                                                                          <  4:34171/147‑1 (MQ=255)                                                               cAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGcc                                                                                       <  7:325871/148‑1 (MQ=255)                                                               cAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCg                                                                                      <  6:220535/149‑1 (MQ=255)                                                                aaGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACATTAGGCCGa                                                                                     >  2:335812/1‑149 (MQ=255)                                                                    aCGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGAcccc                                                                                 >  6:106408/1‑149 (MQ=255)                                                                      gCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagc                                                                                                     <  7:38470/127‑1 (MQ=255)                                                                      gCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagc                                                                                                     >  8:38470/1‑127 (MQ=255)                                                                      gCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCgg                                                                               >  5:162278/1‑149 (MQ=255)                                                                        cGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGc                                                                              >  5:328796/1‑148 (MQ=255)                                                                            ggCGTCGACCGCCGGTATATCCATCATTGGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAgg                                                                         <  5:303154/149‑1 (MQ=255)                                                                            ggCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCACGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGCACCCGGCCAAgg                                                                         >  4:237197/1‑149 (MQ=255)                                                                             gCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCTAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCCGGGGC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCGGGCCAAGGc                                                                        <  6:187713/149‑1 (MQ=255)                                                                               gTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGa                                                                      <  2:104700/149‑1 (MQ=255)                                                                                  gACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCa                                                                   <  7:87847/149‑1 (MQ=255)                                                                                       ccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCa                                                             <  3:294926/149‑1 (MQ=255)                                                                                               tctatcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGc                                                      <  4:123971/145‑1 (MQ=255)                                                                                                ccATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCCGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTCAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACg                                                       >  2:241230/1‑146 (MQ=255)                                                                                                ccATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACg                                                       >  6:275609/1‑146 (MQ=255)                                                                                                 cctctttCGCTGCAGCGAGATACGGCAGCGTTCTTAAGCCAGGCGCTGGGGTCGCCAGCGCC‑TGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCag                                                    <  3:233784/143‑1 (MQ=255)                                                                                                 catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCaga                                                   <  7:171029/149‑1 (MQ=255)                                                                                                 catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCag                                                    >  4:45984/1‑148 (MQ=255)                                                                                                 catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCag                                                    >  3:171321/1‑148 (MQ=255)                                                                                                  atcatTCGCTGCAGCGAGAAAGGGCAGGGTTCTTCAGCCAGGCGGTGGGGTCACCAGCGCC‑GGGTGAGGGTTCAGCTGGGACGGGCGCGCGGCGCCAGCAACAACGTTAGGCCTACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCagaa                                                  <  2:136136/149‑1 (MQ=255)                                                                                                      ttCGCTGCAGCGAGGAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCCGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGAGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGAt                                               >  7:320270/1‑148 (MQ=255)                                                                                                         gCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGa                                            <  3:329836/148‑1 (MQ=255)                                                                                                         gCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAg                                           <  3:45984/149‑1 (MQ=255)                                                                                                                 gagaAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATAccc                                   >  6:311060/1‑149 (MQ=255)                                                                                                                      aCGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGgt                              <  7:83649/149‑1 (MQ=255)                                                                                                                           aGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAAATCGAGGCATACCCCAGGTGTAcc                         <  3:212605/149‑1 (MQ=255)                                                                                                                                  ttcAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACgg                  >  4:149601/1‑149 (MQ=255)                                                                                                                                   tcAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAcc                         <  3:110066/141‑1 (MQ=255)                                                                                                                                   tcAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAcc                         >  4:110066/1‑141 (MQ=255)                                                                                                                                    cAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACgg                  <  1:89111/148‑1 (MQ=255)                                                                                                                                    cAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGTAGATCGAGGCATACCCCAGGTGTACCGCCACGgc                 <  6:328796/148‑1 (MQ=255)                                                                                                                                      gCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGcc              >  5:134389/1‑149 (MQ=255)                                                                                                                                      gCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCCTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGCTGTACCGCCACGGCGcc              >  1:252073/1‑149 (MQ=255)                                                                                                                                       ccAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAg                                           <  2:286579/119‑1 (MQ=255)                                                                                                                                        cAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa            >  4:127565/1‑149 (MQ=255)                                                                                                                                          ggCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCAt           >  4:110798/1‑148 (MQ=255)                                                                                                                                                 gggTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGc   <  6:31215/149‑1 (MQ=255)                                                                                                                                                  ggTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGcc  <  5:106408/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCC  >  NC_002947/3337637‑3337930

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