Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,256,790	G→A	6.5%	intergenic (+37/+74)	mrp → / ← PP_mr14	ATP‑dependent Fe‑S cluster transferase/cspA regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,256,790	0	G	A	6.5%	178.0 / 2.1	61	intergenic (+37/+74)	mrp/PP_mr14	ATP‑dependent Fe‑S cluster transferase/cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base A (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAA  >  NC_002947/1256646‑1256934                                                                                                                                                 |                                                                                                                                                 cAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCg                                                                                                                                               >  6:54175/1‑148 (MQ=255)   gATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTAc                                                                                                                                            <  1:67317/149‑1 (MQ=255)        ccATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGaa                                                                                                                                       <  6:137991/149‑1 (MQ=255)         cATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAg                                                                                                                                      >  5:14985/1‑149 (MQ=255)         cATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAg                                                                                                                                      >  7:117363/1‑149 (MQ=255)         cATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAg                                                                                                                                      >  3:95486/1‑149 (MQ=255)                 aCCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCtgtg                                                                                                                                >  7:20212/1‑147 (MQ=255)                  ccAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAg                                                                                                                             <  5:54175/149‑1 (MQ=255)                   cAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGt                                                                                                                            >  1:313692/1‑149 (MQ=255)                     ggAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTcc                                                                                                                          <  6:173368/149‑1 (MQ=255)                        gCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCtg                                                                                                                       >  8:200229/1‑149 (MQ=255)                            gCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTgg                                                                                                                    <  3:209874/148‑1 (MQ=255)                              ccGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAg                                                                                                                 >  5:88450/1‑149 (MQ=255)                                gCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCg                                                                                                               >  3:211375/1‑149 (MQ=255)                                     gTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGc                                                                                                            >  5:306250/1‑147 (MQ=255)                                     gTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGc                                                                                                            >  2:9260/1‑147 (MQ=255)                                     gTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGa                                                                                                          <  2:202054/149‑1 (MQ=255)                                         gcgcgcGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGccc                                                                                                      >  8:41966/1‑149 (MQ=255)                                                 aTCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAagag                                                                                              <  8:117363/149‑1 (MQ=255)                                                 aTCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAagag                                                                                              <  2:152466/149‑1 (MQ=255)                                                  tCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAgg                                                                                             <  1:9260/149‑1 (MQ=255)                                                    gTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGc                                                                                            <  4:211375/148‑1 (MQ=255)                                                     tgctgcAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTAATTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGc                                                                                            >  3:167975/1‑147 (MQ=255)                                                       ctgcAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAgc                                                                                        <  7:120128/149‑1 (MQ=255)                                                            ggAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCTAAGACTGACGTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGGGAGCCCGCGAAGAGGCCAGCGCCTa                                                                                   <  4:95486/149‑1 (MQ=255)                                                               aGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTAttt                                                                                <  7:296092/149‑1 (MQ=255)                                                                 cTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTc                                                                               >  8:213609/1‑148 (MQ=255)                                                                    cggcgACGTCAATGCCCAGCATTACAATCAGCGAAGACTGACGTGCCATACGGAGTTTGCCCGCACAGAGCGACTGGGCCGTACTTGAAGGCTGTGCAGTCCCTGTGGTAGCTGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGa                                                                             <  8:112394/147‑1 (MQ=38)                                                                          cGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGc                                                                                                                >  5:93611/1‑106 (MQ=255)                                                                          cGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGc                                                                                                                <  6:93611/106‑1 (MQ=255)                                                                          cGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAt                                                                     <  3:157730/149‑1 (MQ=255)                                                                              aaTGCCCAGCATTACCATCAGCGAAGACTGCCCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATaaa                                                                  >  8:62998/1‑148 (MQ=255)                                                                                            ccATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAg                                                   >  3:247799/1‑149 (MQ=255)                                                                                             cATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAgg                                                  <  7:211196/149‑1 (MQ=255)                                                                                                 aGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                               <  8:98694/148‑1 (MQ=255)                                                                                                  gCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAg                                                                         <  8:246821/121‑1 (MQ=255)                                                                                                  gCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAg                                                                         >  7:246821/1‑121 (MQ=255)                                                                                                    gAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGttt                                           <  7:200229/149‑1 (MQ=255)                                                                                                     aaGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGtttt                                          >  5:269504/1‑149 (MQ=255)                                                                                                         cTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttttt                                       >  1:297888/1‑148 (MQ=255)                                                                                                            aCCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGttttt                                         >  3:63110/1‑143 (MQ=255)                                                                                                            aCCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGttttt                                         <  4:63110/143‑1 (MQ=255)                                                                                                            aCCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCACCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttt                                         <  3:73651/143‑1 (MQ=255)                                                                                                            aCCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCACCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttt                                         <  3:77937/143‑1 (MQ=255)                                                                                                            aCCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCACCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttt                                         >  4:77937/1‑143 (MQ=255)                                                                                                            aCCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCACCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttt                                         >  4:73651/1‑143 (MQ=255)                                                                                                             ccTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGc                                                 >  1:124231/1‑134 (MQ=255)                                                                                                                gCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTa                               >  8:195328/1‑149 (MQ=255)                                                                                                                 ccATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttttt                                       >  5:192214/1‑140 (MQ=255)                                                                                                                 ccATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttttt                                       <  6:192214/140‑1 (MQ=255)                                                                                                                 ccATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAg                              >  3:182074/1‑149 (MQ=255)                                                                                                                      ccGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                         <  7:213609/149‑1 (MQ=255)                                                                                                                       cGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                         >  7:194556/1‑148 (MQ=255)                                                                                                                             tGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATc                   <  7:62998/148‑1 (MQ=255)                                                                                                                                ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGCCt               <  6:306250/149‑1 (MQ=255)                                                                                                                                ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACt               <  7:195328/149‑1 (MQ=255)                                                                                                                                 cGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACtt              <  4:182074/149‑1 (MQ=255)                                                                                                                                  gCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTa             <  6:269504/149‑1 (MQ=255)                                                                                                                                    acacaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAgg           <  4:167975/149‑1 (MQ=255)                                                                                                                                        aGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGAt                    >  4:166853/1‑136 (MQ=255)                                                                                                                                        aGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGAt                    <  3:166853/136‑1 (MQ=255)                                                                                                                                                 gcCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGaa  <  6:204715/145‑1 (MQ=255)                                                                                                                                                 gcCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGaa  >  5:204715/1‑145 (MQ=255)                                                                                                                                                 |                                                                                                                                                 CAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAA  >  NC_002947/1256646‑1256934

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