Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	2,087,676	G→C	14.3%	intergenic (+489/+429)	PP_1865 → / ← aroH	transposase/Phospho‑2‑dehydro‑3‑deoxyheptonate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	2,087,676	0	G	C	14.3%	165.4 / 22.1	56	intergenic (+489/+429)	PP_1865/aroH	transposase/Phospho‑2‑dehydro‑3‑deoxyheptonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (28/20);  new base C (4/4);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

ATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTT  >  AE015451/2087478‑2087906                                                                                                                                                                                                       |                                                                                                                                                                                                                                         aTCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGGTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGTGGAGGGAGGTGACGGGcct                                                                                                                                                                                                                  <  1:643509/223‑3 (MQ=255)     gccggcACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTACTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAgcg                                                                                                                                                                                                  <  1:562512/235‑1 (MQ=255)          cACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCggagggag                                                                                                                                                                                                                             <  1:292286/203‑1 (MQ=255)                                  tCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGAt                                                                                                                                                                                                              >  2:303665/1‑194 (MQ=255)                                     gCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTg‑gg                                                                                                                                                                                            >  2:1207271/1‑208 (MQ=255)                                            cATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCgg                                                                                                                                                         <  1:1179186/235‑1 (MQ=255)                                              tCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATgg                                                                                                                                                                                                                <  2:684883/180‑1 (MQ=255)                                              tCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATgg                                                                                                                                                                                                                >  1:684886/1‑180 (MQ=255)                                              tCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATg                                                                                                                                                                                                                 >  1:1053373/1‑179 (MQ=255)                                              tCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGGCTGCAgcg                                                                                                                                                                     >  1:849487/1‑221 (MQ=255)                                                 ggATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAgcg                                                                                                                                                                                                  <  1:28787/191‑1 (MQ=255)                                                        gATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa                                                                                                                                            >  1:541889/1‑236 (MQ=255)                                                           atCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAg                                                                                                                                           <  2:541886/234‑1 (MQ=255)                                                             cGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGTGCTTGCGGAGGGAGGTGa                                                                                                                                                                                                                         >  1:132985/1‑156 (MQ=255)                                                                cGGCAGGGCTGCCCTCGGCCTGGTGTACCGGGGAGAACGCCGTTGCTCTGGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTg                                                                                                                                                                                                                                      <  2:1147212/140‑1 (MQ=255)                                                                     tGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCt                                                                                                                                                                           >  2:184781/1‑192 (MQ=255)                                                                         tGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGag                                                                                                                           <  1:184781/236‑1 (MQ=255)                                                                           cccTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGAGGGAGATGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGcgcgtgctccc                                                                                                                                                            <  2:821251/203‑8 (MQ=14)                                                                           cccTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGAGGGAGATGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGcgcgtgctccc                                                                                                                                                            >  1:821254/1‑196 (MQ=14)                                                                             atCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTATTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATgg                                                                                                                                                                                                                <  1:344861/148‑1 (MQ=255)                                                                              tCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGccc                                                                                                                                                                                         >  1:1044378/1‑170 (MQ=255)                                                                                ggCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa                                                                                                                                            >  1:1188950/1‑212 (MQ=255)                                                                                   cTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGGGTGGAGGCATTGCGGAGGGAGATGACGGGCATGGATGccc                                                                                                                                                                                                          >  1:989333/1‑149 (MQ=17)                                                                                         tACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTg                                                                                                                                                                                               >  2:927381/1‑154 (MQ=255)                                                                                              ggAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGGGAGGGGACCCCGGAGCGTAGCGTAggg                                                                                                               >  2:888717/1‑227 (MQ=255)                                                                                                 agaaCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGccc                                                                                                                                                                                         >  2:408549/1‑151 (MQ=255)                                                                                                   aaCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATgg                                                                                                                                                                                                                >  1:132198/1‑127 (MQ=255)                                                                                                   aaCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                                                                                                   <  1:1018303/140‑1 (MQ=255)                                                                                                   aaCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                                                                                                   >  2:1018300/1‑140 (MQ=255)                                                                                                    aCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGcccaa                                                                                                                                                                                       >  2:942536/1‑148 (MQ=255)                                                                                                       ccGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCgg                                                                                                                                                                                                                                   <  1:788218/104‑1 (MQ=17)                                                                                                        cGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgt                                                                                                                        >  2:700757/1‑208 (MQ=40)                                                                                                        cGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTCGCGTAGGGGCCGGATGATGGGAGcc                                                                                              >  2:974440/1‑234 (MQ=255)                                                                                                        cGTTGCTCTTGCTCTTGCGCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGcgcg                                                                                                                                                                   <  1:976436/165‑1 (MQ=255)                                                                                                         gTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGAc                                                                                           <  1:555439/236‑1 (MQ=40)                                                                                                          ttgctcttgctcttgctcTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAg                                                                                                                                                                                                    <  1:632401/132‑1 (MQ=255)                                                                                                           tgctcttgctcttgctctTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGTCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGACTGCAGCGCGGTCCTccc                                                                                                                                                           >  1:97878/1‑170 (MQ=255)                                                                                                            gctcttgctcttgctcttGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATgggg                                                                                                 >  2:555436/1‑226 (MQ=39)                                                                                                               ccgaattcgcttctcgctTTTGATCTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGtttt                                                                                     >  1:882625/12‑236 (MQ=11)                                                                                                                cgaattcgcttctcgctTTTGATCTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGa                                                                                                 <  2:882622/213‑1 (MQ=11)                                                                                                                           cttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGTGCGTGCggagc                                                                                                                                                                                                                                >  2:800481/1‑86 (MQ=255)                                                                                                                                     cTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑g                                                                                                                                                                                     <  2:1277322/118‑1 (MQ=32)                                                                                                                                     cTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑g                                                                                                                                                                                     >  1:1277325/1‑118 (MQ=32)                                                                                                                                           cTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAggg                                                                                                                                     >  1:614367/1‑160 (MQ=31)                                                                                                                                             tCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgtagc                                                                                                                     <  1:303666/174‑1 (MQ=31)                                                                                                                                                                     aCGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGCTTTGGTTACt                                                                              >  1:865315/1‑189 (MQ=14)                                                                                                                                                                       gCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCTGGATGATGGGa                                                                                                 >  2:711190/1‑168 (MQ=11)                                                                                                                                                                       gCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCTGGATGATGGGa                                                                                                 <  1:711193/168‑1 (MQ=11)                                                                                                                                                                        cAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGAGGGAGATGACGGGCATCCATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                   >  1:574159/1‑229 (MQ=255)                                                                                                                                                                         aGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGcc                                                                                              <  2:1153418/169‑1 (MQ=32)                                                                                                                                                                         aGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGcc                                                                                              >  1:1153421/1‑169 (MQ=31)                                                                                                                                                                           aTTGCGCCTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGa                                   <  1:960227/226‑1 (MQ=17)                                                                                                                                                                              gCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGAGGGAGATGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGtttt                                                                                     <  1:1147215/173‑1 (MQ=255)                                                                                                                                                                                     cAGGAGGCCGGGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt                         <  1:952890/226‑1 (MQ=18)                                                                                                                                                                                                gcgcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGa                                                  >  2:768147/1‑190 (MQ=33)                                                                                                                                                                                                  gcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACtt  >  1:1205201/1‑236 (MQ=33)                                                                                                                                                                                                       |                                                                                                                                                                                                                                         ATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTACTTTGGCCACGACCAAAGTAACCCGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTACATCGTGAATGAATGTCTACTT  >  AE015451/2087478‑2087906

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