Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,893,834	A→C	7.7%	F409V (TTC→GTC)	yegH ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,893,834	0	A	C	7.7%	232.6 / 7.5	75	F409V (TTC→GTC)	yegH	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base C (3/3);  total (41/34)
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

CGCGCCTGCACCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCC  >  NC_002947/3893686‑3893982                                                                                                                                                      |                                                                                                                                                     cgcgCCTGCACCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGaa                                                                                                                                                      <  2:85169/149‑1 (MQ=255)      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                                                                                                                                       >  1:228948/1‑146 (MQ=255)                   aGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGtt                                                                                                                                       <  2:228948/146‑1 (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                                                                                    >  6:156485/1‑142 (MQ=255)                                                                          cTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCAtt                                                                                    <  5:156485/142‑1 (MQ=255)                                                                                  tCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGATGATGCTGAAGct                                                                     >  3:250107/1‑149 (MQ=255)                                                                                     cTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCtctc                                                                  <  8:235419/149‑1 (MQ=255)                                                                                     cTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCtctc                                                                  <  1:344419/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                                                                          <  5:215686/138‑1 (MQ=255)                                                                                        gCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTg                                                                          >  6:215686/1‑138 (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              <  7:245320/149‑1 (MQ=255)                                                                                                                                         aC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGccc              <  8:286899/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)                                                                                                                                                      |                                                                                                                                                     CGCGCCTGCACCTGGCTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCC  >  NC_002947/3893686‑3893982

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