Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	405,538	(G)5→4	coding (107/810 nt)	proC ←	pyrroline‑5‑carboxylate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	405,538	0	G	.	89.8%	174.0 / 9.7	49	coding (107/810 nt)	proC	pyrroline‑5‑carboxylate reductase
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base . (27/17);  total (32/17)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATGCCGGAAGGCAAAAGCGCGACACAATAGAGGATTACCCAACAAAGGATGACTTTATGACAATTTGGGTG  >  NC_000913/405303‑405778                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 tttCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTg                                                                                                                                                                                                                                  >  2:194045/1‑251 (MQ=255)     cGGTCATGGACCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCt                                                                                                                                                                                                                              >  2:361884/1‑251 (MQ=255)             gCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGCAGCACTAGCAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCAc                                                                                                                                                                                                                      >  2:565464/1‑251 (MQ=255)                    gcgcgGGCAAGCTGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCg                                                                                                                                                                                                               <  1:46131/251‑1 (MQ=255)                                tGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATc                                                                                                                                                                                                      >  1:572680/1‑248 (MQ=255)                                        cGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGAccgccg                                                                                                                                                                                             >  1:30042/1‑249 (MQ=255)                                        cGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGG‑CGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGAccgccg                                                                                                                                                                                             >  2:8166/1‑249 (MQ=255)                                                   gagcaATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCAGGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGAAATCAGACCGCCGAGAATGGCTTTTc                                                                                                                                                                                <  1:79169/250‑1 (MQ=255)                                                   aacaaaTAGAAACGACCAGAGAGTCGTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTGAATCATGATGCCAGTTTTAACGGCAGCAAAAAGGAGGTCGGCGATTTGCGCCACTGCTGGCGCCGATTCTGCGGCGTGGATGCCGAACGGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCAGGCAATCAGACCGCCGAGAATGGCTTTTc                                                                                                                                                                                <  2:278603/247‑1 (MQ=255)                                                                gACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGc                                                                                                                                                                    >  1:456270/1‑250 (MQ=255)                                                                caCCAGAGAAACTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGGGCCACTTCTTGCGCCGATTCTGCGGCGTTGCTGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGG‑GTATACCCAGATTTGCCCTGGAAGCCCCTGACCGCTGGCAATCAGACCGACGAGAATGGCTGTTACCATATTGCCGc                                                                                                                                                                    >  2:443013/2‑250 (MQ=255)                                                                   cAGAGAGTCTTTATTCAGGCTAGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCTCTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGcc                                                                                                                                                                <  2:16989/251‑1 (MQ=255)                                                                      agagTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAAt                                                                                                                                                             <  1:533760/251‑1 (MQ=255)                                                                       gagTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAAt                                                                                                                                                             >  2:16285/1‑250 (MQ=255)                                                                        agTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCAAATa                                                                                                                                                            >  2:174660/1‑250 (MQ=255)                                                                             ttATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAAc                                                                                                                                                        <  1:35441/249‑1 (MQ=255)                                                                                atCAGGCTGGATGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGCTTTGCCCTGGAAGACCCTGACCGCTGGGAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGAtt                                                                                                                                                   >  2:347375/2‑251 (MQ=255)                                                                                  cAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGAtttt                                                                                                                                                 >  1:598075/1‑251 (MQ=255)                                                                                               aTTTCGCTAAGCACTTTAATCATGATGCCAGGTTTACCGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGCACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TATATACCCAGATTTGCCCTGGACGCACCTGACCGCTGGCAATCAGACCGCCGCGAATGGCTTTTCCCATATTGCCGCAGCCAATCAAACCGATTTTCTTTTCCATTGcc                                                                                                                                    >  1:252534/1‑251 (MQ=255)                                                                                                  tCGCTAAGCACTTTAATCATGATTCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCCCTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTc                                                                                                                                  <  1:203372/250‑1 (MQ=255)                                                                                                    gCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACt                                                                                                                               >  1:120128/1‑251 (MQ=255)                                                                                                      tAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGAGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTcc                                                                                                                             >  1:294418/1‑251 (MQ=255)                                                                                                      aaaGCAATTTATTCATGATGCCAGGTTTAACGTCAGCAAAAATGATGTCGGAGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTTCCCTGGCAGCACCTGACCGCTAGCAATCAGACCGCCGAGACTGGCTTTTCCCATATTGCCGAAGCCACTAAAACCGATTTTCTTTTCCCTTGCCTCACTcc                                                                                                                             >  2:56205/2‑251 (MQ=255)                                                                                                       aaGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGAGGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTGGATGCCGAACTGGTCATGCAGGGCGGCGATTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCGGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCt                                                                                                                            <  2:244941/251‑1 (MQ=255)                                                                                                          cACTTTAATCATGATGACAGGTTTAACGGCAGAAAAAATGATGGCAGCGATTTCCGCCACTTCTTGCGCCGATTCTGCGGCGTTGCTGCCGAACTGGTCATGCAGGGCGGCGCCTTTATACGGG‑AAGGGGTGTATACCCAGATTTGCCCTGGAAGCAACTGACCGCTGGCAATCAGACCGCCGAGAATGGATTTTCCAATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGcc                                                                                                                         >  1:51122/1‑251 (MQ=255)                                                                                                                          ccAGGTTTACCGACAGCAACAATGATGTCGGCGATTGGCGCAATTTCTGGCGCCGATTCTGCGGTGTTGAGGCCGAACTGGTCATGCATGGCGGCTCCTTTAGCCGGGGAGGGG‑TGTATTCCCAGTTTTGCCCTGGAAGCTCCTGACCGCTGGCAATCAGCCCGCCGAGAATGGCTTTTCCCATAGTGCCGCAGCCAATCAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATCCATTGttt                                                                                                         <  1:499788/251‑1 (MQ=255)                                                                                                                            aGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTg                                                                                                       >  2:212137/1‑251 (MQ=255)                                                                                                                             ggTTTAACGGCAGCAAAAAGGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTGGATGCCGAACTGGGCATGCAGGGCGGCGACTTTATCCGGGGAGGGT‑TTTATCCCCAGATTTGCCCTGAAAGCCCCTGGCCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCTCAGCCAATAAAACCGAGTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGa                                                                                                      <  2:66218/251‑1 (MQ=255)                                                                                                                                 tAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCCCTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATaa                                                                                                   <  1:228943/250‑1 (MQ=255)                                                                                                                                    cGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGACTTTCTTTTCCATTGCCTCACTCCTGCCGTAAAATTCATTGTTTTGATAATCGc                                                                                               >  1:278596/1‑251 (MQ=255)                                                                                                                                          caAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGc                                                                                           >  1:180423/1‑249 (MQ=255)                                                                                                                                             aaaTGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGa                                                                                        >  2:134865/1‑249 (MQ=255)                                                                                                                                               atgatgTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGCTTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAGTCGCTGGCAGAAGc                                                                                     >  1:630135/1‑250 (MQ=255)                                                                                                                                                  atgTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCCTTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATa                                                                                  >  1:216412/1‑250 (MQ=255)                                                                                                                                                         cGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGTCTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAAc                                                                          <  1:298207/251‑1 (MQ=255)                                                                                                                                                          gATTTGCGCCACTTCTTGCGCCGTTTCTGCGGCATTGATGCCGAACTGGTCATGCAGCGCGGCTACTTTATCCGGGGAGGGGG‑GTATACCCAGATTTACCCTGGCCGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACt                                                                         >  1:190501/1‑251 (MQ=255)                                                                                                                                                                  acACTTATTGCGCCGCTTCTGCGGCGTTGATGCCGAACTGGTCATGCAAGGCGGCGACTTTATCCGTGG‑CGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAAACCGCCGAGAATGGCTTTTCCCCTATTGCCACAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCCGACGCATACACAGAACTaatccgg                                                                  >  2:301/2‑244 (MQ=255)                                                                                                                                                                    aCTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATGCCGGa                                                                 >  1:452912/1‑249 (MQ=255)                                                                                                                                                                        cttGCGCCGCCTCTGCGGCATTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGCGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGAttttc                                                                                                                                                >  2:427899/1‑166 (MQ=255)                                                                                                                                                                        cttGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGTTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGAttttc                                                                                                                                                <  1:427899/166‑1 (MQ=255)                                                                                                                                                                          tGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATGCCGGAAGGCaa                                                           <  1:134865/249‑1 (MQ=255)                                                                                                                                                                                     tGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATGCCGGAAGGCAAAAGCGCGAcaca                                               >  2:517037/1‑250 (MQ=255)                                                                                                                                                                                        ggcgTTGATGCCGATCTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATCCCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATGCCGGAAGGCAAAAGCGCGACACAATag                                           <  2:171042/251‑1 (MQ=255)                                                                                                                                                                                          cgTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATGCCGGAAGGCAAAAGCGCGACACAATagcg                                         >  2:418576/1‑249 (MQ=255)                                                                                                                                                                                                ggCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATCCCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATGCCGGAAGGCAAAAGCGCGACACAATAGAGGATTAc                                   <  1:194045/250‑1 (MQ=255)                                                                                                                                                                                                    gAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATGCCGGAAGGCAAAAGCGCGACACAATAGAGGATTACCca                                >  1:85156/1‑250 (MQ=255)                                                                                                                                                                                                                              tATCCGGGGAGGGG‑TGTATCCCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATGCCGGAAGGCAAAAGCGCGACACAATAGAGGATTACCCAACAAAGGATGACTTTATGACAATTTg      <  1:16285/250‑1 (MQ=255)                                                                                                                                                                                                                               aTCCGGGGCGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATACCGGAAGGCAAAAGCGCGACACAATAGAGGATTACCCAACCAAGGATGACTTTATGACAATTTggg    >  2:145219/1‑251 (MQ=255)                                                                                                                                                                                                                                 ccGGGGAGGGG‑TGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATGCCGGAAGGCAAAAGCGCGACACAATAGAGGATTACCCAACAAAGGATGACTTTATGACAATTTGGGTg  >  1:395018/1‑251 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 TTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGCAGCAATAGAAACGACCAGAGAGTCTTTATTCAGGCTGGAGGTGATTTCGCTAAGCACTTTAATCATGATGCCAGGTTTAACGGCAGCAAAAATGATGTCGGCGATTTGCGCCACTTCTTGCGCCGATTCTGCGGCGTTGATGCCGAACTGGTCATGCAGGGCGGCGACTTTATCCGGGGAGGGGGTGTATACCCAGATTTGCCCTGGAAGCACCTGACCGCTGGCAATCAGACCGCCGAGAATGGCTTTTCCCATATTGCCGCAGCCAATAAAACCGATTTTCTTTTCCATTGCCTCACTCCTGCCGTGAAATTCATTGTTTTGATAATCGCTGGCAGAAGCATAAACAGAACTATGCCGGAAGGCAAAAGCGCGACACAATAGAGGATTACCCAACAAAGGATGACTTTATGACAATTTGGGTG  >  NC_000913/405303‑405778

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