Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,151,329	T→G	23.8%	A590A (GCA→GCC)	ppc ←	phosphoenolpyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,151,329	0	T	G	23.8%	101.1 / 2.7	67	A590A (GCA→GCC)	ppc	phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (34/12);  new base G (11/7);  total (46/21)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTT  >  NC_000913/4151094‑4151565                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             ggTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCACTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCCTGAGGAGGtct                                                                                                                                                                                                                                            >  1:126694/1‑236 (MQ=255)     gACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTTG‑GACCACAGCGCCCCCTGCGCAGGGGggccg                                                                                                                                                                                                                                        >  2:1303247/1‑238 (MQ=255)      aCATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGCGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCTCCTTTCAGGCTTCCTGGCGGTTTT‑GAAAGGACCGCCGCATGAGCAGGGGggccgc                                                                                                                                                                                                                                       >  1:1340019/1‑238 (MQ=255)       cATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGcgccgccg                                                                                                                                                                                                                                     <  1:946423/239‑1 (MQ=255)       cATCGCAGGAGATGACTGACAGTTCATCAATAATGCGAGGCCAGCTCTCTGTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTGAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCGTTCAGGCTGCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGcgccgcc                                                                                                                                                                                                                                      <  2:58767/238‑1 (MQ=255)         tCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATCTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCcgcg                                                                                                                                                                                                                                   <  1:371563/239‑1 (MQ=255)          cGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGGGCGCCGCcgcg                                                                                                                                                                                                                                   >  1:101826/1‑238 (MQ=255)          cGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAG‑‑CTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCACGAGCCGGGGCGCCGCcgcgga                                                                                                                                                                                                                                 >  1:1054388/1‑236 (MQ=255)               gagaTGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATAATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCACCCCCGCCTGAGCCGTGGCGCCGCCGGGACCAAt                                                                                                                                                                                                                             >  1:202000/1‑239 (MQ=255)                         acAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTTG‑TACAGCAACGCCGCCTGAGCAGGTGCGCCGCAGCCGACAAAGGAAccgcc                                                                                                                                                                                                                    >  2:1228706/1‑238 (MQ=255)                             ttCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAAGCCGCCTTTCAGGCTTCCTGGCGGCTGC‑TACAACACCGCCGCACGAGAAGATGCGTCCCCCCCACCACTAGAACcgccgccg                                                                                                                                                                                                                >  2:786143/1‑236 (MQ=255)                               cATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTTT‑GACAGCAGCGCCGCATGAGCCAGGGCGCCGCCGCGACCAATGGAACCGCCGCGcccc                                                                                                                                                                                                             >  1:1237757/1‑238 (MQ=255)                                     tAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGACTTTCAGGCTTCCTGGCGGGTGT‑GACAGCCGGGGCGGATGAGCCGGTGGGCCGCCGCCGACAAAGGAACCGCCGCGAGCGGGGaaa                                                                                                                                                                                                       >  2:113904/1‑238 (MQ=255)                                        tGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGGGGaacaac                                                                                                                                                                                                    >  1:1171719/1‑239 (MQ=255)                                            aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCACCGCCGCATGCGCAGGTGCGCCGCCGCGACCACTGGACACGCCGCGGCCGCGGCAACACGTCa                                                                                                                                                                                                >  2:29693/1‑239 (MQ=255)                                                        tttCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGGATGAGCAGGTGCGCCGCCGCGACCAATGGAAACGCCCCCACCCGGGAACAACGTCCGCCCAAAAccc                                                                                                                                                                                     >  2:1378109/1‑238 (MQ=255)                                                                       gggcggtgggTTGGCGTCCAGAGTCCCCCCGTTATAAAGGGGCGGTCTGCGGACGTGTTTTTCTG‑GCAGACCATTTTTAATGGGGATCACCTCGCCCTGGTCGGTTACGCGCGGGCCGCCTTTCAGGTTTCCTGGGGGTTGT‑GACAGCGGCGCCGCATGAGCAGGTGCCCCGCCGCGCCCATTGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGtt                                                                                                                                                                      <  2:596800/230‑1 (MQ=255)                                                                               tgTTGGCTCCCAGAGCCGCCCGGCTAAAAAGCGACGGGGTGGGGTGGGTGTTTTCTG‑GCAGCCCCTATTAAAAGCTGATCATCTCGCCCTGTGCGGTTGCGCGCAGCCCCCCTTCCAGGCTTCCTGGCGGTGGC‑GACCGCAGCGCCGCAGGAGCGGGGGCGCCGCCGCGACCAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAAt                                                                                                                                                             <  2:889630/238‑1 (MQ=255)                                                                                        ccAGAAGCCCCCCGTTATAAAGCGATAGGCGGGTGAGGGTGTTTTCGG‑GCGGCCCTTATTTAAAGCGCTTCATCCCGCCTTGTTCGGTTCCGCGCAGGCCGCCTTTCAGGCTTCGGGGGGGTTGT‑GACAGCAGCGCCCCATGGGCAGGGGCGCCGCCGCCACCAATGAAGCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCtgt                                                                                                                                                    <  2:1051564/239‑1 (MQ=255)                                                                                          aGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCCGCGCCGCCTGGGCAGGGGCGCCGCCCCGACCAAAGGAACCGCCGCGGCCCTGGGACAACGCCAGCTCACTACCCCCCTTTTCCCAGGGTTTGTTTAATTCCTCCtgtg                                                                                                                                                   >  2:72028/1‑238 (MQ=255)                                                                                                           aaggcacgggcggcTGA‑GGTTATTTCGGTGCAGAACTTATTCAAAGCGCCCCATCCCGCCTTTGCCGATTCCGCGCGCGCCGCCTTTCGGGTTGCCGGGCGGTAGT‑GCGAGCAGCCCCGTATGGGCGGGTGCCCGCCCGCACCCAATGAACCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAgg                                                                                                                                 <  2:938457/231‑1 (MQ=255)                                                                                                                aCAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTCATATTGCGCCCCggca                                                                                                                               >  1:1394129/1‑234 (MQ=255)                                                                                                                 cgggCTGCGGGCGTTGTTTTCGG‑GCAGCCTTTTTTTAAGGGGTACCCTCGCCCCCTTTTGGTTTCCGCCCAGCCCGCCTTCCGGGCTTCCTGGGGTTTGT‑GACAGCGGCGCCGCATGAGCGGGTCCGCCCCCGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCt                                                                                                                            <  2:1126347/236‑1 (MQ=255)                                                                                                                 cAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTTT‑GACAGCAGCGCCGCATGAGCAGGGGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCt                                                                                                                            >  1:1334210/1‑238 (MQ=255)                                                                                                                 cAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCt                                                                                                                            >  1:194610/1‑238 (MQ=255)                                                                                                                 cAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAAACGCCGCGCACGTGGAACAACGGCAGCTCAATACCCGCTTTTTCCCAAGTTTTGATTAATACATCCTGTGCCCTGTATTTCGCCCAGGAAGCt                                                                                                                            >  2:1207071/1‑238 (MQ=255)                                                                                                                      tgctgACGGTGGTTTCTG‑GCAGACCATATTTAAAGCGTATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCCACCAATGGAACCGCCGCGCCCCTGGAACAACGTCAGCCCCATACCCCCCTTTTCGCAGGTTTTCATTAATGCATCCTGGGCCTGGTATTTCTGCCAGGAAGCGTCCAt                                                                                                                       >  2:1160170/1‑238 (MQ=255)                                                                                                                        ctgACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATc                                                                                                                      >  1:1199554/1‑237 (MQ=255)                                                                                                                               ttttttCTG‑GCCGACAATTTTCAAAGCGGCTCCCCCCGCCCTGTTGGTTTGCGCGCGCGCCCCTTTTCGGTCTTCGTGGGGGTTGT‑GCCACCGCCCCCACTGGAGGAGGGCCCCCCCCCCGACAAAGGTAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGc                                                                                                              <  2:1334210/235‑1 (MQ=255)                                                                                                                               tgttttCTG‑GCATACCATATTAAAAGCGTCCCCCCCCGCCCTGTTCGTTTCCGCGCGGGCCCCCTTTCAGCCTGCCTGGCGGTTGT‑GACAGCAGCGCCGCAGGAGCGGGGGCGCCGCCGCGACCAATGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGc                                                                                                              <  2:194610/235‑1 (MQ=255)                                                                                                                                  tttCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCCCCGCGACCGTGGGACAACGTCACCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCACGAAGCTGCCATCACTCCCGCATCt                                                                                                          >  2:1287046/1‑239 (MQ=255)                                                                                                                                       g‑gCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCtttt                                                                                                       >  1:198583/1‑237 (MQ=255)                                                                                                                                          cAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATg                                                                                                                                                            >  1:1030215/1‑182 (MQ=255)                                                                                                                                            gACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGTCGGTTGT‑GACAGCAGCGCCGCCTGAGCAGGGGCGCCGCCGCGACCAATGGAACCGCCGCGCCCGTGGAACAACGTCCGCTCAATACCCGCTTTTTCGCATGTTTTGCTTAATCCCTCCTGCTCCTGATATTTCGTCCAGGGCACTTCCATCACTCCCGCATCTTTTTCTCAg                                                                                                 >  2:327242/1‑239 (MQ=255)                                                                                                                                                            gATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGGGCCGCCGCGACCAATGGAACCGCCGCGGCCGTGGGACAACGTCAGCCCAATACCCGCTTTTTCGCAGGGATTGATTAATGCATCCTGTGCCTGATATTGCGACCAGGAAGCTGCAATCACTCCCCCATCTTTTGTCGGGTACGAATCGGCAatca                                                                                 >  2:325215/1‑239 (MQ=255)                                                                                                                                                            gATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGGACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTACCGCATCTTTTGCTGAGTCGGAATAGCCAatca                                                                                 >  2:120491/1‑239 (MQ=255)                                                                                                                                                             atcatcTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGGAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCCGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatca                                                                                 >  1:380892/1‑238 (MQ=255)                                                                                                                                                             atcatcTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCCGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCCGGAAGCTGCCATCACCCCCGCCTCTTTTGCTGCTGCGGAATACCCAatca                                                                                 >  2:1039839/1‑238 (MQ=255)                                                                                                                                                             atcatcTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatca                                                                                 >  1:784876/1‑238 (MQ=255)                                                                                                                                                             atcatcTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatca                                                                                 >  1:786479/1‑238 (MQ=255)                                                                                                                                                                     gCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCGGCCATCACTCCCGCCTCTTTTGGTGAGTCGGGATAGCCAATCATCACCa                                                                           >  2:398153/1‑236 (MQ=255)                                                                                                                                                                        cTGTTCGGTTACGCGCACGCCGCCTTTCGGGCTTCCGGGGGGTGGT‑AAGAGCGGCGCCGAAGGGGGAGGTGCGCCGCCGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGtt                                                                     <  1:261934/239‑1 (MQ=255)                                                                                                                                                                          gTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGGGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGGTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGttt                                                                    >  2:394447/1‑238 (MQ=255)                                                                                                                                                                           ttCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTg                                                                   >  2:1293107/1‑238 (MQ=255)                                                                                                                                                                           ttCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTg                                                                   >  2:253957/1‑238 (MQ=255)                                                                                                                                                                                ttACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGGTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGATGAGTCGGAATAGCCCATCATCCCCATCTGTTTGCCCTgc                                                             >  1:1014027/1‑238 (MQ=255)                                                                                                                                                                                   cgcgCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGGGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTCCTGAGTCGGAAAAGCCAATCATCACCATCTGTTTGCCCTGCATc                                                          >  2:60029/1‑239 (MQ=255)                                                                                                                                                                                    cccgcggCCGCCTCCCAGGCGTCGGGGCGGTTGC‑CACAGCGCCGCCGAGTGAGCGGGCGCGCCGCCGCACCAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCa                                                         <  2:1054388/234‑1 (MQ=255)                                                                                                                                                                                     cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCCGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGGTTGCCCTGAATCAg                                                        >  2:705671/1‑239 (MQ=255)                                                                                                                                                                                     cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTTCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAAACAg                                                        >  2:1271521/1‑239 (MQ=255)                                                                                                                                                                                     cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGcc                                                                 >  1:977951/1‑230 (MQ=255)                                                                                                                                                                                     cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGCGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCCGCTCAATACCCGCTTTTTCGCAGGGTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGACATCACTCCCGCCTCTTTTGCTGAGTCGGAATATGCCATAATCACCATATGTTTGCCCTGAAtacg                                                        >  2:749779/1‑236 (MQ=255)                                                                                                                                                                                     cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGG‑GACAGCAGCGCCGCATGAGCATGGGCGCCGCCGCGAACAATGGAACCGACGCGACCCTGGAACAACGTCAGCTCAATAACCGCTTTTTCGAAGGTTTTGATTAATGCAACCTGTGCCTGATATTGCGCCCAGGAAGATGCCATCACTCCCCCCTGTTTTGCTGCGTGAGGATAGCCAAACATCACCATCTGTTTGCCCTGAATCAg                                                        >  2:274557/1‑239 (MQ=255)                                                                                                                                                                                        gggCCCCTTTCCAGCCCTCCGGGCGTTTGG‑GACGGCAGCCCCGCGAGAGCGGGTGCCCCGCCGCGCCCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGc                                                      <  1:749779/237‑1 (MQ=255)                                                                                                                                                                                          gccgccTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCAc                                                   >  1:1075699/1‑239 (MQ=255)                                                                                                                                                                                              ccTTTCAGGCTCCCGGGCGTTTGA‑GACCGGCGCCCCGCGTGAGCAGGTGCCCCGCCGAGACCAATGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa                                               <  1:673273/239‑1 (MQ=255)                                                                                                                                                                                                ttccAGGCTCCCTGGCGGTTGT‑AACGGCGGCGCCCCATGGGAGGGGCCGCCGCCGCGCCCAAGGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATAcc                                             <  1:60029/236‑1 (MQ=255)                                                                                                                                                                                                   ccggCTTCTGGGGGTTTGT‑GACAGCACCCCCGAAGGGGCGGGCGCCCCGCCGCGCAAAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                          <  1:786143/237‑1 (MQ=255)                                                                                                                                                                                                    aGGCTCCCTGG‑GGTTGTGGACAGCGCCCCCGCAGGGGCGGGTCCCCCGCCCAGCACATGGAAACCGCCGCGACCGTGGAACAGCGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                          <  1:113904/238‑1 (MQ=255)                                                                                                                                                                                                         tCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCCGTCAatat                                    >  1:368268/1‑239 (MQ=255)                                                                                                                                                                                                         tCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCA‑CTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATAtt                                   >  1:1171110/1‑239 (MQ=255)                                                                                                                                                                                                                  ttGT‑GACAGCCGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACCTCAGCTCAATACCCGCTTTTTTGCAGGGTTTGATTAATGCATCCTGTGCCTGATATTTGGCCCAGGAAGCGGCCATCACTCCCGCATCTTTTGCTGAGTCGGACTAGCCAATCATCACCATCTGTTTGCCCTTAATCAGGCCACGATACCAGGCAATATTGcgca                              >  2:446570/1‑236 (MQ=255)                                                                                                                                                                                                                     ggaaCAGCGGCGCCACATGGGCGGGGGCCCCGCCGCGACCAATGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTgg                         <  1:1378109/236‑1 (MQ=255)                                                                                                                                                                                                                             accgccgcAGGGCCGGGGGCCCCGCGGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCAt                    <  1:1243984/234‑1 (MQ=255)                                                                                                                                                                                                                              gcgcCGCATGAGGGGGTGCGCCGCCGCCACAAAGGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGAc                 <  1:398153/238‑1 (MQ=255)                                                                                                                                                                                                                                        ggggggTGCCCCGCCGTGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCg       <  1:1207071/234‑1 (MQ=255)                                                                                                                                                                                                                                         gCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCgtt     >  1:1025865/1‑239 (MQ=255)                                                                                                                                                                                                                                            ggTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCgttgt   >  2:1409572/1‑238 (MQ=255)                                                                                                                                                                                                                                              tGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCAt                                                                                                                       >  2:758479/1‑120 (MQ=255)                                                                                                                                                                                                                                              tGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCgttgtt  <  1:1160170/237‑1 (MQ=255)                                                                                                                                                                                                                                              tGCGCCGCCGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCAt                                                                                                                       <  1:758479/120‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             GGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTG‑GCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGT‑GACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTT  >  NC_000913/4151094‑4151565

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