Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,152:1	+C	11.6%	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	11.6%	233.6 / 29.0	69	coding (92/714 nt)	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand):  ref base . (28/33);  new base C (3/5);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-01

TGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCC  >  NC_002947/554005‑554292                                                                                                                                                     |                                                                                                                                             tGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGGCTGGGAATCGTCAAGGAGCACACCTCCGTCTg‑g                                                                                                                                             <  7:279033/149‑1 (MQ=255)    cAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCACCGTCTgcgta                                                                                                                                           >  6:114526/1‑145 (MQ=255)       aGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTgcgtacca                                                                                                                                        <  6:311514/149‑8 (MQ=255)         gCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCTTATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGa                                                                                                                                     >  4:43009/1‑149 (MQ=255)             tGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACgga                                                                                                                                 <  7:115860/149‑2 (MQ=255)               cTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCACCGTCTGCGTAC‑CAGACggagg                                                                                                                               <  5:114526/149‑4 (MQ=255)                  aTCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACggaggtg                                                                                                                             <  7:275271/148‑6 (MQ=255)                       aTCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGAGATGTTACGGAGTCGTTGTGCTGGGGCAGAAAGTACATCCCACTGGCATTCGCCTGTGATTCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTAc                                                                                                                        <  6:262153/148‑1 (MQ=255)                       aTCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGc                                                                                                                                        <  8:97700/132‑1 (MQ=255)                       aTCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGc                                                                                                                                        >  7:97700/1‑132 (MQ=255)                           tCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACg                                                                                                                       >  5:144836/1‑145 (MQ=255)                                cGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTAcc                                                                                                              <  6:91744/149‑1 (MQ=255)                                 ggTCTTGCCTTATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCt                                                                                                             <  3:43009/149‑1 (MQ=255)                                       gCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTg                                                                                                         >  6:15443/1‑147 (MQ=255)                                       gCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTg                                                                                                         >  7:104936/1‑147 (MQ=255)                                        ccATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGaa                                                                                                       <  5:170800/148‑1 (MQ=255)                                         cATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACgg                                                                                                                                  <  2:74827/120‑1 (MQ=255)                                         cATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACgg                                                                                                                                  >  1:74827/1‑120 (MQ=255)                                         cATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAg                                                                                                      >  3:257501/1‑148 (MQ=255)                                          atatCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACgg                                                                                                                                  >  7:340985/1‑119 (MQ=255)                                          atatCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACgg                                                                                                                                  <  8:340985/119‑1 (MQ=255)                                          atatCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAgg                                                                                                     >  4:175779/1‑148 (MQ=255)                                             tCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATc                                                                                                  <  2:222038/148‑1 (MQ=255)                                               aCTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTg                                                                                                <  5:288579/148‑1 (MQ=255)                                                 tGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAg                                                                                                      >  7:234638/1‑140 (MQ=255)                                                 tGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAg                                                                                                      <  8:234638/140‑1 (MQ=255)                                                  gTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTAcc                                                                                                              >  4:362171/1‑131 (MQ=255)                                                  gTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTAcc                                                                                                              <  3:362171/131‑1 (MQ=255)                                                   tCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAAc                                                                                           >  3:169559/1‑149 (MQ=255)                                                     aaGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAAcgc                                                                                         >  7:243485/1‑149 (MQ=255)                                                      tggTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGAGCTGAAAAcgcg                                                                                        <  7:347948/148‑1 (MQ=255)                                                       ggTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGt                                                                                       <  8:104936/149‑1 (MQ=255)                                                         ttGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGa                                                                                     <  2:10412/149‑1 (MQ=255)                                                             ggACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAAcgcg                                                                                        <  2:253667/142‑1 (MQ=255)                                                             ggACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAAcgcg                                                                                        >  1:253667/1‑142 (MQ=255)                                                                cAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTc                                                                              <  3:175779/149‑1 (MQ=255)                                                                   gTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCaa                                                                           <  2:99211/149‑1 (MQ=255)                                                                    tAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTTTTCAAGGATATGAAAACGGGTGAGTAcgaccat                                                                           >  7:251374/1‑142 (MQ=255)                                                                    tAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTcc                                                                             >  7:265645/1‑146 (MQ=255)                                                                         gAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACaaa                                                                     <  8:101155/149‑1 (MQ=255)                                                                              gATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGaa                                                                                                       >  8:330740/1‑110 (MQ=255)                                                                              gATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGaa                                                                                                       <  7:330740/110‑1 (MQ=255)                                                                              gATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTaaa                                                                <  1:158007/149‑1 (MQ=255)                                                                               aTCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTaaa                                                                >  8:103540/1‑148 (MQ=255)                                                                                tCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTaaaaa                                                              <  6:251758/149‑1 (MQ=255)                                                                                tCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTaaaaa                                                              <  7:103540/149‑1 (MQ=255)                                                                                tCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTaaa                                                                >  4:111972/1‑147 (MQ=255)                                                                                  aGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAgc                                                            <  4:169559/149‑1 (MQ=255)                                                                                    aTGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAgcgc                                                          <  4:287728/149‑1 (MQ=255)                                                                                        gTCAGAAAGTACATCCCACTGTCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTcc                                                      <  4:257501/149‑1 (MQ=255)                                                                                            gAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTaa                                                  <  8:243485/149‑1 (MQ=255)                                                                                                   cATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTAt                                           >  7:199930/1‑149 (MQ=255)                                                                                                              cATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGGAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCATCCGTAAGCCGTATCGATATTCaac                                <  8:265645/149‑3 (MQ=255)                                                                                                               aTTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAACGCGCGTCCGTAAGCCGTATCGATATTCATCg                               >  2:147554/1‑149 (MQ=255)                                                                                                                     cTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGc                         >  1:173030/1‑149 (MQ=255)                                                                                                                       gggAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCt                        >  6:194398/1‑148 (MQ=255)                                                                                                                        ggAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTc                       >  2:44544/1‑148 (MQ=255)                                                                                                                               tCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCATATTACCTGTTGAAGGATCTGAAACCGCGTGAGTACCTCCAAGACAAACTAAAAAgc                                                            <  8:302808/104‑1 (MQ=255)                                                                                                                               tCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAgc                                                            >  7:302808/1‑104 (MQ=255)                                                                                                                                 aaGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCAc              <  3:111972/148‑1 (MQ=255)                                                                                                                                    gAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCAt          >  1:357841/1‑149 (MQ=255)                                                                                                                                     aGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATc         <  2:357841/149‑1 (MQ=255)                                                                                                                                     aGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATc         <  6:144836/149‑1 (MQ=255)                                                                                                                                      gCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCa        <  8:199930/149‑1 (MQ=255)                                                                                                                                       cacaCCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGCCAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCAc       >  6:10320/1‑149 (MQ=255)                                                                                                                                          acCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCAt    >  5:238029/1‑149 (MQ=255)                                                                                                                                            caccGTCTGCGTAC‑CAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATc   >  8:119438/3‑148 (MQ=255)                                                                                                                                            cTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCaa                                                                           >  5:253327/1‑76 (MQ=255)                                                                                                                                            cTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCaa                                                                           <  6:253327/76‑1 (MQ=255)                                                                                                                                            cTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCAcc      >  5:93106/1‑145 (MQ=255)                                                                                                                                            cTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCAcc      <  6:93106/145‑1 (MQ=255)                                                                                                                                            cTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATcc  <  5:15443/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             TGCCAAAGGCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCC  >  NC_002947/554005‑554292

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