Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,366,216	T→G	19.8%	T187P (ACT→CCT)	dcuA ←	C4‑dicarboxylate antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,366,216	0	T	G	19.8%	125.4 / 2.3	67	T187P (ACT→CCT)	dcuA	C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base T (39/12);  new base G (4/11);  total (43/24)
Fisher's exact test for biased strand distribution p-value = 1.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGATTGCGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAG  >  NC_000913/4365985‑4366454                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                tCAGGATTGCGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCagg                                                                                                                                                                                                                                          >  2:226264/1‑239 (MQ=255)  cAGGATTGCGTTGGGGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCagga                                                                                                                                                                                                                                         <  2:49001/239‑1 (MQ=255)     gATTGCGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAAACGGATCGGCAGAAAGTTTGGAGTTG‑AACAACAAAATAGCCagga‑agg                                                                                                                                                                                                                                     >  1:274682/1‑239 (MQ=255)        tGCGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACa                                                                                                                                                                                                                                  >  2:589839/1‑239 (MQ=255)          cGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACCCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GGGTAAAATCGGATCGTCAGAAAGTTTGGAGTTG‑AACAGACTAGTGAACggga‑aggaaaac                                                                                                                                                                                                                                >  1:572885/1‑235 (MQ=255)               gtggtgTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGCTAAATCGGATCGTCAGAAAGTTTGGGGTTG‑AACAGCATAGTGACCAGGA‑AGGAAATCACCAg                                                                                                                                                                                                                           >  2:637655/1‑239 (MQ=255)               gtggtgTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATAACCAg                                                                                                                                                                                                                           >  2:552408/1‑239 (MQ=255)               gtggtgTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGAACTCACCAg                                                                                                                                                                                                                           >  2:678369/1‑239 (MQ=255)                     ttCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGGGTTG‑AACAGCATAAGGCACAGGA‑AGGACATCACCAGAACCg                                                                                                                                                                                                                      >  2:248603/1‑238 (MQ=255)                             cGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGTTAAATCGGATCGTCAGAGGGTTTGGAGTTG‑ACAAGCATAGGGACCAGGA‑AGGAAATCACCAAAACCGCAAGAAggg                                                                                                                                                                                                             >  1:83350/1‑239 (MQ=255)                                 ttttCACCCACCCCCCTGTTGGGGGTGTGGATAGTTGCAGAGATACCCCCGCCAACTCGC‑CCCAGAAGGAACGGCCAGACGGACGTTTTTCCACGGGATTTGATTTA‑AATCTG‑TTTTCCACC‑GCGC‑ATTCCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAAGCGGATCGTAAGAGAGTTGGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGa                                                                                                                                                                                                         <  1:226264/239‑1 (MQ=255)                                                    tggggtttttGATGATTGCAGAGATAACCCGCCCAAATACC‑CCCAGAGGGAACAGCCAGACGGACTTTTTGGCCCCGGTTTTGATTTC‑AATCTG‑CTTTTCCCC‑GCGC‑AGTTCACCCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGagag                                                                                                                                                                                       <  2:635899/230‑1 (MQ=255)                                                      gggCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑ATTTTCACC‑GCGC‑AGTTCAACAAGGCCC‑TCTCACAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCATCCGAACCGCCAGCAGGGAGGACGGGATGACCACGGAGAGCAg                                                                                                                                                                                    >  2:61817/1‑239 (MQ=255)                                                           gTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCCCCAAAACCGCCAGCAGGGGGGACGGGATAACAAAGGAGGGAAGAgggg                                                                                                                                                                               >  1:309217/1‑238 (MQ=255)                                                            ttGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGCGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGAAAATgggg                                                                                                                                                                              >  1:358014/1‑237 (MQ=255)                                                            ttGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAAAACCCCCAACAGGGGGGCAGGGAGTACCACGGagaggaaagggga                                                                                                                                                                              >  1:284592/1‑230 (MQ=255)                                                            ttGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGGTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAAAACACCCAGAAGGGTGAACGGGATGACCACGAAGAGGAGAggggg                                                                                                                                                                              >  1:401055/1‑237 (MQ=255)                                                                  atTGGATATATACCCACGCCGACTGCG‑CCAATCATGTACACCCGGCCGAGTCTTTGTGCACCGGATGTAATTTC‑ATTCGG‑CTTTTC‑CC‑CCGCGAGTTCACCCAGTCCG‑TCTTCCGGACGCTG‑GCAAAAAACGGGACCGTAAGAGACTTGGAATTGG‑AACGTTATAGGGCCCAGAA‑AGGACACCACCGAAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCt                                                                                                                                                                        <  1:251934/239‑1 (MQ=255)                                                                  atTGAATAGATACCCCCG‑CAACTTCGCCCCAGAAGGACAGGCCAGCCGGACTTTTTGGCCCCGG‑TTTGAATTCAAATCTG‑TTTTCCCCC‑GGGA‑AGTTAAACAAGGCCT‑TTTTCAAGACGCTT‑GCGATAAATCGGATCGTAAGAGAGTTGGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCt                                                                                                                                                                        <  2:309217/239‑1 (MQ=255)                                                                    tGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGGATGAGAGTTTGGAGTTG‑AAAAGCATAGGGACCAGGA‑AGGACATCACCAGAACCGCCAACAGGGTGGACGGGGTGACCACGGAAAGCAAATGGAGGCAgggaa                                                                                                                                                                      >  1:423069/1‑235 (MQ=255)                                                                      cATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGGAGCTGAt                                                                                                                                                                     >  2:185176/1‑238 (MQ=255)                                                                              aCCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGAACCTGATGCCAGGa                                                                                                                                                              >  2:142937/1‑237 (MQ=255)                                                                                 cccgccactctcccccAGCAGGAACCACAAGACGGACTTTTGTGCACCGTATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑ATTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATGGGATCGTCAGAGAGTTTGGAGTTA‑AACAGCATAGTGACCAGGA‑AGGACATCCCCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCtt                                                                                                                                                          <  2:466581/226‑1 (MQ=255)                                                                                 accccAACTCCC‑CCCGGCAGAAACACCCAGCCGGACGTTTTTGCACGGGATTTGATTCC‑AATCTG‑CTTTCCCCC‑GCGC‑AGTTCAACCAGGCCT‑TCTTCCAGCCGCTG‑GCAATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTc                                                                                                                                                         <  2:83350/236‑1 (MQ=255)                                                                                  cGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTcc                                                                                                                                                        >  2:610847/1‑239 (MQ=255)                                                                                     cAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCAta                                                                                                                                                     >  2:284449/1‑238 (MQ=255)                                                                                     cAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTACCTGATGCCATGACCTTCCATc                                                                                                                                                     >  2:58135/1‑239 (MQ=255)                                                                                     cAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACAACGGAGAGCAGATGGAGGTAGCTGATGCCAGGACCTTCCATc                                                                                                                                                     >  2:863649/1‑239 (MQ=255)                                                                                             gccaagaagaaccAGCCAGGCGGACTTTTTTGCACGGGTTTTGATTTC‑AATCTG‑CTTTTCCCC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAAGGGGATCGTAAGAGAGTTGGGAGTTG‑AACAGCATAGGGACCAGGA‑AGGACATCACCAAAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAg                                                                                                                                              <  2:426701/227‑1 (MQ=255)                                                                                             cccccggagGCACAGCCGGAGGGACTTTTTGGCCCGGAATGTTTTTTC‑ATTCTG‑CTGTGCCCG‑GCGC‑AGTGAAACCAGGCCT‑TCTGCGAGCGGCTA‑GCGAAAAATGGGTGCGTGAGAGGGTTTGGAGTGG‑AAAAGACGAGGGACAAGGA‑AGGCAATCCCAAGAGCCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAg                                                                                                                                              <  1:147064/232‑1 (MQ=255)                                                                                             ccccAGCAGGAACCGCCAGACGGCCTTTTTTGCACGGGATTTGATTTC‑AATCTG‑CTTTCCCCC‑GCGC‑AGTTCAACCAGGCCC‑CCTTCCAGACGCTT‑GCGATAAACCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAg                                                                                                                                              <  2:572885/238‑1 (MQ=255)                                                                                             ccccAGCAGAAACAGCCAGACGGTCTTTTTTGCCGCGGTTTTTATTTC‑‑ATGTGTCTTTTCACC‑GGGC‑AGTCCAACCAGGCCC‑TCTCCCAGACGCTT‑GCGAAAAATCGGATCGACAGAGAGTTTGGATTTG‑ACCGGAATAGGGACCAGGA‑AGGACATCCCCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAg                                                                                                                                              <  1:694991/238‑1 (MQ=255)                                                                                                   cagGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑CGGACATCCCCCGAACCCCCCGCAGGCGGGACGGGATGACCCCGGCGACCAAATGGAGGTAGCTGATGCCCTGACCTTCCACCAAGGAAGACAttt                                                                                                                                        >  1:7790/1‑237 (MQ=255)                                                                                                     ggAACAGCCAGACGGACGTTTTTGCCCCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACGAGGCCC‑TCTTCCAGCCGCTT‑GCGATAAAGCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTa                                                                                                                                       <  2:643752/238‑1 (MQ=255)                                                                                                      gAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGAAGCTGATGCCATACCCATCAACCACGGAAAAAAAGTaaa                                                                                                                                     >  1:426701/1‑239 (MQ=255)                                                                                                            ccAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACAAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGCAGGTAGCTGATGCAATGACCTTCCATCACGGAAGACCTGTAAACCACt                                                                                                                                >  2:447666/1‑238 (MQ=255)                                                                                                             cAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGCTGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGaaaaaaaatg                                                                                                                               >  1:235554/1‑232 (MQ=255)                                                                                                             cAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGGGACCAGGA‑AGGACATCACCCGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGCAGCTGATGCCATGACCTTCCATAACGGacaacagttaaaacaagg                                                                                                                               >  1:151912/1‑221 (MQ=255)                                                                                                                acggacgTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGAAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAAcaccagcccc                                                                                                                           >  1:178989/1‑230 (MQ=255)                                                                                                                      tttttgGCACGGGTTTTGATGTC‑ATTCTG‑TTTTCCACG‑GGGA‑AGTACAACGAGGCCC‑TTTTCAGGAGGTTT‑GCGATAAATCGGATCGGAAGAGGGTTGGGAGTTG‑AACAGCATGGGGACGAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGAt                                                                                                                     <  2:401055/238‑1 (MQ=255)                                                                                                                               ccGTATTTGATTAC‑AAGTGT‑CTTCTCCCC‑GCGC‑ATTTCAACAGGGCCT‑TCTTCCAAACGCTT‑GCAATAAAGGGGACGGTAAGGGAGTTGGGAGTGG‑AACAGCATAGTGACCAGGA‑AGGAAATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGc                                                                                                            <  2:875428/239‑1 (MQ=255)                                                                                                                                 ggttttgttttA‑AATGTG‑TTTTCCCCC‑GCGC‑ATTCCAACCGGGCCT‑CCTGCCAGCCGTTG‑GCGATAAATCGGACCGTCAGAGAGTGTGGAGTGA‑AAAAGCATGGGGACCAGGA‑GGACCATCCCGAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCgg                                                                                                          <  1:185176/236‑1 (MQ=255)                                                                                                                                  gatttgatttC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGGAACCCAATGCCGCTGAGATTGGCGATGCGGt                                                                                                         >  2:420151/1‑239 (MQ=255)                                                                                                                                  cgttttcctctcatgcTG‑CTTTCCCCC‑GAGA‑ATCTCACCCAGGCCC‑TCTACCAGACCCGT‑GCATTAAATCGGCTCGTAAGAGAGTTTGGATGTA‑AATGGAAGAGTGACAAGGA‑AGGAAGCCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCgg                                                                                                          <  1:408950/224‑1 (MQ=255)                                                                                                                                                  tctgttttGACCGGCGC‑ATTTCAACGAGGCCC‑TTTTCCAACCGCTT‑GCAAAAAACGGGATCGTAAGGGAGTTTGGAGTTG‑AACAGCATAGGGCCAGGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTgc                                                                                            <  2:358014/235‑1 (MQ=255)                                                                                                                                                   ctgttttttccccggtactttCACCCGGGCCC‑GTTCCCAGCCGTTT‑CCGATAAACGGGACCGGCAGAAGGTTGGGAG‑TG‑AAAAGGATGAGGAACAGGAGAGGACTCCACGGGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTgc                                                                                            <  1:850531/220‑1 (MQ=255)                                                                                                                                                              cccgcgcatctaaCCCGGCCCTTCTCCCAGGCCCTT‑GCGTAAAAGCGGATCGGAAAAGAGTTTGGAGTTA‑AACAGCATAGTAACAAGGA‑AGGACATCCCCAGGACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCa                                                                                  <  1:678369/228‑1 (MQ=255)                                                                                                                                                                gtgcttttCAACCAGGCCC‑TGTTCCAGACGTTT‑GCAATAAAGAGGACCGTCAAAGAGTTTGGAGTTG‑AACAGCATAGTGGCCAGGA‑AGGATATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTg                                                                               <  2:656784/233‑1 (MQ=255)                                                                                                                                                                cgc‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGCGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGCTGACCCCGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGACGACATTGAAACCACTGCCGCTGGACTTGACCATGCGGGGATCCCAATCTTCACGGAAACCCCTGc                                                                              >  1:559226/1‑239 (MQ=255)                                                                                                                                                                cgc‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAGGACATGTAAACCACTGCCGCTGAGATGGGCGATTCGTTGATCGCAATCTGCGCCGATACCACTGc                                                                              >  2:745862/1‑239 (MQ=255)                                                                                                                                                                cgc‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGCGATTGGc                                                                                                                 >  2:256255/1‑204 (MQ=255)                                                                                                                                                                cgc‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCAGCTGAGATTGGCGATGCGGTGATCGCAATCTGCCCGGATACCACTGc                                                                              >  2:311423/1‑239 (MQ=255)                                                                                                                                                                        taaCCAGGCCC‑TTTTCCAGGTGCTTGGCGAAAAAGCGGAGCGTCGGGGGGTTGGGGTTTG‑AACAGCATAGGCACAAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGt                                                                           <  1:142937/235‑1 (MQ=255)                                                                                                                                                                         cgccgtggccttttctcgtgctcTT‑GCGGTAAGTCGGATCGTGAGAG‑GGTTGGAGGTGCAACGGAAGTGAGACCAGGA‑AGGCAATACCCGGAGCCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGa                                                                        <  1:410971/216‑1 (MQ=255)                                                                                                                                                                              ggCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCgg                                                                  >  2:674842/1‑238 (MQ=255)                                                                                                                                                                                ccc‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCAcccacgtagcaagcggccg                                                               >  2:633247/1‑221 (MQ=255)                                                                                                                                                                                ccc‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACg                                                               >  2:402023/1‑239 (MQ=255)                                                                                                                                                                                             cGCTT‑GGAATAACTCGTCACGTCAAATGGTGTGGAGTTA‑AACAGAAAATGGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTaa                                                    <  2:274682/238‑1 (MQ=255)                                                                                                                                                                                                tt‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGCTGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACAGCCGCTGAGATTGGCGATGCGGGGAGCGCAACCTGCGCGGATACCAATCCAGGAGACAACGGTAGGCAAGGTTTAAAGcc                                                >  1:774312/1‑239 (MQ=255)                                                                                                                                                                                                                     ggagagTGTGGAGTAG‑AACAGCATGGGGACCGGGA‑AGGACATCACAAGAACCGCCAGCAGGGCGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAg                            <  2:7790/238‑1 (MQ=255)                                                                                                                                                                                                                           tttGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGCCAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATa                       >  2:564087/1‑238 (MQ=255)                                                                                                                                                                                                                           tttGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATATCCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATa                       >  2:380887/1‑238 (MQ=255)                                                                                                                                                                                                                             tGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGaga                       >  2:582887/1‑234 (MQ=255)                                                                                                                                                                                                                             tGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGAACCCACTGCAGACGAAAGCGGGCGGCAAGGTTTACCCCCTTGTTCCTTCGCAACTCAAGCggaaacg                    >  1:276610/1‑233 (MQ=255)                                                                                                                                                                                                                                   tG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAgt              >  2:589055/1‑239 (MQ=255)                                                                                                                                                                                                                                           cATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTGGCCa       >  2:530606/1‑239 (MQ=255)                                                                                                                                                                                                                                           cATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACAACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCCACTTCACCGGTAACTGGCAGTGTCGCCa       >  2:888876/1‑239 (MQ=255)                                                                                                                                                                                                                                                tGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTc                                                                                                                                                         >  1:800423/1‑88 (MQ=255)                                                                                                                                                                                                                                                tGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGGGTCGCCagagag  >  2:37784/1‑239 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                TCAGGATTGCGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACG‑CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTC‑AATCTG‑CTTTTCACC‑GCGC‑AGTTCAACCAGGCCC‑TCTTCCAGACGCTT‑GCGATAAATCGGATCGTCAGAGAGTTTGGAGTTG‑AACAGCATAGTGACCAGGA‑AGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAG  >  NC_000913/4365985‑4366454

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