Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,893,964	G→C	10.1%	S365R (AGC→AGG) ‡	yegH ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,893,964	0	G	C	10.1%	205.1 / 12.4	69	S365R (AGC→AGG) ‡	yegH	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base C (3/4);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

TGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAGGT  >  NC_002947/3893816‑3894109                                                                                                                                                     |                                                                                                                                                  tGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctgg                                                                                                                                                   >  7:170786/1‑149 (MQ=255) tGGTCAACAGCCCGGTGAAGTCAACGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctgg                                                                                                                                                   <  6:356065/149‑1 (MQ=255)   gTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctggct                                                                                                                                                 >  6:166784/1‑149 (MQ=255)   gTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctggct                                                                                                                                                 >  2:40441/1‑149 (MQ=255)   gTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctggc                                                                                                                                                  <  8:238781/148‑1 (MQ=255)         aGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGTGCCATGCTTTCCAGGTCGGGCTGGCTGCCCga                                                                                                                                           <  1:40441/149‑1 (MQ=255)           cccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCAGCCAGCCCgacc                                                                                                                                         <  5:109451/149‑3 (MQ=255)           cccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCAGCCAGCCCgacc                                                                                                                                         <  6:101674/149‑3 (MQ=255)            ccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGagag                                                                                                                                        <  3:135060/149‑1 (MQ=255)            ccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCAGCCAGCCCgacct                                                                                                                                        <  5:5698/149‑4 (MQ=255)                       aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt                                                                                                                               >  2:213901/1‑147 (MQ=255)                       aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt                                                                                                                               >  3:139334/1‑147 (MQ=255)                       aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt                                                                                                                               >  4:14175/1‑147 (MQ=255)                        ccGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcag                                                                                                                            <  8:170786/149‑1 (MQ=255)                         cGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcag                                                                                                                            >  8:241060/1‑148 (MQ=255)                          gAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGACagcagttccttgag                                                                                                                            >  1:192510/1‑147 (MQ=255)                                ggtGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagcagttcc                                                                                                                     <  6:382033/146‑1 (MQ=255)                                gTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagcagttcc                                                                                                                     <  5:252742/148‑1 (MQ=255)                                gTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagcagttcc                                                                                                                     <  3:196034/148‑1 (MQ=255)                                   gACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagcagttccttc                                                                                                                  <  6:197116/148‑1 (MQ=255)                                      cACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTtgt                                                                                                              <  5:269848/149‑1 (MQ=255)                                       aCGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTtgtg                                                                                                             <  4:139334/149‑1 (MQ=255)                                        cGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGc                                                                                                            <  2:156472/149‑1 (MQ=255)                                                 tatGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGcc                                                                                                    >  2:351378/1‑148 (MQ=255)                                                  atGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGccc                                                                                                   >  8:325651/1‑148 (MQ=255)                                                     ggTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAg                                                                                                <  6:347992/148‑1 (MQ=255)                                                     ggTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAg                                                                                                <  3:366766/148‑1 (MQ=255)                                                              aCCGCGCATTTGCTCCAGGGCAGTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCa                                                                                      >  8:81306/1‑149 (MQ=255)                                                                 gcgcATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCAcgc                                                                                   <  3:100514/149‑1 (MQ=255)                                                                 gcgcATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCAcgc                                                                                   >  5:334298/1‑149 (MQ=255)                                                                    cATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGAc                                                                                <  2:130395/149‑1 (MQ=255)                                                                    cATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGGTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCGGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGAc                                                                                <  6:276790/149‑1 (MQ=255)                                                                     aTTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGa                                                                                 <  2:276708/147‑1 (MQ=255)                                                                     aTTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTACAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGa                                                                                 >  1:276708/1‑147 (MQ=255)                                                                         gCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCGGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATc                                                                            <  4:145388/148‑1 (MQ=255)                                                                          cTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgc                                                                          <  5:125973/149‑1 (MQ=255)                                                                           tCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                         >  6:144623/1‑149 (MQ=255)                                                                           tCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                         >  5:305523/1‑149 (MQ=255)                                                                             cAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGa                                                                       >  4:252167/1‑149 (MQ=255)                                                                                  cATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGc                                                                  <  5:166784/149‑1 (MQ=255)                                                                                    ttGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCgg                                                                >  5:13759/1‑149 (MQ=255)                                                                                      gAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCgg                                                                >  5:103106/1‑147 (MQ=255)                                                                                       aGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCa                                                              <  6:334298/148‑1 (MQ=255)                                                                                        ggATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAgg                                                            >  7:7212/1‑149 (MQ=255)                                                                                        ggATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAgg                                                            >  1:99770/1‑149 (MQ=255)                                                                                             cTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                       <  1:213901/149‑1 (MQ=255)                                                                                             cTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                       <  7:241060/149‑1 (MQ=255)                                                                                             cTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                       <  3:252167/149‑1 (MQ=255)                                                                                              tGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAg                                                      >  6:180706/1‑149 (MQ=255)                                                                                                aaGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTa                                                    >  1:55747/1‑149 (MQ=255)                                                                                                    tctcCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGc                                                >  3:228032/1‑149 (MQ=255)                                                                                                      tcCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGa                                                                       <  3:366347/124‑1 (MQ=255)                                                                                                      tcCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGa                                                                       >  4:366347/1‑124 (MQ=255)                                                                                                        cagcagGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGaa                                             <  3:14175/148‑1 (MQ=255)                                                                                                        cagcagGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAAt                                            >  3:323611/1‑149 (MQ=255)                                                                                                            agGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGc                                        <  4:228032/149‑1 (MQ=255)                                                                                                              gTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCa                                      <  8:7212/149‑1 (MQ=255)                                                                                                                tCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAgg                                    >  7:345457/1‑149 (MQ=255)                                                                                                                 cAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAgg                                    <  1:351378/148‑1 (MQ=255)                                                                                                                   aCGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGcc                                 >  1:85131/1‑149 (MQ=255)                                                                                                                   aCGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCGTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGcc                                 >  6:134717/1‑149 (MQ=255)                                                                                                                     ggCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTg                               >  2:61878/1‑149 (MQ=255)                                                                                                                          ccGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCaa                          >  4:114671/1‑149 (MQ=255)                                                                                                                            gggCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATg                        <  7:72787/149‑1 (MQ=255)                                                                                                                                     tttCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGcc               >  7:104479/1‑149 (MQ=255)                                                                                                                                           tctCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGGGAATTGGCCAGGGCCTGGGCAATGGCGTCGGGCGGTTgc         >  4:310943/3‑149 (MQ=255)                                                                                                                                           ggTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgc         <  1:61878/149‑1 (MQ=255)                                                                                                                                           ggTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgc         <  5:180706/149‑1 (MQ=255)                                                                                                                                           ggTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgc         <  8:1163/149‑1 (MQ=255)                                                                                                                                            ctCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcg        >  5:343820/2‑149 (MQ=255)                                                                                                                                             tCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcgc       <  2:55747/149‑1 (MQ=255)                                                                                                                                             tCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGc                                                                  >  4:213015/1‑90 (MQ=255)                                                                                                                                             tCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGc                                                                  <  3:213015/90‑1 (MQ=255)                                                                                                                                              cGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGGAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcgc       >  6:182900/1‑148 (MQ=255)                                                                                                                                                  ctggctgCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAGGt  <  6:305523/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAGGT  >  NC_002947/3893816‑3894109

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