Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,979,658	G→A	7.5%	A110T (GCC→ACC)	rpsA →	30S ribosomal protein S1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,979,658	0	G	A	7.5%	193.9 / 3.2	66	A110T (GCC→ACC)	rpsA	30S ribosomal protein S1
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base A (2/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGGGCGTCATCCCGCTCGAGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGTGCGTTCCTGCCGGGCTCCCTGGTTGATGTCCGCC  >  NC_002947/1979516‑1979802                                                                                                                                               |                                                                                                                                                 cGAGGGCGTCATCCCGCTCGAGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGt                                                                                                                                            >  1:53801/1‑149 (MQ=255) cGAGGGCGTCATCCCGCTCGAGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGt                                                                                                                                            >  6:134818/1‑149 (MQ=255)    gggCGTCATCCCGCTCGAGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCt                                                                                                                                         >  4:283231/1‑149 (MQ=255)    gggCGTCATCCCGCTCGAGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCt                                                                                                                                         <  3:308983/149‑1 (MQ=255)             cccGCTCGAGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCACCGAGTGCTGGATTGTTc                                                                                                                                >  6:18410/1‑149 (MQ=255)             cccGCTCGAGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCACCGAGTGCTGGATTGTTc                                                                                                                                >  2:76266/1‑149 (MQ=255)              ccGCTCGAGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTTGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTc                                                                                                                                <  3:303157/148‑1 (MQ=255)                    gAGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAag                                                                                                                          <  8:231094/148‑1 (MQ=255)                    gAGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAag                                                                                                                          >  2:79202/1‑148 (MQ=255)                     aGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAagca                                                                                                                        <  2:340196/149‑1 (MQ=255)                          ttCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAagcag                                                                                                                       <  4:10205/145‑1 (MQ=255)                          ttCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAagcag                                                                                                                       >  3:10205/1‑145 (MQ=255)                             tACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAagcagc                                                                                                                      <  8:3141/143‑1 (MQ=255)                             tACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAagcagc                                                                                                                      >  7:3141/1‑143 (MQ=255)                              aCAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTcgcc                                                                                                               >  6:283212/1‑149 (MQ=255)                                aaCGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTcgccgc                                                                                                             >  5:369659/1‑149 (MQ=255)                                 aCGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCgccgcc                                                                                                            >  2:29580/1‑149 (MQ=255)                                   gAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGccgccg                                                                                                           >  7:305180/1‑148 (MQ=255)                                     aGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCACCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCgaag                                                                                                        <  5:18410/149‑1 (MQ=255)                                     aGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCACCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCgaag                                                                                                        <  1:76266/149‑1 (MQ=255)                                        tGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGaagaag                                                                                                     >  3:100613/1‑149 (MQ=255)                                            gAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCt                                                                                                                               <  3:343920/119‑1 (MQ=255)                                            gAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCt                                                                                                                               >  4:343920/1‑119 (MQ=255)                                               cTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCaa                                                                                              <  5:134818/149‑1 (MQ=255)                                                tGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAg                                                                                             >  2:351803/1‑149 (MQ=255)                                                  aCCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAgg                                                                                            <  6:369659/148‑1 (MQ=255)                                                     aTCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGt                                                                                        >  7:70279/1‑149 (MQ=255)                                                      tCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAgg                                                                                            >  6:411462/1‑144 (MQ=255)                                                      tCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAgg                                                                                            <  5:411462/144‑1 (MQ=255)                                                      tCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGtt                                                                                       <  1:289420/149‑1 (MQ=255)                                                         aGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAgg                                                                                            <  6:179334/141‑1 (MQ=255)                                                         aGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAgg                                                                                            >  5:179334/1‑141 (MQ=255)                                                              ggTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGTCGAAGAAGTGGTCAAGGGCGTTATCAACg                                                                                <  1:29580/148‑1 (MQ=255)                                                              ggTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACgg                                                                               >  1:311041/1‑149 (MQ=255)                                                                             gTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCgaag                                                                                                        <  8:159430/109‑1 (MQ=255)                                                                             gTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCgaag                                                                                                        >  7:159430/1‑109 (MQ=255)                                                                                 cgcTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCtt                                                             <  4:100613/148‑1 (MQ=255)                                                                                  gcTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCa                                                           <  3:283231/149‑1 (MQ=255)                                                                                      gACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACt                                                         <  1:79202/147‑1 (MQ=255)                                                                                      gACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGt                                                       >  3:152212/1‑149 (MQ=255)                                                                                       aCGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTc                                                      <  8:70279/149‑1 (MQ=255)                                                                                       aCGCGGTCGAAGAAGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGAACCGAGTGATGGATTGTTCTGTAAGCAGCTTTCGGCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTc                                                      <  2:122489/149‑1 (MQ=255)                                                                                         gcgGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCACCGCCgaag                                                                                                        <  7:204145/97‑1 (MQ=255)                                                                                         gcgGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCACCGCCgaag                                                                                                        >  8:204145/1‑97 (MQ=255)                                                                                          cgGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGAc                                                   <  8:305180/149‑1 (MQ=255)                                                                                             tCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGtt                                                <  2:53801/149‑1 (MQ=255)                                                                                               gAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTaa                                              >  6:230321/1‑149 (MQ=255)                                                                                               gAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTaa                                              >  8:322796/1‑149 (MQ=255)                                                                                                   aCGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCTGCTTCa                                                           <  3:152544/132‑1 (MQ=255)                                                                                                   aCGGCTTTGGCGAAACCAACCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCa                                                           >  4:152544/1‑132 (MQ=255)                                                                                                                  ccAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGTGCGTTCCTg                           >  2:54207/1‑149 (MQ=255)                                                                                                                   cAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCCTCCGTGCGTTCCTGc                          >  7:323010/1‑149 (MQ=255)                                                                                                                    aaGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGTGCGTTCCTGcc                         >  5:379873/1‑149 (MQ=255)                                                                                                                     aGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGTGCGTTCCTGCCg                        <  6:100090/149‑1 (MQ=255)                                                                                                                      gCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAg                                                                     <  5:105795/103‑1 (MQ=255)                                                                                                                      gCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAg                                                                     >  6:105795/1‑103 (MQ=255)                                                                                                                      gCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGTGCGTTCCTGCCgg                       <  7:322796/149‑1 (MQ=255)                                                                                                                      gCTGCCCCGTGAAAAAGCCACGCGTGCAGAGTTCGTGATTGTTCTGGAAGTAGCTTTCTCCGCCGAAGAAGTGGTAAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGTGCTTTCCTGCCgg                       <  4:400453/149‑1 (MQ=255)                                                                                                                            ccGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGt                                    >  8:158505/1‑130 (MQ=255)                                                                                                                            ccGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGt                                    <  7:158505/130‑1 (MQ=255)                                                                                                                             cGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGTGCGTTCCTGCCGGGCTCCCt                >  6:368779/1‑149 (MQ=255)                                                                                                                                 aaaaaGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGTAGTAGTCGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGTGCGTTCCTGCCGGGCTCCCTGGt             <  5:283212/148‑1 (MQ=255)                                                                                                                                   aaaGCCAAGCGAGCAGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTAAAGGGAGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATACGTGCGTTCCTGCATGGCTCCCTGGTTGa          >  2:243565/1‑149 (MQ=255)                                                                                                                                      gCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGTCGCCGAAGAAGTGGTCAAGGGCGTTATCaa                                                                                  <  8:240808/74‑1 (MQ=255)                                                                                                                                      gCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGTCGCCGAAGAAGTGGTCAAGGGCGTTATCaa                                                                                  >  7:240808/1‑74 (MQ=255)                                                                                                                                        cAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGTGCGTTCCTGCCGGGCTCCCTGGTTGATGTcc     <  5:230321/149‑1 (MQ=255)                                                                                                                                           gcgcgcCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGTGCGTTCCTGCCGGGCTCCCTGGTTGATGTCCGcc  >  2:411801/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                 CGAGGGCGTCATCCCGCTCGAGCAGTTCTACAACGAAGCTGGCGAGCTGACCATCAAGGTCGGTGACGAAGTTCACGTTGCGCTGGACGCGGTCGAAGACGGCTTTGGCGAAACCAAGCTGTCCCGTGAAAAAGCCAAGCGCGCCGAGTGCTGGATTGTTCTGGAAGCAGCTTTCGCCGCCGAAGAAGTGGTCAAGGGCGTTATCAACGGTAAGGTTAAGGGCGGCTTCACTGTCGACGTTAACGGCATCCGTGCGTTCCTGCCGGGCTCCCTGGTTGATGTCCGCC  >  NC_002947/1979516‑1979802

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