Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,256,882	G→T	7.4%	noncoding (433/451 nt)	PP_mr14 ←	cspA regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,256,882	0	G	T	7.4%	210.4 / 3.3	66	noncoding (433/451 nt)	PP_mr14	cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base G (34/27);  new base T (3/2);  total (37/29)
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

CATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCT  >  NC_002947/1256738‑1257029                                                                                                                                                 |                                                                                                                                                     cATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAgg                                                                                                                                                  >  7:108117/1‑149 (MQ=255)         aaGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGtttt                                                                                                                                          >  5:200867/1‑149 (MQ=255)                     ccATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAg                                                                                                                              >  3:144774/1‑149 (MQ=255)                        tACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAgg                                                                                                                             >  5:337818/1‑147 (MQ=255)                         aCCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCGTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAg                                                                                                                          >  1:283779/1‑149 (MQ=255)                           cGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCTGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCt                                                                                                                        <  4:308993/149‑1 (MQ=255)                                   cccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGAc                                                                                                                <  6:105626/149‑1 (MQ=255)                                   cccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGAc                                                                                                                <  6:100373/149‑1 (MQ=255)                                          acaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCGCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGAAAGACTTAGGc                                                                                                          >  6:350870/1‑148 (MQ=255)                                          acaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGCCTAAgtc                                                                                                          >  2:284249/1‑146 (MQ=255)                                          acaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGc                                                                                                          >  2:14132/1‑148 (MQ=255)                                                cctatctcGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTg                                                                                                    <  5:166591/141‑1 (MQ=255)                                                 cTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGCCTAAGTCGATCtagc                                                                                                  <  1:284249/149‑4 (MQ=255)                                                    cgcCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAAc                                                                                                 <  6:200867/147‑1 (MQ=255)                                                     accGTACTTGAAGCCTGTGCAGTCCCTCTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAACGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCCATCTGAActac                                                                                              >  6:70870/2‑147 (MQ=255)                                                        gTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTg                                                                                                    >  6:166591/1‑140 (MQ=255)                                                         tACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttttt                                                                                                                                       >  1:133692/1‑104 (MQ=255)                                                         tACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttttt                                                                                                                                       <  2:133692/104‑1 (MQ=255)                                                              gAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTg                                                                                     >  2:49762/1‑149 (MQ=255)                                                               aaGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGttttt                                                                                                                                         <  6:243433/96‑1 (MQ=25)                                                               aaGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGttttt                                                                                                                                         >  5:243433/1‑96 (MQ=25)                                                               aaGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGc                                                                                    <  1:190522/149‑1 (MQ=255)                                                                   cTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTTAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGc                                                                                >  1:164926/1‑149 (MQ=255)                                                                   cTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGc                                                                                <  3:142141/149‑1 (MQ=255)                                                                   cTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGc                                                                                <  3:332552/149‑1 (MQ=255)                                                                         aGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCt                                                                          >  5:153291/1‑149 (MQ=255)                                                                          gTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTg                                                                         <  8:108117/149‑1 (MQ=255)                                                                           tCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTgg                                                                        >  5:293465/1‑149 (MQ=255)                                                                              cTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCt                                                                          <  8:47374/144‑1 (MQ=255)                                                                              cTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCt                                                                          >  7:47374/1‑144 (MQ=255)                                                                                 tgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCgg                                                                  >  7:61160/1‑149 (MQ=255)                                                                                    gAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTa                                                               <  1:190476/149‑1 (MQ=255)                                                                                         ggCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGcgac                                                          <  4:144774/149‑1 (MQ=255)                                                                                              gCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAg                                                     >  3:96808/1‑149 (MQ=255)                                                                                                 cgcgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGCCCGCGGGTAGCGACGAAAGTaa                                                  >  8:147504/1‑149 (MQ=255)                                                                                                  gcgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAAc                                                 <  1:14132/149‑1 (MQ=255)                                                                                                         ggCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTgg                                          <  6:337818/149‑1 (MQ=255)                                                                                                          gCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑c                                        <  2:164926/149‑1 (MQ=255)                                                                                                           ccAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGaa                                                                                                  >  4:8476/1‑91 (MQ=255)                                                                                                           ccAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGaa                                                                                                  <  3:8476/91‑1 (MQ=255)                                                                                                             aGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGCGTTTTTTTGAAAGCTAGGAAGCTAGAGCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTATCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTCG‑Cctt                                     >  4:306980/1‑149 (MQ=255)                                                                                                                 ccTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCttt                                 <  2:283779/149‑1 (MQ=255)                                                                                                                   tATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTaa                                                  >  3:322121/1‑131 (MQ=255)                                                                                                                   tATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTaa                                                  <  4:322121/131‑1 (MQ=255)                                                                                                                   tATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTc                                >  3:294637/1‑148 (MQ=255)                                                                                                                     tttCGAGCAGGCATAAAAAAACCCGCCAGAATGCCGGTTTTTTTGAAAGCTATGAAGCTAGATCGACTTTGGCGATCTTAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGTGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCa                               <  5:350870/147‑1 (MQ=255)                                                                                                                     tttCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTc                                <  4:294637/146‑1 (MQ=255)                                                                                                                           gCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGtttt                        >  7:240312/1‑148 (MQ=255)                                                                                                                             aGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTTCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGaa                     <  7:147504/149‑1 (MQ=255)                                                                                                                              ggCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAg                    >  6:15284/1‑149 (MQ=255)                                                                                                                               gCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTgg                                                                        >  6:309/1‑97 (MQ=255)                                                                                                                                 tcaaaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAGCCg                 >  2:117850/3‑149 (MQ=255)                                                                                                                                 aTAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAGCCg                 >  4:83906/1‑149 (MQ=255)                                                                                                                                    aaataaCCCGCCAGTAGGCGTGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCTCGGGTAGCGACGATAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAGCCGTc               <  5:70870/148‑1 (MQ=255)                                                                                                                                    aaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAg                    >  6:134354/1‑143 (MQ=255)                                                                                                                                    aaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAg                    <  5:134354/143‑1 (MQ=255)                                                                                                                                    aaaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTc                                >  6:352342/1‑131 (MQ=255)                                                                                                                                    aaaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTc                                <  5:352342/131‑1 (MQ=255)                                                                                                                                    aaaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAGCCGt                >  1:265192/1‑147 (MQ=255)                                                                                                                                     aaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAGCCGTCAg             >  3:192997/1‑149 (MQ=255)                                                                                                                                        aaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAGCCGTCAGCtt          >  8:355311/1‑149 (MQ=255)                                                                                                                                         aCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTg         >  3:243619/1‑149 (MQ=255)                                                                                                                                            cGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGa       <  1:49762/149‑1 (MQ=255)                                                                                                                                            cGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGAt      <  4:96808/149‑1 (MQ=255)                                                                                                                                              ccAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATgg    >  5:118585/1‑149 (MQ=255)                                                                                                                                                aGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑ACTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCt  >  5:185754/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGG‑CCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCT  >  NC_002947/1256738‑1257029

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