Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,337,782	G→C	8.7%	G329A (GGG→GCG)	PP_2934 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,337,782	0	G	C	8.7%	194.5 / 7.1	68	G329A (GGG→GCG)	PP_2934	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base C (2/4);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

CCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCC  >  NC_002947/3337635‑3337930                                                                                                                                                    |                                                                                                                                                       ccGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTgggg                                                                                                                                                       >  8:83649/1‑149 (MQ=255)   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                                                                                                                                        >  4:329836/1‑149 (MQ=255)               ttCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGc                                                                                                                                        >  7:336933/1‑149 (MQ=255)                 cATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCCGGCGCTGGCGACCCCAGCGcctgg                                                                                                                                         >  3:76208/1‑144 (MQ=255)                 cATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCCGGCGCTGGCGACCCCAGCGcctgg                                                                                                                                         <  4:76208/147‑4 (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                                                    >  3:171321/1‑148 (MQ=255)                                                                                                   catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCag                                                    >  4:45984/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                         >  4:110066/1‑141 (MQ=255)                                                                                                                                     tcAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAcc                         <  3:110066/141‑1 (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)                                                                                                                                                    |                                                                                                                                                       CCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCA‑GGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCC  >  NC_002947/3337635‑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: ‑