Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	433,254	C→G	5.6%	G649A (GGC→GCC)	glcB ←	malate synthase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	433,254	0	C	G	5.6%	229.2 / 2.8	71	G649A (GGC→GCC)	glcB	malate synthase G
Reads supporting (aligned to +/- strand):  ref base C (36/31);  new base G (2/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TACCAGTTCCACAGCCGCCTGGAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATCCAGCGCACCACGTAGCCGAGAATGCCCTGGG  >  NC_002947/433109‑433401                                                                                                                                                    |                                                                                                                                                    tACCAGTTCCACAGCCGCCTGGAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGt                                                                                                                                                  >  6:302355/1‑149 (MQ=255)   ccAGTTCCACAGCCGCCTGGAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTg                                                                                                                                                 <  2:281256/148‑1 (MQ=255)      gatCCACAGGCGCCTGGAACGCCACGTTGTCGTCGAAGTTGTGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGGTTG‑AGGCTTTCCAGCACGTGTTCCTGGCTG‑ACCACGCCGTGgcgc                                                                                                                                             <  2:234539/147‑1 (MQ=255)         ccACAGCCGCCTGGAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTGCGCCACGGCGTggtcagc                                                                                                                                           >  8:202981/1‑144 (MQ=255)           acaGCCGCCTGGAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTGCGCCACGGCGTggtcagcca                                                                                                                                         <  7:202981/149‑8 (MQ=255)            caGCCGCCTGGAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAg                                                                                                                                       >  8:42741/1‑149 (MQ=255)              gccgccTGGAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGtt                                                                                                                                     >  5:351784/1‑149 (MQ=255)               ccgccTGGAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTg                                                                                                                                    <  7:42741/149‑1 (MQ=255)                      gAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTGCAGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCcag                                                                                                                               <  1:104141/148‑1 (MQ=255)                      gAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCcagca                                                                                                                             <  5:102433/149‑1 (MQ=255)                      gAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCcagc                                                                                                                              <  5:242626/148‑1 (MQ=255)                        aCGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCagcagc                                                                                                                           <  4:49522/149‑1 (MQ=255)                            cACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTggg                                                                                                                       >  4:240194/1‑149 (MQ=255)                            cACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTggg                                                                                                                       >  6:251245/1‑149 (MQ=255)                            cACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGATggg                                                                                                                       <  3:229120/149‑1 (MQ=255)                                     gtcgAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATgc                                                                                                              >  2:124966/1‑149 (MQ=255)                                              gggCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTg                                                                                                      <  7:43994/148‑1 (MQ=255)                                               ggCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTgg                                                                                                     >  1:284557/1‑148 (MQ=255)                                                  gcCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACg                                                                                                 >  4:140600/1‑149 (MQ=255)                                                    cATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGt                                                                                               >  6:299020/1‑149 (MQ=255)                                                    cATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGt                                                                                               >  2:200778/1‑149 (MQ=255)                                                    cATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTAGCACGGt                                                                                               >  3:189641/1‑149 (MQ=255)                                                     aTCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATgcgc                                                                                                            >  6:162376/1‑135 (MQ=255)                                                     aTCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATgcgc                                                                                                            <  5:162376/135‑1 (MQ=255)                                                            ggTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCAt                                                                                        >  8:57613/1‑148 (MQ=255)                                                                 aGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGcc                                                                                  <  1:124966/149‑1 (MQ=255)                                                                       tCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACgtt                                                                            <  5:302355/149‑1 (MQ=255)                                                                       tCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACgtt                                                                            >  5:286072/1‑149 (MQ=255)                                                                          ccAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACgttgtt                                                                         <  3:240194/149‑1 (MQ=255)                                                                          ccAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACgttgtt                                                                         <  5:251245/149‑1 (MQ=255)                                                                           cAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGttgttg                                                                        <  3:140600/149‑1 (MQ=255)                                                                           cAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGttgttg                                                                        <  7:85183/149‑1 (MQ=255)                                                                             gCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGAt                                                                      <  1:200778/149‑1 (MQ=255)                                                                             gCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTCCCATCAGGCCGACGTTGTTGAt                                                                      <  7:57613/149‑1 (MQ=255)                                                                                ttCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTc                                                                   <  2:284557/149‑1 (MQ=255)                                                                                     ctgATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCa                                                              >  2:329899/1‑149 (MQ=255)                                                                                     ctgATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTAGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCgg                                                                <  4:189641/147‑1 (MQ=255)                                                                                      tgATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCAc                                                             >  7:225295/1‑149 (MQ=255)                                                                                           gACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTc                                                         >  5:359605/1‑148 (MQ=255)                                                                                           gACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTc                                                         >  5:56510/1‑148 (MQ=255)                                                                                            accaccACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCATGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTc                                                         >  6:145077/1‑147 (MQ=255)                                                                                            accaccACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGtt                                                                                                                                     >  1:202136/1‑71 (MQ=255)                                                                                             ccaccaCGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTgg                                                                                                                        >  3:361688/1‑83 (MQ=255)                                                                                             ccaccaCGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTgg                                                                                                                        <  4:361688/83‑1 (MQ=255)                                                                                             ccaccaCGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATgcg                                                                                                             <  5:156314/94‑1 (MQ=255)                                                                                             ccaccaCGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATgcg                                                                                                             >  6:156314/1‑94 (MQ=255)                                                                                               accacGGCCATGCGCTTG‑AGGCTTTCCAGCACTTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGc                                                     >  3:150433/1‑148 (MQ=255)                                                                                                ccacGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGc                                                               >  6:193850/1‑137 (MQ=255)                                                                                                ccacGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGc                                                               <  5:193850/137‑1 (MQ=255)                                                                                                ccacGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGc                                                     >  6:3188/1‑147 (MQ=255)                                                                                                ccacGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGc                                                     <  5:3188/147‑1 (MQ=255)                                                                                                  acGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAAcc                                                 >  4:263729/1‑149 (MQ=255)                                                                                                         tGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGttgttg                                                                        <  1:379000/119‑1 (MQ=255)                                                                                                         tGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGttgttg                                                                        >  2:379000/1‑119 (MQ=255)                                                                                                         tGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGccc                                          <  1:329899/149‑1 (MQ=255)                                                                                                             cTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGt                                      <  6:215176/149‑1 (MQ=255)                                                                                                                g‑AGGCTTTCCAGCACCTGTTCCGGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCg                                    >  7:241173/1‑148 (MQ=255)                                                                                                                    gCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCa                                                    <  6:370706/129‑1 (MQ=255)                                                                                                                    gCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCa                                                    >  5:370706/1‑129 (MQ=255)                                                                                                                                cTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATCCAGCGCACCACg                     >  5:109756/1‑148 (MQ=255)                                                                                                                                   ttCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATc                                 <  6:92/133‑1 (MQ=255)                                                                                                                                   ttCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATc                                 >  5:92/1‑133 (MQ=255)                                                                                                                                         gCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATCCAGCGCACCACGTAGCCgaga            <  6:351784/148‑1 (MQ=255)                                                                                                                                         gCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATCCAGCGCACCACGTAGCCGAGaa           >  6:53366/1‑149 (MQ=255)                                                                                                                                           tG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATCCAGCGCACCACGTAGCCGAGAATg         >  8:179198/1‑149 (MQ=255)                                                                                                                                            gcgccACGGCGTGG‑TCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATCCAGCGCACCACGTAGc                 >  6:180795/4‑140 (MQ=255)                                                                                                                                            gcgccACGGCGTGG‑TCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATCCAGCGCACCACGTAGc                 <  5:180795/137‑1 (MQ=255)                                                                                                                                                cACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATCCAGCGCACCACGTAGCCGAGaa           >  7:118177/1‑143 (MQ=255)                                                                                                                                                cACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATCCAGCGCACCACGTAGCCGAGaa           <  8:118177/143‑1 (MQ=255)                                                                                                                                                cACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATCCAGCGCACCACGTAGCCGAGAATGCCCt     >  4:348342/1‑149 (MQ=255)                                                                                                                                                   gCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGCGCAACCCACGCCCTGGTCGATCCAGCGCACCACGTAGCCGAGAATGCCCTggg  >  3:250520/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TACCAGTTCCACAGCCGCCTGGAACGCCACGTTGTCGTCGAAGTTGGGCGCCATCGGGCGGTACAGGGGGTCGCCAGCGTTCTGCTGATCGACCACCACGGCCATGCGCTTG‑AGGCTTTCCAGCACCTGTTCCTGGCTG‑ACCACGCCGTGGCGCAGCCAGTTGGCCAGCAGCTGGGCAGAAATGCGCAGGGTGGCACGGTCTTCCATCAGGCCGACGTTGTTGATGTCCGGCACCTTCGAGCAACCCACGCCCTGGTCGATCCAGCGCACCACGTAGCCGAGAATGCCCTGGG  >  NC_002947/433109‑433401

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