Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,743,363	G→A	16.0%	intergenic (‑196/+51)	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,363	0	G	A	16.0%	135.7 / 15.3	57	intergenic (‑196/+51)	fbp/glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (29/18);  new base A (4/5);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01

GAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGC  >  NC_002947/5743218‑5743508                                                                                                                                                    |                                                                                                                                                  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                                                                                                               >  6:94651/1‑149 (MQ=255)                                  cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcg                                                                                                               >  4:335220/1‑149 (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                                                   <  4:26942/142‑1 (MQ=255)                                                                                                     tCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                   >  3:26942/1‑142 (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                                                                                           >  4:261930/1‑93 (MQ=255)                                                                                                              gggCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                           <  3:261930/93‑1 (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                                                            <  5:320522/112‑1 (MQ=255)                                                                                                                            tcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcg                                                            >  6:320522/1‑112 (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)                                                                                                                                      gtctcgcAAGCTTCAAGGCTTGCGCGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCaa             >  7:104062/5‑149 (MQ=38)                                                                                                                                      gtctcgcAAGCTTCAAGGCTTGCGCGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCaa             >  7:98887/5‑149 (MQ=38)                                                                                                                                      aCAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGCGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCACTGCCGTGTTCaa             >  5:47/1‑149 (MQ=255)                                                                                                                                          cgcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATAcc         >  5:58828/1‑149 (MQ=255)                                                                                                                                           gcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCg        >  4:313838/1‑149 (MQ=255)                                                                                                                                                 cTTGACGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgc  <  3:252593/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGC  >  NC_002947/5743218‑5743508

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