Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	1,617,782	C→G	11.3%	intergenic (‑37/+136)	PP_1417 ← / ← PP_1418	putative Tricarboxylate transport protein TctB/putative Tricarboxylate transport protein TctC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,617,782	0	C	G	11.3%	129.2 / 8.5	43	intergenic (‑37/+136)	PP_1417/PP_1418	putative Tricarboxylate transport protein TctB/putative Tricarboxylate transport protein TctC
Reads supporting (aligned to +/- strand):  ref base C (26/12);  new base G (2/3);  total (28/15)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

ACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTCGCCGGTCATGGC  >  AE015451/1617555‑1618007                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  acCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGgc                                                                                                                                                                                                                           >  1:59369/1‑236 (MQ=255)          agcagGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGTGTCGGGATCGGGCGCGAAGcg                                                                                                                                                                                                                   <  1:484705/235‑2 (MQ=255)               gCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGTGTCGGGATCGGGCGCGAAGCGGcccct                                                                                                                                                                                                             <  1:502355/236‑4 (MQ=255)                             ggCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGGGATCGGGCGCGaa                                                                                                                                                                                                                      >  2:455720/1‑213 (MQ=255)                              gCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGTGTCGGGATCGGGCGCGaa                                                                                                                                                                                                                      <  1:916587/212‑1 (MQ=255)                                     aCCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGa                                                                                                                                                                                                                       <  2:572670/204‑1 (MQ=255)                                     aCCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGa                                                                                                                                                                                                                       >  1:572673/1‑204 (MQ=255)                                             cGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATg                                                                                                                                                                                                <  2:620744/219‑1 (MQ=255)                                                      gcgcgGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGc                                                                                                                                                                                               >  2:1058948/1‑211 (MQ=255)                                                                         tAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCCCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCCCGAACGGGCGCGAAGCTGCCGCAAAACCGGCCGac                                                                                                                                                                                                 >  2:66614/1‑189 (MQ=255)                                                                                   cAGCCGGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACAGCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTAt                                                                                                                                                                          <  2:23849/203‑1 (MQ=255)                                                                                            aGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTAc                                                                                                                                                                                      >  1:336679/1‑182 (MQ=255)                                                                                             gggCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAg                                                                                                                                                              >  1:283785/1‑205 (MQ=255)                                                                                              ggccaggccaTCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGc                                                                                                                              <  1:394721/236‑1 (MQ=255)                                                                                              ggccaggccaTCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGc                                                                                                                              <  1:960375/236‑1 (MQ=255)                                                                                                 caggccaTCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGcc                                                                                                                                                                                                    >  1:146221/1‑163 (MQ=255)                                                                                                 caggccaTCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCATCGTTGGGATTGCTCAAGTGAAATGCAAATCCGGTGGGAGCGGCCTTGCGTcgcg                                                                                                                                                                                                                                  >  2:625186/1‑133 (MQ=255)                                                                                                   ggccaTCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCt                                                                                                                                                                                                             >  2:513479/1‑152 (MQ=255)                                                                                                   ggccaTCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCg                                                                                                                                                                        >  2:413226/1‑189 (MQ=255)                                                                                                       aTCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGCTCggg                                                                                                                                                                                                                            >  2:198128/1‑133 (MQ=255)                                                                                                       aTCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATTGGGCGCGAAGCTGCCGCTAGACCTGCCg                                                                                                                                                                                                   >  1:1014878/1‑158 (MQ=255)                                                                                                       aTCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTg                                                                                                                                                                                         >  1:951841/1‑168 (MQ=255)                                                                                                                                         ccAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAggg                                                                                                                                                            >  1:858742/1‑163 (MQ=255)                                                                                                                                         ccAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGGGc                                                                                                                                                                 >  2:260109/1‑158 (MQ=255)                                                                                                                                            aGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGa                                                                                                                                            <  1:480123/176‑1 (MQ=255)                                                                                                                                                           gAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAg                                                                                                                                                                                  >  1:981431/1‑123 (MQ=255)                                                                                                                                                           gAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGc                                                                 >  1:1012163/1‑236 (MQ=255)                                                                                                                                                                               tCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGc                                                                                                                              >  2:864354/1‑155 (MQ=255)                                                                                                                                                                               tCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGc                                                                                                                              <  1:864357/155‑1 (MQ=255)                                                                                                                                                                               tCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATGCCTTGGCCAgcgc                                                           >  1:1077257/1‑222 (MQ=255)                                                                                                                                                                                cGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCAc                                            >  1:593828/1‑236 (MQ=255)                                                                                                                                                                                 gTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCGGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGGGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAgc                                                             <  1:533444/218‑1 (MQ=255)                                                                                                                                                                                  ttGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTccc                                                                                                                       >  1:51888/1‑159 (MQ=255)                                                                                                                                                                                      gATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGctt                                      >  1:611496/1‑236 (MQ=255)                                                                                                                                                                                            tCAAGTGAAAGGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCGGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGc                                                                              <  2:659722/190‑1 (MQ=255)                                                                                                                                                                                             cAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGGGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCcaca                                                                                                                    <  2:745226/151‑1 (MQ=255)                                                                                                                                                                                                        cAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATccc                                                                                                                                        >  2:52277/1‑120 (MQ=255)                                                                                                                                                                                                                   ggAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTCg             >  2:975634/1‑232 (MQ=255)                                                                                                                                                                                                                    gAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACCCAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGcc                                                          >  2:962974/1‑186 (MQ=255)                                                                                                                                                                                                                      gCGGCCTTGTGTCGGGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGAACTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGGTCCCGACACCAGGCCGCCTCCTCAGGCgtggg                                                                                                            >  2:603435/1‑133 (MQ=255)                                                                                                                                                                                                                      gCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTCGCCgg         >  2:691340/1‑233 (MQ=255)                                                                                                                                                                                                                      gCGGCCTTGCCTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAgtg                                                                                                              <  2:358521/132‑1 (MQ=255)                                                                                                                                                                                                                          ccTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCCCTTCGCCGGTCATGGc  >  1:1026279/1‑236 (MQ=255)                                                                                                                                                                                                                                gTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGc                                                                                                                                                                                               <  2:194699/41‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  ACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTCGCCGGTCATGGC  >  AE015451/1617555‑1618007

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