Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	BA000036	3,208,522	+T	coding (243/1179 nt)	Cgl3007 →	Malic enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BA000036	3,208,520	1	.	T	98.3%	191.7 / ‑3.5	58	coding (241/1179 nt)	Cgl3007	Malic enzyme
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (28/29);  total (29/29)
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.41e-01

TGAGGTTCTACCCTTAACAGAGCTTCCCGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑G‑TTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAGCGTCTCGATATTCCAGTTATGCACGATGACCAGCACGGCACCGCTGTGGTTATCCTCGCTGCGC  >  BA000036/3208221‑3208814                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                       tGAGGTTCTACCCTTAACAGAGCTTCCCGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTAACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAATCAATCAAGGAAGATCCAGAGGCTGCGCGCACCCACACGGGCATTGGAAACATCGTCGCGGTTATTTCCGACGGCACCGCT‑g                                                                                                                                                                                                                                                                                                         >  1:257181/1‑300 (MQ=255)           cccTTAACAGAGCTTCCCGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTcttggc                                                                                                                                                                                                                                                                                                <  1:213723/299‑1 (MQ=255)            ccTTAACAGAGCTTCCCGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCCCCCGGGCATTGGCAACACCGTCGCGGTTATCTCCGACGGCACACCT‑GTTTcttggcct                                                                                                                                                                                                                                                                                              >  2:534666/1‑300 (MQ=255)               tAACAGAGCTTCCCGCAAAAACACCGATTAACAAGGCTAAATGATATGTCCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTcttggccttgg                                                                                                                                                                                                                                                                                           >  2:450726/1‑300 (MQ=255)                aaCAGAGCTTCCCGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTcttggccttggc                                                                                                                                                                                                                                                                                          >  1:382696/1‑300 (MQ=255)                aaCAGAGCTTCCCGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATGTCCGACGGCACCGCT‑GTTTcttggccttg                                                                                                                                                                                                                                                                                            <  1:260468/298‑1 (MQ=255)                  cAGAGCTTCCCGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCACGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTTGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGa                                                                                                                                                                                                                                                                                        >  1:31916/1‑300 (MQ=255)                     agCTTCCCGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCTGG‑TTCTTGGCCTTGGCGatat                                                                                                                                                                                                                                                                                     >  2:286936/1‑300 (MQ=255)                     agCTTCCCGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGatat                                                                                                                                                                                                                                                                                     <  1:143802/300‑1 (MQ=255)                          cccGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCgg                                                                                                                                                                                                                                                                                  <  2:110497/298‑1 (MQ=255)                          cccGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGAc                                                                                                                                                                                                                                                                                >  2:58513/1‑300 (MQ=255)                                              cAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAAGCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCACTTCCCGTc                                                                                                                                                                                                                                                            <  1:277626/300‑1 (MQ=255)                                                   cTACATGATATTACCCTCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGTCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTCTTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATCTCGGCACTTCGGCCTCCTTTCCCGTCCTggg                                                                                                                                                                                                                                                       >  2:43438/1‑299 (MQ=255)                                                      aagaatcttgccATCTACCTTCAGCCTTCCACCCAAAAACCCACCCAGAATGCATCCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTATTTCCACCCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTGGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAggg                                                                                                                                                                                                                                                    <  2:189414/290‑1 (MQ=255)                                                                  aTCGACCTGCATCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGCAAGCTCTCCATTATTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCt                                                                                                                                                                                                                                        >  1:299706/1‑300 (MQ=255)                                                                  aTCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCt                                                                                                                                                                                                                                        <  1:58513/300‑1 (MQ=255)                                                                   tCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGGCATGGAGGGCAAGGCTCAGCTg                                                                                                                                                                                                                                       >  1:408070/1‑300 (MQ=255)                                                                                ttCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGGTTAGC‑CTTTccg                                                                                                                                                                                                                         >  2:125786/1‑298 (MQ=255)                                                                                   cACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCTTGGAGCGCAAGGCTCAGCTGTTTAGCTCTTTCGCTg                                                                                                                                                                                                                       >  2:169715/1‑300 (MQ=255)                                                                                        aaaaCCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCt                                                                                                                                                                                                                   >  1:404336/1‑299 (MQ=255)                                                                                                     tGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTa                                                                                                                                                                                                     <  2:397714/300‑1 (MQ=255)                                                                                                     tGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTa                                                                                                                                                                                                     >  1:50626/1‑300 (MQ=255)                                                                                                          aaaTCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGtt                                                                                                                                                                                                <  1:261275/300‑1 (MQ=255)                                                                                                            accgtCGAGGCACACTGGGCCTGAAAGCTCTCCATTCGTTCCACTCCTCCTCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGtttt                                                                                                                                                                                              <  2:506330/296‑1 (MQ=255)                                                                                                                            aGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGtt                                                                                                                                                                              <  1:25587/300‑1 (MQ=255)                                                                                                                               gCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGAc                                                                                                                                                                           >  1:342094/1‑300 (MQ=255)                                                                                                                                cGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACg                                                                                                                                                                          <  1:193635/300‑1 (MQ=255)                                                                                                                                            cATCAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGGGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTgaga                                                                                                                                                              <  2:382696/300‑1 (MQ=255)                                                                                                                                               tAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTATTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCa                                                                                                                                                           >  1:271404/1‑300 (MQ=255)                                                                                                                                                    ccACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTGTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTGCACGATGTTGACGCTTTGGTTGAGACCATCGc                                                                                                                                                       <  2:135880/299‑1 (MQ=255)                                                                                                                                                     cACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCAAGATGTTGACGCTTTGGGTGAGACCATCGCAg                                                                                                                                                     >  1:189548/1‑300 (MQ=255)                                                                                                                                                                tccgcGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTCTCGTTTTGGACGTTCACGATTTTGTCGCTTTGGTTGAGACCATCGCACCCATCgcgcat                                                                                                                                          >  2:20553/1‑298 (MQ=255)                                                                                                                                                                 ccgcGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCctt                                                                                                                                         <  1:450726/300‑1 (MQ=255)                                                                                                                                                                  cgcgACGTGCGCATTCTTGCCCTGGCTTGCACCCCTGGTTTTGCTCGGGTTTGTGAAGCATTCACGGAGGTTCCAGCGGTTGCGCGCACCCACCCGGGCATGGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCcttc                                                                                                                                        <  2:189548/300‑1 (MQ=255)                                                                                                                                                                     gACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCttt                                                                                                                                     <  1:443602/300‑1 (MQ=255)                                                                                                                                                                                   ttCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCCGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTCCACGATGTTGACGCTTTGGTTGAGACCATCTCAGCCATCGCGCCTTCTTTCGTTGCTATCAACt                                                                                                                       >  1:521697/1‑300 (MQ=255)                                                                                                                                                                                     cccTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGTCGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTGTCGGTGCTATCAACtt                                                                                                                      <  1:193871/299‑1 (MQ=255)                                                                                                                                                                                       ctggCTTACACCCCTGTTGTTGCTTAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTGGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTgga                                                                                                                   <  1:311052/297‑1 (MQ=255)                                                                                                                                                                                                acCCCTGGTGTGGCTCAAGTTGGGGAAGCAATCCAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTCTCGGTGCTATCAACTTGGAGGACAtctc                                                                                                          <  2:252263/300‑1 (MQ=255)                                                                                                                                                                                                         gtTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGTTTGAGACCATCGCAGCCATCGCCCCTTCTTTCgg                                                                                                                                  >  2:299328/1‑267 (MQ=255)                                                                                                                                                                                                          ttGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGt                                                                                                                                 <  1:299328/267‑1 (MQ=255)                                                                                                                                                                                                          ttGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCCACTTGGAGGACATCTCCGATCCTCGt                                                                                                >  2:535753/1‑300 (MQ=255)                                                                                                                                                                                                           tGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCCGCCATCGCGCCTTCTTTCAGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGtt                                                                                               >  2:310750/1‑300 (MQ=255)                                                                                                                                                                                                              tCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCt                                                                                            >  2:275352/1‑300 (MQ=255)                                                                                                                                                                                                                      gtgAAGCACTCAAGGAAGTTCCGGCGGTTGCGCGCCCCACCGCGTGCCTTGGACACACCGTCGCGGTTATTTCCGACGGCCCCGCT‑GTTTCTTGGCCTTGGCGATATCGGCCCTCAGGCCTCCCTTCCCGTCATGGAGTGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCGCGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGt                                                                                     <  2:519746/299‑1 (MQ=255)                                                                                                                                                                                                                      gtgAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTGTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGGTGACGCTTGGGTGGAGACCATCGCAGCCATGGCGCCTTCTTTCGGGGATATCAACAGGGAGGACCACACAGGTCACAGTTGCTTCGGGGTg                                                                                    >  2:234917/1‑300 (MQ=255)                                                                                                                                                                                                                       tgAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGTgg                                                                                   <  1:434928/300‑1 (MQ=255)                                                                                                                                                                                                                                        agccAGAGGTTGCGAGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAg                                                                  <  2:137539/298‑1 (MQ=255)                                                                                                                                                                                                                                            agagGTTGCGCGCACCCACACGGGCAGTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAGCGtc                                                              <  1:123580/300‑1 (MQ=255)                                                                                                                                                                                                                                                                    cATTAGAAACACGGTCGCGGTTATTGCCGCCGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGGTGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAGCGTCTCGATATTCCAGTTATGCACGATg                                      <  2:408070/300‑1 (MQ=255)                                                                                                                                                                                                                                                                     aTTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGATGACATCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAGCGTCTCGATATTCCAGTTATGCACGATGa                                     >  1:254013/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                 gTCGCGGTTATTTCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGCAGGCCATCCCTATCGCTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCCCCGCTCCTCGTTCCTTCCAGGCGGAGCGTCGACTCCTCGCGCGCCCCCATATTCCTGCTATGCACGATGAACAGCACGGCAc                         >  2:160740/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                     cgGTTATTTCCGCCGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAGCGTCTCGATATTCCAGTTATGCACGATGACCAGCACGGCACCGCt                     <  2:404336/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                       gTTATTACCGACGGCTCCGCG‑GTTTCTTGGCCGTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAGCGTCTCGATATTCCAGTTATGCACGATGACCAGCACGGCACCGCtgt                   <  2:367571/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                            ttCCGACGGCACCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCCACTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAGCGTCTCGATATTCCAGTTATGCACGATGACCAGCCCGGCACCGCCGTGGCTa              >  2:196340/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                   ggCCCCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCGCCCTTCCCGTCAGGGAGGGCCAGGCTCTGCTGTTTCGCTCTTTCGCTGGCCTGACGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAGCGTCTCGATATTCCAGTTATGCACGATGACCAGCACGGCACCGCTGTGGTTATCCTcgc       <  2:13991/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                      aCCGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAGCGTCTCGATATTCCAGTTATGCACGATGACCAGCACGGCACCGCTGTGGTTATCCTcgctgc    >  1:514820/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                        cGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCCGCTGTTTAGCTCTTTCGCTGGCCTGAAGTCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACCTCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCCAATCGCGCGTCTCGATATTCCAGTTATGCACGATGACCAGCCCGGCACCGCTGTGGTTATCCTcgctgcgc  >  2:388183/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                        cGCT‑GTTTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAGCGTCTCGATATTCCAGTTATGCACGATGACCAGCACGGCACCGCTGTGGTTATCCTcgctgcgc  >  2:81107/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                       TGAGGTTCTACCCTTAACAGAGCTTCCCGCAAAAACACCGATTAACAAGGCTAAATGATATGACCATCGACCTGCAGCGTTCCACCCAAAACCTCACCCATGAGGAAATCTTCGAGGCACACGAGGGCGGAAAGCTCTCCATTAGTTCCACTCGTCCGCTCCGCGACATGCGCGATCTTTCCCTTGCTTACACCCCTGGTGTTGCTCAGGTTTGTGAAGCAATCAAGGAAGATCCAGAGGTTGCGCGCACCCACACGGGCATTGGAAACACCGTCGCGGTTATTTCCGACGGCACCGCT‑G‑TTCTTGGCCTTGGCGATATCGGACCTCAGGCCTCCCTTCCCGTCATGGAGGGCAAGGCTCAGCTGTTTAGCTCTTTCGCTGGCCTGAAGGCTATCCCTATCGTTTTGGACGTTCACGATGTTGACGCTTTGGTTGAGACCATCGCAGCCATCGCGCCTTCTTTCGGTGCTATCAACTTGGAGGACATCTCCGCTCCTCGTTGCTTCGAGGTGGAGCGCCGCCTCATCGAGCGTCTCGATATTCCAGTTATGCACGATGACCAGCACGGCACCGCTGTGGTTATCCTCGCTGCGC  >  BA000036/3208221‑3208814

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑