Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	554,157	Δ1 bp	11.6%	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	.	11.6%	207.4 / 21.9	69	coding (97/714 nt)	rpsC	30S ribosomal protein S3
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base . (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 = 2.54e-01

GCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCC  >  NC_002947/554013‑554292                                                                                                                                                  |                                                                                                                                        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                                                                                                                                        >  7:97700/1‑132 (MQ=255)               aTCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGc                                                                                                                                        <  8:97700/132‑1 (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                                                                                                         >  7:104936/1‑147 (MQ=255)                               gCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTg                                                                                                         >  6:15443/1‑147 (MQ=255)                                ccATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGaa                                                                                                       <  5:170800/148‑1 (MQ=255)                                 cATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACgg                                                                                                                                  >  1:74827/1‑120 (MQ=255)                                 cATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTGCGTAC‑CAGACgg                                                                                                                                  <  2:74827/120‑1 (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                                                                                                      <  8:234638/140‑1 (MQ=255)                                         tGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAg                                                                                                      >  7:234638/1‑140 (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                                                                                        >  1:253667/1‑142 (MQ=255)                                                     ggACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAAcgcg                                                                                        <  2:253667/142‑1 (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)                                                                                                                                                  |                                                                                                                                        GCCGTGCTGATCGCATCGTCAAGCGGTCTTGCCATATCACTGTCAAGGTTGCGGACAAGTAACGGAGTCGATCAGATGGGTCAGAAAGTACATCCCACTGGCATTCGCCTGGGAATCGTCAAGGAGCACACCTCCGTCTG‑GTACGCAGACGGTGCTACTTACGCAGATTACCTGTTGAAGGATCTGAAAACGCGTGAGTACCTCCAAGACAAACTAAAAAGCGCGTCCGTAAGCCGTATCGATATTCATCGTCCGGCTCAAACTGCACGCATCACCATCC  >  NC_002947/554013‑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: ‑