Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,157	Δ1 bp	6.5%	coding (97/714 nt)	rpsC →	30S ribosomal protein S3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	554,157	0	G	.	6.5%	299.0 / 12.8	92	coding (97/714 nt)	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand):  ref base G (45/41);  new base . (4/2);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCGCTC  >  NC_002947/554009‑554301                                                                                                                                                      |                                                                                                                                                    aaaGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACg                                                                                                                                                     >  8:38671/1‑149 (MQ=255) aaaGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACg                                                                                                                                                     >  5:235887/1‑149 (MQ=255) aaaGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACg                                                                                                                                                     >  2:275840/1‑149 (MQ=255)  aaGGCCGTGCTAATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGc                                                                                                                                                    <  6:226472/149‑1 (MQ=255)         tGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACg                                                                                                                                               >  6:239072/1‑147 (MQ=255)           cTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACggag                                                                                                                                            <  5:239072/148‑3 (MQ=255)            tGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGt                                                                                                                                             <  5:342712/146‑1 (MQ=255)            tGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGt                                                                                                                                             >  6:342712/1‑146 (MQ=255)             gATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTa                                                                                                                                         >  3:290593/1‑149 (MQ=255)              aTCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCt                                                                                                                                          <  3:16945/147‑1 (MQ=255)                  cATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTAc                                                                                                                                    >  3:366194/1‑149 (MQ=255)                   aTCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACg                                                                                                                                   >  3:348158/1‑149 (MQ=255)                   aTCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGAACAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACg                                                                                                                                   >  4:310136/1‑149 (MQ=255)                         aaGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGAtt                                                                                                                             <  7:306296/149‑1 (MQ=255)                         aaGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGAt                                                                                                                              >  2:192108/1‑148 (MQ=255)                         aaGCGGTCTTGCCATATCACTGACAACCTTGCGGACAAGTAACGCAGTCTATCAGATGGCTAAGAAAGTACATCCCACTAGCATTCGCATGCCAATCGTCAAGGAGCAAACCTCCGTCTG‑GTACGAAGACGGAGCTACTTAAGCAGaa                                                                                                                              >  5:234301/1‑147 (MQ=255)                           gCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTAc                                                                                                                           <  8:331661/149‑1 (MQ=255)                               tCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGt                                                                                                                       <  3:191542/149‑1 (MQ=255)                                 ttGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGtt                                                                                                                      >  6:297943/1‑148 (MQ=255)                                   gCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGtt                                                                                                                      <  4:198740/146‑1 (MQ=255)                                   gCCAAATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTAGGCAGATTACCTGTTg                                                                                                                     <  1:359140/147‑1 (MQ=255)                                     cATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAg                                                                                                                  >  7:284265/1‑148 (MQ=255)                                      atatCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGa                                                                                                                >  1:161753/1‑149 (MQ=255)                                      atatCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGa                                                                                                                >  3:220013/1‑149 (MQ=255)                                      atatCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCAG‑GTACGCAGACGGTGCAACTTACGCAGATTACCTGTTGAAGGa                                                                                                                >  2:200760/1‑149 (MQ=255)                                       tatCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGAt                                                                                                               <  8:284265/149‑1 (MQ=255)                                       tatCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTAACGCAGATTACCTGTTGAAGGAt                                                                                                               <  8:122569/149‑1 (MQ=255)                                         tCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCt                                                                                                             <  7:235171/149‑1 (MQ=255)                                          cACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCAGTCTG‑GTAAGAAGACGGTGCTACTTCCGCAGATGACCTGTTGAAGGAtag                                                                                                             >  4:51292/1‑146 (MQ=255)                                            cTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCACCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGc                                                                                                                                           <  2:191410/116‑1 (MQ=255)                                            cTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCACCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGc                                                                                                                                           >  1:191410/1‑116 (MQ=255)                                             tGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCa                                                                                                                                 <  2:379284/125‑1 (MQ=255)                                             tGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCa                                                                                                                                 >  1:379284/1‑125 (MQ=255)                                              ttCAAGGTTGCGGACAAGTAACGGAGTAGATCAGATGGGTCACAAAGTACATCCCACTGGCATTAGCCTGGGAATCGTCAAGGAGCACAACTCCGTCTC‑GTACGCAGACCGAGCTACTTACGCAGATTACCTCTTGAAGGATCTGaaaa                                                                                                        >  4:350982/2‑149 (MQ=255)                                               tCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAAc                                                                                                       >  2:85337/1‑149 (MQ=255)                                               tCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAAc                                                                                                       >  1:316366/1‑149 (MQ=255)                                               tCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACATCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAAc                                                                                                       >  4:33348/1‑149 (MQ=255)                                                  aGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCGGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAAcgcg                                                                                                    >  8:231902/1‑149 (MQ=255)                                                   ggTTGCGGATAAGTAACGGAGTCGATCAGTTGGGTCAGAAAGTACATCCCACTGTCATTCGCCTGGGATTCGTCAAGGAGCACAGCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCATATTACCTGTTGAAGGATCTGAAAACGCGt                                                                                                   <  3:315970/149‑1 (MQ=255)                                                   ggTTGCGGACAAGTAACGGAGTCGTTCAGATGGGTCAGATTGTACATCCGTCTGGCATTCGCCTGGGAATCGTCAAGGAGCAGACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGTTTAGGTGTTGAAGTATCTGAAAACGCGt                                                                                                   <  6:234301/149‑1 (MQ=255)                                                   ggTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGt                                                                                                   <  6:174246/149‑1 (MQ=255)                                                   ggTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGt                                                                                                   <  6:235887/149‑1 (MQ=255)                                                     ttGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGa                                                                                                 <  7:38671/149‑1 (MQ=255)                                                             aaGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTc                                                                                          >  8:133009/1‑148 (MQ=255)                                                             aaGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTc                                                                                          >  6:20283/1‑148 (MQ=255)                                                               gTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTcc                                                                                         >  6:347250/1‑147 (MQ=255)                                                                     gAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAACACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACaaa                                                                                 <  1:188823/149‑1 (MQ=255)                                                                      aGTCGAACAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAAc                                                                                <  3:310136/149‑1 (MQ=255)                                                                        tCGATAAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACt                                                                               >  2:80244/1‑148 (MQ=255)                                                                            tCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTaaaaa                                                                          <  4:290593/149‑1 (MQ=255)                                                                            tCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTaaaaa                                                                          <  8:152303/149‑1 (MQ=255)                                                                              aGATGGGTCAGAAAGTCCATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAg‑c                                                                       <  4:273271/149‑1 (MQ=255)                                                                              aGATGGGTCAGAAAGTACATCACACTTGCATTCGCCTGGGAATCGTCAAGGAGCACAACTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTTAAGTATCTGAAAACGCGTGAGAACCTCCAAGACAAACTAAAAAg‑c                                                                       <  4:46777/149‑1 (MQ=255)                                                                              aGATGGGGCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAg‑c                                                                       <  6:285190/149‑1 (MQ=255)                                                                               gATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAg‑cg                                                                      >  8:114014/1‑149 (MQ=255)                                                                               gATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAg‑cg                                                                      <  5:297943/149‑1 (MQ=255)                                                                               gATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAg‑cg                                                                      >  8:205463/1‑149 (MQ=255)                                                                                aTGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAg‑cgc                                                                     <  2:161753/149‑1 (MQ=255)                                                                                aTGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAg‑cgc                                                                     <  4:348158/149‑1 (MQ=255)                                                                                   ggTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTc                                                                  <  2:316366/149‑1 (MQ=255)                                                                                    gtgaGAAAGTACATCCCACTGGCATTCGCCTGGTAATCGTAAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTcc                                                                 <  3:350982/146‑1 (MQ=255)                                                                                     tCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCg                                                                >  5:6311/1‑149 (MQ=255)                                                                                          aaGTACATCCCACTGGCATTCGCCTGGGAATCGTCAATGAGCAGACCTGCGTCTT‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑Ag                                                           <  1:200760/148‑1 (MQ=255)                                                                                             tACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCg                                                        >  5:214542/1‑148 (MQ=255)                                                                                                  cccACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑Ga                                                 >  8:207546/1‑149 (MQ=255)                                                                                                  cccACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑Ga                                                 >  3:59250/1‑149 (MQ=255)                                                                                                         gCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCTTCCGTA‑AGCCGTATC‑GATATTCAt                                          >  1:200390/1‑149 (MQ=255)                                                                                                         gCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCa                                           <  4:59250/148‑1 (MQ=255)                                                                                                              cGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTc                                      <  5:20283/148‑1 (MQ=255)                                                                                                                 cTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGGTCTGAAAACGCGTGa                                                                                                 <  5:392729/89‑1 (MQ=255)                                                                                                                 cTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGGTCTGAAAACGCGTGa                                                                                                 >  6:392729/1‑89 (MQ=255)                                                                                                                      aaTCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCaaa                             <  4:366194/149‑1 (MQ=255)                                                                                                                        tCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACt                           >  3:157054/1‑149 (MQ=255)                                                                                                                           tCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATc                                         >  8:70617/1‑132 (MQ=255)                                                                                                                           tCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATc                                         <  7:70617/132‑1 (MQ=255)                                                                                                                            cAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCAc                       >  7:87911/1‑149 (MQ=255)                                                                                                                             aaGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCCGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCAc                       <  1:85337/149‑1 (MQ=255)                                                                                                                             aaGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCTTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCAc                       <  2:200390/148‑1 (MQ=255)                                                                                                                               ggAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCa                    >  2:103154/1‑149 (MQ=255)                                                                                                                                  gCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCa                 >  4:324199/1‑149 (MQ=255)                                                                                                                                   cacaCCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTc                                            >  4:236233/1‑121 (MQ=255)                                                                                                                                   cacaCCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTc                                            <  3:236233/121‑1 (MQ=255)                                                                                                                                   cacaCCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCAc                >  3:46423/1‑149 (MQ=255)                                                                                                                                     cacCTCCGTGTG‑GTACGCATACGGTGCTACTTACGCATAATACCTGTTGAATGATCTGAAAACGCGTGAGTACCTCCAACACAAACTAAAAAG‑CGCGTCCGTACAGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCAcc               <  1:275840/149‑1 (MQ=255)                                                                                                                                      acCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCAt             >  3:292074/1‑149 (MQ=255)                                                                                                                                      acCTCCGCCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCAt             >  7:953/1‑149 (MQ=255)                                                                                                                                      acCTCCGCCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCAt             >  7:6585/1‑149 (MQ=255)                                                                                                                                        caccGTCTGCGTAC‑CAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCAc                       >  5:48024/3‑137 (MQ=255)                                                                                                                                        caccGTCTGCGTAC‑CAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCAc                       <  6:48024/135‑1 (MQ=255)                                                                                                                                        cTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCATc            >  4:252220/1‑148 (MQ=255)                                                                                                                                          ccGTCTGCGTAC‑CAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCca          >  8:208818/1‑148 (MQ=255)                                                                                                                                          ccGTCTGCGTAC‑CAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCca          >  8:214375/1‑148 (MQ=255)                                                                                                                                             tCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATCGGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCacac       <  1:80244/149‑1 (MQ=255)                                                                                                                                                g‑gTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCGCt   >  5:273306/1‑149 (MQ=255)                                                                                                                                                  gTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCGCTc  >  2:33566/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCGCTC  >  NC_002947/554009‑554301

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