Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,893,966	T→A	10.1%	S365C (AGC→TGC) ‡	yegH ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,893,966	0	T	A	10.1%	186.7 / 10.8	69	S365C (AGC→TGC) ‡	yegH	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base A (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 = 7.21e-01

GTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAGGT  >  NC_002947/3893818‑3894109                                                                                                                                                     |                                                                                                                                                gTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctggct                                                                                                                                                 >  6:166784/1‑149 (MQ=255) gTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctggct                                                                                                                                                 >  2:40441/1‑149 (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                                                                                                                               >  4:14175/1‑147 (MQ=255)                     aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt                                                                                                                               >  2:213901/1‑147 (MQ=255)                     aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt                                                                                                                               >  3:139334/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                                                                                   >  5:334298/1‑149 (MQ=255)                                                               gcgcATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCAcgc                                                                                   <  3:100514/149‑1 (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                                                                         >  5:305523/1‑149 (MQ=255)                                                                         tCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                         >  6:144623/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                                                            >  1:99770/1‑149 (MQ=255)                                                                                      ggATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAgg                                                            >  7:7212/1‑149 (MQ=255)                                                                                           cTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                       <  7:241060/149‑1 (MQ=255)                                                                                           cTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                       <  1:213901/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                                                                       >  4:366347/1‑124 (MQ=255)                                                                                                    tcCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGa                                                                       <  3:366347/124‑1 (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         <  8:1163/149‑1 (MQ=255)                                                                                                                                         ggTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgc         <  1:61878/149‑1 (MQ=255)                                                                                                                                         ggTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgc         <  5:180706/149‑1 (MQ=255)                                                                                                                                          ctCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcg        >  5:343820/2‑149 (MQ=255)                                                                                                                                           tCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcgc       <  2:55747/149‑1 (MQ=255)                                                                                                                                           tCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGc                                                                  <  3:213015/90‑1 (MQ=255)                                                                                                                                           tCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGc                                                                  >  4:213015/1‑90 (MQ=255)                                                                                                                                            cGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGGAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcgc       >  6:182900/1‑148 (MQ=255)                                                                                                                                                ctggctgCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAGGt  <  6:305523/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                GTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAGGT  >  NC_002947/3893818‑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: ‑