Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,743,351	A→C	10.0%	intergenic (‑184/+63)	fbp ← / ← glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,743,351	0	A	C	10.0%	167.6 / 10.3	60	intergenic (‑184/+63)	fbp/glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (33/21);  new base C (2/4);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

ACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAA  >  NC_002947/5743204‑5743497                                                                                                                                                      |                                                                                                                                                   aCGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAg                                                                                                                                                   >  8:157202/1‑149 (MQ=255) aCGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAg                                                                                                                                                   >  2:334642/1‑149 (MQ=255)   gCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAgcg                                                                                                                                                 >  5:257736/1‑149 (MQ=255)   gCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAgcg                                                                                                                                                 >  5:257741/1‑149 (MQ=255)      gTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCaa                                                                                                                                              <  2:232570/149‑1 (MQ=255)            aGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCtt                                                                                                                                         >  5:129047/1‑148 (MQ=255)               gagaGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAg                                                                                                                                     >  5:259818/1‑149 (MQ=255)                agagTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGc                                                                                                                                    <  4:167805/149‑1 (MQ=255)                  agTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCtt                                                                                                                                  <  1:334642/149‑1 (MQ=255)                      aCGGGCTAGGGAGTGATGAGAGTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGTCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCTCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCg                                                                                                                               <  1:163810/148‑1 (MQ=255)                      aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCTTCAAGGCTTGCgc                                                                                                                              <  1:69064/149‑2 (MQ=255)                           cTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCgctgtac                                                                                                                         >  4:292091/1‑143 (MQ=25)                                gAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                                    <  8:166313/149‑1 (MQ=255)                                 aGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCGTTCGCGGGCACG‑‑CCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                   <  3:365565/149‑1 (MQ=18)                                 aGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                   <  3:292091/149‑1 (MQ=21)                                   tgatgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCACG‑‑CCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGa                                                                                                                 <  3:57614/149‑1 (MQ=18)                                     atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAg                                                                                                                <  6:129047/148‑1 (MQ=255)                                     atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGc                                                                                                               >  2:342238/1‑149 (MQ=255)                                      tgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCGTGAAGCTTGCGAGACACCTGTGGGAGCg                                                                                                              >  4:252593/1‑149 (MQ=255)                                        agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCgg                                                                                                             <  6:10619/148‑1 (MQ=255)                                             cGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCCCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGagcgagcgc                                                                                                         <  1:105725/149‑2 (MQ=17)                                               ccACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGcac                                                                                                     >  3:121427/1‑147 (MQ=255)                                                cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcg                                                                                                    >  4:335220/1‑149 (MQ=255)                                                cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcg                                                                                                    >  6:94651/1‑149 (MQ=255)                                                                     tAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCaa                                                                                                                                              >  4:99510/1‑86 (MQ=255)                                                                     tAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCaa                                                                                                                                              <  3:99510/86‑1 (MQ=255)                                                                        gTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCGTTGAAGCTGGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGCGCGGCCTCTCGTCAATCCAAt                                                                            >  2:158704/1‑149 (MQ=38)                                                                            gTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCt                                                                        >  8:105267/1‑149 (MQ=255)                                                                               gTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTc                                                                      <  7:105267/148‑1 (MQ=255)                                                                                         cGCCTGGCGGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑AGCGCTCCCCAAGGTACAGCGCAAGCCTTTAAGCTTGCGAGACACCTGTGGGAGCGGGGGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCgg                                                            <  1:342238/148‑1 (MQ=255)                                                                                          gCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGTAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                          <  3:335220/149‑1 (MQ=255)                                                                                          gCCAGGCCGGTGGGCTGGTTGGGGGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                          <  6:259818/149‑1 (MQ=255)                                                                                              ggCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATCCCTTCCAGATCTCGGTGGCat                                                      >  5:62783/1‑149 (MQ=255)                                                                                                  ggtggtCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTc                                                  <  7:157202/149‑1 (MQ=255)                                                                                                      gtCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGGGAGGCACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCTTATTcgcg                                               >  7:76998/1‑148 (MQ=255)                                                                                                           gTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                                         >  3:133239/1‑149 (MQ=255)                                                                                                                gtgtCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGCCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGt                                    >  8:116087/1‑149 (MQ=255)                                                                                                                   tCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                        >  3:26942/1‑142 (MQ=255)                                                                                                                   tCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                        <  4:26942/142‑1 (MQ=255)                                                                                                                     cATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                               <  5:94651/149‑1 (MQ=255)                                                                                                                      aTTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCAGCTCAATCCAATGCCTTCCAGATCTCGGTGGGATATTCGCGAATGGTGCGGTCcgac                               >  5:123946/1‑148 (MQ=255)                                                                                                                          gcgGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGaa                          >  1:40275/1‑149 (MQ=255)                                                                                                                            gggCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                <  3:261930/93‑1 (MQ=255)                                                                                                                            gggCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                >  4:261930/1‑93 (MQ=255)                                                                                                                             gcggcg‑‑cTCCCTCCCACACGTACACCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCa                       <  8:76998/146‑1 (MQ=255)                                                                                                                              gcacg‑‑cccgctCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGAGCGAGCGCGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAg                      >  6:161439/4‑149 (MQ=17)                                                                                                                                gcg‑‑cTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAAc                         >  5:66774/1‑144 (MQ=255)                                                                                                                                gcg‑‑cTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAAc                         <  6:66774/144‑1 (MQ=255)                                                                                                                                       cgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTCCGAGACACATGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATcc               <  3:72106/149‑1 (MQ=255)                                                                                                                                        gctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccg              >  3:107289/1‑149 (MQ=255)                                                                                                                                          tcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcg                                                 >  6:320522/1‑112 (MQ=255)                                                                                                                                          tcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcg                                                 <  5:320522/112‑1 (MQ=255)                                                                                                                                          tcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtg            >  3:170462/1‑149 (MQ=255)                                                                                                                                            ccACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgcc          >  5:379820/1‑149 (MQ=255)                                                                                                                                            ccACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgcc          >  8:186032/1‑149 (MQ=255)                                                                                                                                              acaGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccgt        >  5:379083/1‑149 (MQ=255)                                                                                                                                              acaGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCACATccgtgccgt        >  4:83179/1‑149 (MQ=255)                                                                                                                                                  gTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt     >  3:210430/1‑148 (MQ=255)                                                                                                                                                   tACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTc    >  5:306867/1‑148 (MQ=255)                                                                                                                                                    aCAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGCGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCACTGCCGTGTTCaa  >  5:47/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                   ACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAA  >  NC_002947/5743204‑5743497

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