Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	4,907,927	Δ1 bp	9.6%	coding (98/198 nt)	PP_4317 ←	conserved protein of unknown function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	4,907,927	0	A	.	9.6%	136.0 / 7.8	42	coding (98/198 nt)	PP_4317	conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base A (25/11);  new base . (2/2);  total (28/14)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

GGCGAAGGGGCACGGATCAGGTGGTAACCGGCCTGCTCCAGCATGGGCAGGTAATCGTCGACGGTTTCCACCAGTACCAGGACTGTCTTGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTT  >  AE015451/4907694‑4908154                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    ggcgAAGGGGCACGGATCAGGTGGTAACCGGCCTGCTCCAGCATGGGCAGGTAATCGTCGACGGTTTCCACCAGTACCAGGACTGTCTTGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTGGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGttgatc                                                                                                                                                                                                                                   >  1:158554/1‑236 (MQ=255)                  cAGGTGGTAACCGGCCTGCTCCAGCATGGGCAGGTAATCGTCGACGGTTTCCACCAGTACCAGGACTGTCTTGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGttgatcttgatcttgatctt                                                                                                                                                                                                                     >  2:305094/1‑233 (MQ=255)                      tggtAACCGGCCTGCTCCAGCATGGGCAGGTAATCGTCGACGGTTTCCACCAGTACCAGGACTGTCTTGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTATTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGttgatcttgatcttgatcttgatctt                                                                                                                                                                                                               >  1:126873/1‑235 (MQ=255)                                                            gACGGTTTCCACCAGTACCAGGACTGTCTTGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGttgatcttgatcttgatcttgatcttgatcttgatcttgatcttgatc                                                                                                                                                                                         >  2:786805/1‑219 (MQ=255)                                                                gTTTCCACCAGTACCAGGACTGTCTTGCTCATGCCACCTCCGGGTGAATTCGAGGCGGCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCG‑‑‑‑‑‑TTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAg                                                                                                                                                              >  1:224008/1‑236 (MQ=255)                                                                              cAGGACTGTCTTGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTcg                                                                                                                                                                                     >  2:652469/1‑205 (MQ=255)                                                                                    tGTCTTGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCa                                                                                                                                                >  1:415165/1‑236 (MQ=255)                                                                                     gTCTTGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAg                                                                                                                                               >  1:701305/1‑236 (MQ=255)                                                                                        ttGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGCTCTTGATCTTGATCTTGCTCTTGATCTTCGCGACTTCag                                                                                                                                                                            >  2:38674/1‑204 (MQ=255)                                                                                        ttGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAg                                                                                                                                                              >  1:278734/1‑218 (MQ=255)                                                                                        ttGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGt                                                                                                                          >  1:277075/1‑236 (MQ=21)                                                                                         tGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGTCGCTTCCGTTGCCGttgttcttgatcttgatcttg                                                                                                                                                                                                                    >  2:192838/1‑163 (MQ=255)                                                                                         tGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAgg                                                                                                                                              >  1:189369/1‑233 (MQ=255)                                                                                              aTGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGa                                                                                                                                           >  2:368973/1‑231 (MQ=255)                                                                                                gCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                     <  2:277073/233‑1 (MQ=21)                                                                                                   aCCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGCTCTTGCTCTTGATCTTGATCTTGATCTTGATCTTc                                                                                                                                                                                      >  2:149634/1‑183 (MQ=255)                                                                                                               aTTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                     <  1:221595/236‑1 (MQ=255)                                                                                                                  cGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                     <  2:415163/233‑1 (MQ=255)                                                                                                                              aggCGGAGTCACCGGCCAGGACAGTCAATGGAAATAACAAGGCAAATCTGTACTGTGCTGTTGTTGCTATTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACggg                                                                                                                                       <  2:126871/201‑1 (MQ=255)                                                                                                                                       cACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCg                                                                                                                                                      >  2:335718/1‑179 (MQ=255)                                                                                                                                                    aGTCAATGGAAATAACAAGGCAAAACTTTACTGTGCTGTTGTTGCTGTTGCTGTTTCAGTTGCCGTTGCCGTTTCCGCTGCCGttgatcttgatcttgatcttgatcttgatcttggt                                                                                                                                                                                                      >  2:118477/1‑116 (MQ=21)                                                                                                                                                         aTGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGa                                                                                                                                                                                  >  2:8556/1‑133 (MQ=255)                                                                                                                                                         aTGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGatgcccgtca                                                                                                                        <  2:550704‑M2/191‑11 (MQ=255)                                                                                                                                                          tGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGc                                                                                                >  1:13169/1‑214 (MQ=255)                                                                                                                                                           ggAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCAt                                                                                                                                    >  1:690431/1‑177 (MQ=255)                                                                                                                                                                 aaCAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCG‑‑‑‑‑‑TTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCat                                                                                                        <  2:224006/193‑2 (MQ=17)                                                                                                                                                                  aCAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGatgcccgtcacctcccttcgcaagccctgcgctcggccttctgaagtcgcaatctgcgtcgcct                                                                  >  1:686579‑M2/1‑172 (MQ=255)                                                                                                                                                                         aaaaCTGTACTGTGCTGTTGTTGCTGTTGCGGTTGCCGTTGCCGTTGACGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCag                                                                                                                                                                            >  2:92165/1‑123 (MQ=255)                                                                                                                                                                          aaaCTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGG‑CCCAGGATGGGGccc                                                                                              >  1:94001/1‑198 (MQ=255)                                                                                                                                                                              tgtactgtGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCt                                                                                                                  <  1:652471/176‑1 (MQ=255)                                                                                                                                                                              tgtactgtGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGccc                                                                                                           >  1:661790/1‑183 (MQ=255)                                                                                                                                                                                   tgtgCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCg                                                                                           <  2:383278/194‑1 (MQ=255)                                                                                                                                                                                   tgtgCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGAAGGGCATGGATGCCCGTCAAgcg                                                                                                                    >  2:255016/1‑169 (MQ=255)                                                                                                                                                                                       cTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGAGCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCACGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGag                                              <  1:442799/235‑1 (MQ=255)                                                                                                                                                                                        tgttgtTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGGCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCtca                                                                                                                >  2:369336/1‑166 (MQ=255)                                                                                                                                                                                        tgttgtTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCg                                                                 >  2:280237/1‑215 (MQ=255)                                                                                                                                                                                        tgttgtTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCg                                                                 <  1:280239/215‑1 (MQ=255)                                                                                                                                                                                          ttgttgCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGTAGGTGACGGGCATCCatgcccgtcacctccct                                                                                                                 <  2:603995‑M2/165‑18 (MQ=255)                                                                                                                                                                                             ttgctgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCggt                                                                                                                                                    >  1:51776/1‑126 (MQ=255)                                                                                                                                                                                             ttgctgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGagcgcagc                                        >  1:446640/1‑235 (MQ=255)                                                                                                                                                                                              tgttgttgctgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGccc                                                                                                                            >  2:804790/4‑150 (MQ=14)                                                                                                                                                                                               gctgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAgcgcagcgca                                     >  1:755405/1‑236 (MQ=40)                                                                                                                                                                                                  gttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGttgatcttgatcttgatcttgatctt                                                                                                                                                                                                               >  2:48710/1‑63 (MQ=35)                                                                                                                                                                                                  gttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAg                                                                                                        <  2:436346/166‑1 (MQ=39)                                                                                                                                                                                                  gttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAg                                                                                                        >  1:436348/1‑166 (MQ=40)                                                                                                                                                                                                  gttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAgg                                   <  1:305096/235‑1 (MQ=37)                                                                                                                                                                                                    tgctgttgctgtTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGc                                <  1:298712/236‑1 (MQ=255)                                                                                                                                                                                                                          gttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATgggg                                                                                                 >  2:270857/1‑149 (MQ=37)                                                                                                                                                                                                                             gatcttgatcttgatcttgatctttatcttgatcttgatcttgatcttgatcttgatcttcgcgacttCAGGAGGCCGAACGCAGGCCTTGCGAAGGCAGGTGACGGGCATGGatgcccgtcacctcccttcgcaagccctacgctcggcctcctgaagtcgcaatctgcgtcgcctgaactatcgcgcgcttagaagcaaca                                        >  1:476245‑M2/11‑113 (MQ=255)                                                                                                                                                                                                                                     cgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGtt  >  2:277610/1‑233 (MQ=32)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    GGCGAAGGGGCACGGATCAGGTGGTAACCGGCCTGCTCCAGCATGGGCAGGTAATCGTCGACGGTTTCCACCAGTACCAGGACTGTCTTGCTCATGCCACCTCCGGGTGAATTCGAGGCGTCATTATGCGGAGTCACCGGCCAGTACAGTCAATGGAAATAACAAGGCAAAACTGTACTGTGCTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTT  >  AE015451/4907694‑4908154

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