Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,125,370	G→C	5.6%	G55R (GGT→CGT)	PP_4512 →	DMT superfamily permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,125,370	0	G	C	5.6%	273.2 / 3.8	83	G55R (GGT→CGT)	PP_4512	DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base G (38/40);  new base C (2/3);  total (40/43)
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

CGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGC  >  NC_002947/5125222‑5125518                                                                                                                                                     |                                                                                                                                                      cgcCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCg                                                                                                                                                       <  3:267172/149‑1 (MQ=255)     cTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTcc                                                                                                                                                   >  6:298328/1‑149 (MQ=255)        ggCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCggg                                                                                                                                                <  6:163992/149‑1 (MQ=255)          cgcACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCt                                                                                                                                              >  3:204324/1‑149 (MQ=255)           gcACTGCACATTGGCGCGCTGTTGTTCGGACTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTg                                                                                                                                             <  8:160102/149‑1 (MQ=255)           gcACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTg                                                                                                                                             <  6:38983/149‑1 (MQ=255)           gcACTGCACATTGGCGCGCGGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGCCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCTCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTg                                                                                                                                             <  6:251205/149‑1 (MQ=255)            cACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTgg                                                                                                                                            <  5:214916/149‑1 (MQ=255)                gCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAAc                                                                                                                                        >  4:211718/1‑149 (MQ=255)                  acaTTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGc                                                                                                                                      <  4:173382/149‑1 (MQ=255)                    aTTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCt                                                                                                                                    >  4:315100/1‑149 (MQ=255)                    aTTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCACCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCt                                                                                                                                    <  6:240644/149‑1 (MQ=255)                       ggCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGa                                                                                                                                  <  1:327668/148‑1 (MQ=255)                         cgcgcTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGagca                                                                                                                               <  7:258239/149‑1 (MQ=255)                          gcgcTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTg                                                                                                                                             <  6:354671/134‑1 (MQ=255)                          gcgcTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTg                                                                                                                                             >  5:354671/1‑134 (MQ=255)                            gcTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGcc                                                                                                                            <  5:298328/149‑1 (MQ=255)                                ttcttcGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGc                                                                                                                                           >  8:88199/1‑130 (MQ=255)                                ttcttcGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGc                                                                                                                                           <  7:88199/130‑1 (MQ=255)                                   ttcGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTACGGTAGCGCCCACTTCGCCGTACTCGCGCTGGACCTGTTCGCCAGCCTGGCCCGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATgc                                                                                                                     >  8:169911/1‑149 (MQ=255)                                       gCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATgcgcgc                                                                                                                 <  4:204324/149‑1 (MQ=255)                                            aCCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGcgccgc                                                                                                              <  7:152056/147‑1 (MQ=255)                                            aCCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCACGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCCGAGCAG‑CAGGCAGCGCGCATCct‑g                                                                                                          >  3:227067/1‑149 (MQ=255)                                              cGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCct‑gc                                                                                                         >  8:124007/1‑149 (MQ=255)                                                gCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTGGcgaaca                                                                                                                                      >  6:303807/1‑114 (MQ=38)                                                gCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTGGcgaaca                                                                                                                                      <  5:303807/119‑6 (MQ=38)                                                     ttCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCgg                                                                                                  >  6:184129/1‑149 (MQ=255)                                                       cGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggtg                                                                                                >  8:279686/1‑149 (MQ=255)                                                       cGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCggt                                                                                                 <  3:211718/148‑1 (MQ=255)                                                       cGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAG‑CAGGCGGCGCGCATCCTGGCTGCTcag                                                                                                  <  2:313270/147‑3 (MQ=255)                                                       cGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCACGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCCGAGCAG‑CAGGCAGCGCGCATCCTGGCTGCTcagg                                                                                                 <  4:227067/148‑4 (MQ=255)                                                               tGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCCGGCAACGCCTGAGCAGCCA‑GATGCGCGCCGCCT‑GCTGCTAGGTGGTGtgctgc                                                                                       >  1:117381/1‑149 (MQ=255)                                                                 gCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTc                                                                                                    <  7:144926/135‑1 (MQ=255)                                                                 gCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTc                                                                                                    >  8:144926/1‑135 (MQ=255)                                                                   caccaGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTgg                                                                                    >  4:240366/1‑149 (MQ=255)                                                                     ccaGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTggcc                                                                                  >  6:340683/1‑149 (MQ=255)                                                                     ccaGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTggcc                                                                                  <  3:240366/149‑1 (MQ=255)                                                                      caGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTggccg                                                                                 <  6:101834/149‑1 (MQ=255)                                                                      caGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTggccg                                                                                 <  8:63244/149‑1 (MQ=255)                                                                        gcgcCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTggccggc                                                                               <  5:340683/149‑1 (MQ=255)                                                                          gcCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTggccggcc                                                                              >  7:253068/1‑148 (MQ=255)                                                                          gcCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTggccggcc                                                                              <  8:253068/148‑1 (MQ=255)                                                                           ccAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCAc                                                                            >  4:46574/1‑149 (MQ=255)                                                                            cAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACt                                                                           >  6:110767/1‑149 (MQ=255)                                                                              gCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTgg                                                                         >  4:63034/1‑149 (MQ=255)                                                                               cccGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTggg                                                                        >  2:186654/1‑149 (MQ=255)                                                                                ccGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGt                                                                       <  7:47405/149‑1 (MQ=255)                                                                                   gCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAg                                                                    >  3:3085/1‑149 (MQ=255)                                                                                    cGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGc                                                                   <  7:279686/149‑1 (MQ=255)                                                                                         atcGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCTAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGcttct                                                               >  4:260087/1‑148 (MQ=255)                                                                                            gTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTAGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCa                                                            <  5:184129/148‑1 (MQ=255)                                                                                             tCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATc                                                          <  5:110767/149‑1 (MQ=255)                                                                                               ttCGGTCGCGCCACCTTTGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGc                                                        >  2:243124/1‑149 (MQ=255)                                                                                                   gTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGt                                                     >  5:167033/1‑148 (MQ=255)                                                                                                      gcgcCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTc                                                    >  5:285740/1‑146 (MQ=255)                                                                                                      gcgcCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTc                                                    <  6:285740/146‑1 (MQ=255)                                                                                                      gcgcCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCCGGCAACGCCTGAGCAGCCA‑GATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCaaa                                                 <  2:117381/148‑1 (MQ=255)                                                                                                          cACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTgg                                             <  7:124007/149‑1 (MQ=255)                                                                                                              ttCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGggcgg                                         >  1:114434/1‑149 (MQ=255)                                                                                                                cGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGggcggc                                        <  3:46574/148‑1 (MQ=255)                                                                                                                  ccGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGTGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGgcggcg                                       >  4:250821/1‑147 (MQ=255)                                                                                                                  ccGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGgcggcg                                       >  8:3612/1‑147 (MQ=255)                                                                                                                  ccGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTa                                     <  5:202830/149‑1 (MQ=255)                                                                                                                   cGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAg                                    >  7:42580/1‑149 (MQ=255)                                                                                                                   cGTACTCGCGCGGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCTACGCCTGAGCAGCCCGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAg                                    <  7:186824/149‑1 (MQ=255)                                                                                                                        tCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGt                                               <  8:170012/133‑1 (MQ=255)                                                                                                                        tCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGt                                               >  7:170012/1‑133 (MQ=255)                                                                                                                          gcgcTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCttcttc                                                             >  8:188707/1‑117 (MQ=255)                                                                                                                          gcgcTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCttcttc                                                             <  7:188707/117‑1 (MQ=255)                                                                                                                          gcgcTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATc                               <  3:63034/147‑1 (MQ=255)                                                                                                                                 cccTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTc                      >  2:213296/1‑149 (MQ=255)                                                                                                                                  ccTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGTGGCGATGGCTTCGCCGcccact                     <  1:243124/149‑4 (MQ=255)                                                                                                                                   cTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGGGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCgg                    <  1:320214/149‑1 (MQ=255)                                                                                                                                    tGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGc                   >  6:208672/1‑149 (MQ=255)                                                                                                                                        cGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCg               >  2:272701/1‑149 (MQ=255)                                                                                                                                         gCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGTGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGc              <  3:250821/149‑1 (MQ=255)                                                                                                                                         gCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCg               <  3:315100/148‑1 (MQ=255)                                                                                                                                         gCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGc              >  1:190922/1‑149 (MQ=255)                                                                                                                                         gCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGc              <  6:167033/149‑1 (MQ=255)                                                                                                                                          ccAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCTAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGcc             <  3:260087/149‑1 (MQ=255)                                                                                                                                           cAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGc              >  6:317885/1‑147 (MQ=255)                                                                                                                                                tGGCCGGTACGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATAGCCACCCTCGGCTTCGCCAGCTTc       >  7:170566/1‑149 (MQ=255)                                                                                                                                                 ggCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCCGCTTGTTCATCGCGGTCAAAGTGGGCGGCGTAGCCTTCGCCACCCGCGGGTTCGCCCGCTGcc      >  8:48073/1‑149 (MQ=255)                                                                                                                                                    cGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCgg   >  2:196048/1‑149 (MQ=255)                                                                                                                                                     ggTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGc  <  5:317885/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCT‑GCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGC  >  NC_002947/5125222‑5125518

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