Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,655,631	T→G	8.4%	intergenic (‑86/+66)	epd ← / ← tktA	D‑erythrose 4‑phosphate dehydrogenase/transketolase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,655,631	0	T	G	8.4%	161.2 / 7.1	60	intergenic (‑86/+66)	epd/tktA	D‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base G (3/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

GACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGC  >  NC_002947/5655486‑5655778                                                                                                                                                   |                                                                                                                                                      gACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCAt                                                                                                                                                    <  1:14941/149‑1 (MQ=255) gACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCAt                                                                                                                                                    <  7:267591/149‑1 (MQ=255) gACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCAt                                                                                                                                                    >  4:138237/1‑149 (MQ=255)  acacAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACTGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCAt                                                                                                                                                    <  5:171230/148‑1 (MQ=255)   cacaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGc                                                                                                                                                  <  2:310012/149‑1 (MQ=255)     caGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCaa                                                                                                                                                <  6:177321/149‑1 (MQ=255)      aGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAAc                                                                                                                                               >  8:254163/1‑149 (MQ=255)            cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCtgt                                                                                                                                         <  5:298219/149‑1 (MQ=255)            cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGCGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCtgt                                                                                                                                         <  8:144322/149‑1 (MQ=255)                 cGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGa                                                                                                                                     <  3:138237/148‑1 (MQ=255)                     cATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCggg                                                                                                                                <  6:128199/149‑1 (MQ=255)                                aGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGTAATTTGCGGGGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGGGGGAGCGGGCTTGCCcgcg                                                                                                                      <  8:21141/148‑1 (MQ=255)                                 gTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCgg                                                                                                                                 >  4:315385/1‑136 (MQ=37)                                 gTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCgg                                                                                                                                 <  3:315385/136‑1 (MQ=37)                                    cAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGAGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                                  >  2:300917/1‑148 (MQ=255)                                       cTTTGTAGGGACGGGGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGGGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGTGCGGGCTTGCCCGCGAATG‑AAgg                                                                                                             <  3:214992/149‑1 (MQ=255)                                        tttGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGc                                                                                                            >  3:124452/1‑149 (MQ=255)                                         ttGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGc                                                                                                            <  4:124452/148‑1 (MQ=255)                                            tAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGCCAAGCCCGCTCCCACAG‑GGCTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGaca                                                                                                        <  3:286059/149‑1 (MQ=255)                                               ggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAt                                                                                                                   >  7:282261/1‑136 (MQ=255)                                               ggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAt                                                                                                                   <  8:282261/136‑1 (MQ=255)                                                        ggCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCACCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCCACACCGCCCGCaca                                                                                             >  1:199074/1‑146 (MQ=255)                                                         gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATc                                                                                           <  2:238016/149‑1 (MQ=255)                                                            ttCGGGGCTCACGAGACATTTCACTTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCttt                                                                                        <  1:41028/149‑1 (MQ=255)                                                            ttCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCtt                                                                                         >  5:279229/1‑148 (MQ=255)                                                                gggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCGGTGGGAGCGGGCGTGCCCGCGAATG‑AAGGCGCCCCCGCCCGCAGCTCTTTaaa                                                                                     >  7:265299/1‑148 (MQ=38)                                                                 ggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAAt                                                                                    >  1:326517/1‑148 (MQ=255)                                                                 ggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGAGCTTTTAAt                                                                                    >  7:304581/1‑148 (MQ=255)                                                                  gCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCa                                                                                  <  5:351928/149‑1 (MQ=255)                                                                     cACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTc                                                                               <  3:160827/149‑1 (MQ=255)                                                                        gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACCG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGGCCCCGCCCGCAGCTCGTTAAAtagatcg                                                                              >  4:301847/1‑141 (MQ=37)                                                                        gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTaaa                                                                                     <  3:301847/140‑1 (MQ=255)                                                                           aCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCATGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCca                                                                          <  1:300917/148‑1 (MQ=255)                                                                           aCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCcag                                                                         <  6:343073/149‑1 (MQ=255)                                                                                       tcgtcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑gg                                                             <  1:293880/148‑1 (MQ=255)                                                                                          tcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAATCCCCCTACCACAA‑GGTTTAATTCAACCACTGTAGGAGCGGGCTTTGCCCCGAATG‑AAGGCGACACCGCCCCCCGATCTTTAAATAATTATTACACAAGCTCATA‑GGCAg                                                          >  8:306973/1‑148 (MQ=255)                                                                                               cTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACCCCGCCCGCAGCTCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTAcc                                                      >  5:167987/1‑147 (MQ=255)                                                                                                  cGGAATTTGCAGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAg                                                   <  8:251770/147‑1 (MQ=255)                                                                                                        ttGCGGTGTGGCTTTCGCGTGAAAGCCCGATCCCACAG‑GTTTTCATGAAACCTCTTTGTGAGCGGGCTTGCCCGCGAATG‑AAGGCTACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTtc                                            <  7:306973/148‑1 (MQ=38)                                                                                                          gCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTtctc                                          >  4:265181/1‑148 (MQ=255)                                                                                                          gCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCGCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATAAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGttggc                                          >  8:18773/1‑145 (MQ=255)                                                                                                             gtgtGGCATTCGCGGGAAAGCCCGCTCCCACAT‑GGTTTCATGCAACCTCTGTGGTAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATGAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAgg                                      <  2:199074/149‑1 (MQ=255)                                                                                                              tgtgGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCtt                                                                                         >  1:300669/1‑98 (MQ=255)                                                                                                              tgtgGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCtt                                                                                         <  2:300669/98‑1 (MQ=255)                                                                                                                        gcgGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaaa                           >  7:249436/1‑149 (MQ=255)                                                                                                                        gcgGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaaa                           >  5:312172/1‑149 (MQ=255)                                                                                                                        gcgGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaa                            >  4:232080/1‑148 (MQ=255)                                                                                                                          gggCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAAct                         >  8:73123/1‑149 (MQ=255)                                                                                                                            ccatgCCCCCTCCCACAG‑GGTTTATTGCACCCTCTGTGGTCGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAActct                       <  8:304581/145‑1 (MQ=255)                                                                                                                             cAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCtt                      <  6:279229/149‑1 (MQ=255)                                                                                                                             cAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCtt                      <  2:326517/149‑1 (MQ=255)                                                                                                                              aaGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                     >  1:150586/1‑149 (MQ=255)                                                                                                                              aaGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                     <  3:232080/149‑1 (MQ=255)                                                                                                                                   cGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCCACACCGCAAATTCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa                  >  6:288989/1‑148 (MQ=17)                                                                                                                                     cTCCCACAG‑GGTTTCATGCACCCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCTGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAAt              >  8:4782/1‑149 (MQ=255)                                                                                                                                         cacagaggTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGcc                              <  5:136318/123‑1 (MQ=37)                                                                                                                                         cacagaggTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGcc                              >  6:136318/7‑129 (MQ=37)                                                                                                                                         cacaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCAAGCAGCTCATCTAGCAGTAAAGAGGATGTTATCAAGCGTGAAGACAAACTCTTCGAACaaagaa           >  1:303882/1‑145 (MQ=255)                                                                                                                                          acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg          >  3:338832/1‑148 (MQ=255)                                                                                                                                          acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg          >  8:74202/1‑148 (MQ=255)                                                                                                                                                   tCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgc  <  7:254163/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATG‑AAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTC‑GGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGC  >  NC_002947/5655486‑5655778

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑