Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,893,835	G→T	6.5%	D408E (GAC→GAA)	yegH ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,893,835	0	G	T	6.5%	247.9 / 6.2	77	D408E (GAC→GAA)	yegH	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base T (2/3);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

CTGCACCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCT  >  NC_002947/3893691‑3893983                                                                                                                                                  |                                                                                                                                                     cTGCACCTGGCTGAGGTTCAGCGCACCGATGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCAc                                                                                                                                                  >  1:210515/1‑149 (MQ=255)    cACCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGAGGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCAcc                                                                                                                                                 <  6:93377/147‑1 (MQ=255)     aCCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGaa                                                                                                                                              >  1:240139/1‑149 (MQ=255)       cTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACt                                                                                                                                            >  5:236880/1‑149 (MQ=255)           cTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGtt                                                                                                                                        <  1:6487/149‑1 (MQ=255)            tGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCAcc                                                                                                                                                 >  3:41273/1‑139 (MQ=255)            tGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCAcc                                                                                                                                                 <  4:41273/139‑1 (MQ=255)            tGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATACGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGACATGGTAAAC‑AGCCCGCTCAAGTCACCGAAATCGtt                                                                                                                                        >  7:253618/1‑148 (MQ=255)              aGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGtt                                                                                                                                        <  2:228948/146‑1 (MQ=255)              aGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGtt                                                                                                                                        >  1:228948/1‑146 (MQ=255)                gTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCa                                                                                                                                   >  2:35555/1‑149 (MQ=255)                 ttCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCat                                                                                                                                  <  3:212457/149‑2 (MQ=37)                  tCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAT‑AGCCCGGTGAAGTCACCGAACTCGTTGACCAc                                                                                                                                  >  1:346940/1‑148 (MQ=255)                  tCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCAc                                                                                                                                  >  1:264688/1‑148 (MQ=255)                    aGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCATGa                                                                                                                                <  2:363340/148‑1 (MQ=37)                    aGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGaa                                                                                                                               <  2:99897/149‑1 (MQ=255)                     gcgcACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCAcc                                                                                                                                                 <  4:8132/130‑1 (MQ=255)                     gcgcACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCAcc                                                                                                                                                 >  3:8132/1‑130 (MQ=255)                       gcACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGc                                                                                                                            <  2:240139/149‑1 (MQ=255)                        cACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCat                                                                                                                                  <  3:100155/142‑2 (MQ=37)                        cACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCat                                                                                                                                  >  4:100155/1‑141 (MQ=37)                          ccGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAat                                                                                                                         <  4:240520/149‑1 (MQ=255)                            gCTGACCACAAAGCCTTCCGCCTCCTGAACGATGCCCTGGGCTTCAACCTCACTGGCATGAGGCAGCTCACCGGGAATGGACTCAAGGATGTCGGTCATGGTGTTC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGTAGGCAatat                                                                                                                       <  8:39677/149‑1 (MQ=255)                            gCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAT‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAatat                                                                                                                       <  2:346940/149‑1 (MQ=255)                                     aaaGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATTGGTCGACt                                                                                                              >  3:183640/1‑149 (MQ=255)                                     aaaGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGGCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACt                                                                                                              >  3:267471/1‑149 (MQ=255)                                           ttCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACcgc                                                                                                        >  7:286899/1‑149 (MQ=255)                                               cccTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCAt                                                                                                    >  4:361595/1‑149 (MQ=255)                                                cctcctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCAtt                                                                                                   >  7:240782/1‑149 (MQ=255)                                                cctcctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCAtt                                                                                                   <  8:62886/149‑1 (MQ=255)                                                cctcctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCAtt                                                                                                   <  2:264688/149‑1 (MQ=255)                                                 ctcctGAACGATGCCAGGGCCTTCAACATCAGTGGCATCCGGCAGCGCACCGGCAATCGAGTCAAGCATGTCGGTAATGGTCAAC‑AGCCCGGTGAAGACACCGAACTCGTTGACCACGCAGGCAAGATGGGTCGCCTGACCGCGCAttt                                                                                                  >  1:279723/1‑149 (MQ=255)                                                    ctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCt                                                                                               >  2:344419/1‑149 (MQ=255)                                                    ctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGAATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCt                                                                                               >  7:235419/1‑149 (MQ=255)                                                      gAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTc                                                                                              >  2:73142/1‑148 (MQ=255)                                                          gATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAggg                                                                                         >  4:231071/1‑149 (MQ=255)                                                              cccGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCAtt                                                                                     <  1:226834/149‑1 (MQ=255)                                                                  ggCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAgg                                                                                 <  8:240782/149‑1 (MQ=255)                                                                     cTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCAtt                                                                                     <  5:156485/142‑1 (MQ=255)                                                                     cTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCAtt                                                                                     >  6:156485/1‑142 (MQ=255)                                                                             tCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGATGATGCTGAAGct                                                                      >  3:250107/1‑149 (MQ=255)                                                                                cTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCtctc                                                                   <  1:344419/149‑1 (MQ=255)                                                                                cTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCtctc                                                                   <  8:235419/149‑1 (MQ=255)                                                                                 tGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTc                                                                                              <  4:49309/121‑1 (MQ=255)                                                                                 tGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTc                                                                                              >  3:49309/1‑121 (MQ=255)                                                                                 tGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCtctc                                                                   >  6:109809/1‑148 (MQ=255)                                                                                   gCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTg                                                                           >  6:215686/1‑138 (MQ=255)                                                                                   gCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTg                                                                           <  5:215686/138‑1 (MQ=255)                                                                                     aTCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATTGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCcagca                                                              <  4:183640/149‑1 (MQ=255)                                                                                      tCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGa                                                                                                                >  4:244546/1‑98 (MQ=255)                                                                                      tCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGa                                                                                                                <  3:244546/98‑1 (MQ=255)                                                                                         ggCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGACGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGCAGCTCACCAGCAGGtt                                                          >  2:315850/1‑149 (MQ=255)                                                                                                 aCCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGgcg                                                  >  6:269094/1‑149 (MQ=255)                                                                                                   cGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGcc                                                >  8:245320/1‑149 (MQ=255)                                                                                                       aaTCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCggg                                            >  6:37961/1‑149 (MQ=255)                                                                                                           gACTCAAGGATGTCGGTCATGGTTAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCAt                                        >  3:122675/1‑149 (MQ=255)                                                                                                                aaGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACTACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCttt                                   >  1:176576/1‑149 (MQ=255)                                                                                                                aaGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCttt                                   >  6:150813/1‑149 (MQ=255)                                                                                                                aaGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCttt                                   >  1:305179/1‑149 (MQ=255)                                                                                                                  ggATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTc                                  >  4:356649/1‑148 (MQ=255)                                                                                                                      gTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGt                             >  8:175336/1‑149 (MQ=255)                                                                                                                          ttcgtGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCATGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCGAGCAGGTTCAACGGCGCCCGGGCCATGCTTGCCAGGTCgcg                         >  4:264668/5‑147 (MQ=25)                                                                                                                          gTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCggg                         <  1:35555/149‑1 (MQ=255)                                                                                                                           tCATCGTCAAC‑AGCCCCGTGAAGTCACCGACCTCGTTGACCACGAAGGCAATCTGGGTCGACTGACCGCGCATTTGCGGCTGTGCATTTAGCCTCCTGCAGCTCTCCAGCAGGCTCAACGTCGCCCGGGCCATGCTTTCCAGGTCGggc                        <  2:279723/149‑1 (MQ=255)                                                                                                                              tGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctgg                     <  6:236880/149‑1 (MQ=255)                                                                                                                                gTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctggct                   >  2:162017/1‑149 (MQ=255)                                                                                                                                   aaC‑AGCCCGGTGAAGTCAACGAACTCGTTGACAACGAAGGCAATATTTGTCGACTGACCGCGCATTTGCTCAAGGGCATTGAGGATGCTGATGCTATCCAGCAGGTTAAAAGGAGCCCTGGCCATGCTTTCCAGGTCGGGCTTGCTGcc                >  2:233898/1‑149 (MQ=255)                                                                                                                                    aC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGccc               <  8:286899/149‑1 (MQ=255)                                                                                                                                    aC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGccc               <  7:245320/149‑1 (MQ=255)                                                                                                                                     c‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGccc               >  5:359673/1‑148 (MQ=255)                                                                                                                                       aGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCga             <  1:162017/149‑1 (MQ=255)                                                                                                                                         cccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTCCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCAGCCAGCCCgacc           <  3:231071/149‑3 (MQ=255)                                                                                                                                          ccgtttaaGTCACCGAACTCGTTGACAACGAAGGCAATATGGGGGGACTGACCGCGCATTTGCTCCAGGGCATTGTGGATGCTGAAGCTCTCCTGCAGGTTCAACTGCGACCGGGCCATTCTTTCCAGGTCGGGCTGTCTGCCCGagag          <  8:253618/143‑1 (MQ=255)                                                                                                                                          ccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGagag          <  6:359673/149‑1 (MQ=255)                                                                                                                                            ggTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagca        <  7:175336/149‑1 (MQ=255)                                                                                                                                            ggTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagca        <  3:361595/149‑1 (MQ=255)                                                                                                                                                 aGTCACCGAACTCGTTGACTACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcc   <  2:176576/149‑1 (MQ=255)                                                                                                                                                  gTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcct  >  4:296480/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CTGCACCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCT  >  NC_002947/3893691‑3893983

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