Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,152:1	+C	6.5%	coding (92/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,152	1	.	C	6.5%	366.8 / 19.9	93	coding (92/714 nt)	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand):  ref base . (46/41);  new base C (4/2);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

TGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCGCT  >  NC_002947/554005‑554300                                                                                                                                                     |                                                                                                                                                        tGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGCGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTg‑g                                                                                                                                                        <  6:112677/149‑1 (MQ=255)     aaaGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACg                                                                                                                                                    >  5:235887/1‑149 (MQ=255)     aaaGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACg                                                                                                                                                    >  8:38671/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                                                                                                               >  3:220013/1‑149 (MQ=255)                                          atatCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGa                                                                                                               >  1:161753/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                                                                                                                                          >  1:191410/1‑116 (MQ=255)                                                cTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCACCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGc                                                                                                                                          <  2:191410/116‑1 (MQ=255)                                                 tGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCa                                                                                                                                >  1:379284/1‑125 (MQ=255)                                                 tGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCa                                                                                                                                <  2:379284/125‑1 (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:235887/149‑1 (MQ=255)                                                       ggTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGt                                                                                                  <  6:174246/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                                                                     >  8:205463/1‑149 (MQ=255)                                                                                   gATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAg‑cg                                                                     <  5:297943/149‑1 (MQ=255)                                                                                    aTGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAg‑cgc                                                                    <  4:348158/149‑1 (MQ=255)                                                                                    aTGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAg‑cgc                                                                    <  2:161753/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                                        <  7:70617/132‑1 (MQ=255)                                                                                                                               tCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATc                                        >  8:70617/1‑132 (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                                           <  3:236233/121‑1 (MQ=255)                                                                                                                                       cacaCCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTc                                           >  4:236233/1‑121 (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                      <  6:48024/135‑1 (MQ=255)                                                                                                                                            caccGTCTGCGTAC‑CAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCAc                      >  5:48024/3‑137 (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)                                                                                                                                                     |                                                                                                                                                        TGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAG‑CGCGTCCGTA‑AGCCGTATC‑GATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCCACACCGCT  >  NC_002947/554005‑554300

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