BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,548,216	A→T	100%	intergenic (‑121/+13)	adhP ← / ← maeA	ethanol‑active dehydrogenase/acetaldehyde‑active reductase/malate dehydrogenase, (decarboxylating, NAD‑requiring) (malic enzyme)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,548,216	0	A	T	98.0%	676.5 / 1.0	210	intergenic (‑121/+13)	adhP/maeA	ethanol‑active dehydrogenase/acetaldehyde‑active reductase/malate dehydrogenase, (decarboxylating, NAD‑requiring) (malic enzyme)
Reads supporting (aligned to +/- strand): ref base A (4/1); new base T (112/93); total (116/94)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 210 total aligned reads displayed.

TCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGC  >  CP009273/1548074‑1548360

tCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCggg                                                                                                                                             <  1:851181/148‑1 (MQ=255)

   tAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGgcgc                                                                                                                                          >  1:310795/1‑148 (MQ=255)

   tAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGAAGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGgcgc                                                                                                                                          >  1:2553791/1‑148 (MQ=255)

    aacaacTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCa                                                                                                                                         >  2:2659957/1‑148 (MQ=255)

        acTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCt                                                                                                                                     >  1:687601/1‑148 (MQ=255)

          tGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTa                                                                                                                                   <  2:2937944/148‑1 (MQ=255)

            cAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTAGCCACTACCCAGTATGGATGAGATCTGAAAATGGGAGAGGGAAATATCCAGGTAGCCGTCACTTACGGGCGCAGCCTTAGa                                                                                                                                 >  2:2076354/1‑148 (MQ=255)

               ccTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATgg                                                                                                                              >  2:2298381/1‑148 (MQ=255)

               ccTTCAAAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATgg                                                                                                                              >  2:1058185/1‑148 (MQ=255)

                 ttCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAg                                                                                                                            <  1:1827363/148‑1 (MQ=255)

                   cATAGGTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGt                                                                                                                          <  2:628908/148‑1 (MQ=255)

                      aGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTAcg                                                                                                                       >  1:2909833/1‑148 (MQ=255)

                         tcctccTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTAcggcg                                                                                                                    >  1:1413690/1‑148 (MQ=255)

                         tcctccTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTAcggcg                                                                                                                    >  1:1499576/1‑148 (MQ=255)

                         tcctccTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTAcggcg                                                                                                                    <  1:3291823/148‑1 (MQ=255)

                          cctcctTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTAcggcgg                                                                                                                   <  1:2139193/148‑1 (MQ=255)

                           ctcctTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGt                                                                                                                  <  1:1928050/148‑1 (MQ=255)

                           ctcctTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGt                                                                                                                  >  1:1237203/1‑148 (MQ=255)

                           ctcctTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGt                                                                                                                  <  2:1621828/148‑1 (MQ=255)

                             cctTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAg                                                                                                                <  1:628443/148‑1 (MQ=255)

                              ctTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGt                                                                                                               >  2:2896195/1‑148 (MQ=255)

                                tttCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTcg                                                                                                             >  1:1228442/1‑148 (MQ=255)

                                   cGGATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCgg                                                                                                          <  1:2192501/148‑1 (MQ=255)

                                    ggATGATGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGt                                                                                                         <  1:1652821/148‑1 (MQ=255)

                                     gatgatGTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTa                                                                                                        <  2:2134682/148‑1 (MQ=255)

                                       tgatgTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGGACGGCGGTAGTCGCGGTAtt                                                                                                      <  2:578742/148‑1 (MQ=255)

                                        gatgTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTc                                                                                                     >  2:3089560/1‑148 (MQ=255)

                                        gatgTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTc                                                                                                     >  1:1231495/1‑148 (MQ=255)

                                         atgTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCg                                                                                                    >  2:1930518/1‑148 (MQ=255)

                                          tgTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCgg                                                                                                   >  2:2170421/1‑148 (MQ=255)

                                          tgTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCgg                                                                                                   <  1:2055855/148‑1 (MQ=255)

                                           gTTCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGc                                                                                                  >  2:3113180/1‑148 (MQ=255)

                                             tCTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCtt                                                                                                >  2:3360201/1‑148 (MQ=255)

                                              cTGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTg                                                                                               >  2:2578657/1‑148 (MQ=255)

                                               tGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGc                                                                                              <  2:863698/148‑1 (MQ=255)

                                               tGCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGc                                                                                              <  2:2861511/148‑1 (MQ=255)

                                                gCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGcc                                                                                             >  1:1396958/1‑148 (MQ=255)

                                                gCATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGcc                                                                                             >  2:3481763/1‑148 (MQ=255)

                                                 cATAGCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCa                                                                                            >  2:1061401/1‑148 (MQ=255)

                                                     gCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGaaa                                                                                        >  1:1007724/1‑148 (MQ=255)

                                                     gCAGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGaaa                                                                                        <  1:1931434/148‑1 (MQ=255)

                                                       aGGTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAAtt                                                                                      <  1:1315253/148‑1 (MQ=255)

                                                        ggTGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTa                                                                                     <  2:566669/148‑1 (MQ=255)

                                                          tGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCTGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATc                                                                                   >  1:2649848/1‑148 (MQ=255)

                                                          tGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATc                                                                                   >  2:2552832/1‑148 (MQ=255)

                                                          tGAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATc                                                                                   <  1:776017/148‑1 (MQ=255)

                                                           gAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCg                                                                                  <  2:2975037/148‑1 (MQ=255)

                                                           gAGGCAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCg                                                                                  <  2:1348358/148‑1 (MQ=255)

                                                               cAAATGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCaa                                                                              <  2:1305193/148‑1 (MQ=255)

                                                                aaaTGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAAt                                                                             >  2:1544629/1‑148 (MQ=255)

                                                                aaaTGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAAt                                                                             <  2:1129011/148‑1 (MQ=255)

                                                                  aTGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATgg                                                                           >  2:2967309/1‑148 (MQ=255)

                                                                   tGATATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGc                                                                          <  1:97604/148‑1 (MQ=255)

                                                                     ataTTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCt                                                                        >  1:268795/3‑148 (MQ=255)

                                                                     ataTTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCt                                                                        >  1:2215476/3‑148 (MQ=255)

                                                                      taTTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTg                                                                       >  2:1811197/2‑148 (MQ=255)

                                                                         ttATTCGCCACTACCCAGTAAGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCCGAACTTATCGTCAATGGCCTGGTg                                                                    >  2:2424035/1‑148 (MQ=255)

                                                                           aTTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCa                                                                  >  1:822961/1‑148 (MQ=255)

                                                                             tCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAgg                                                                >  2:2344933/1‑148 (MQ=255)

                                                                               gCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGc                                                              <  2:146361/148‑1 (MQ=255)

                                                                                ccACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCt                                                             >  1:253801/1‑148 (MQ=255)

                                                                                ccACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCt                                                             <  1:547644/148‑1 (MQ=255)

                                                                                 cACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCtt                                                            >  2:1259901/1‑148 (MQ=255)

                                                                                    tACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCgg                                                         <  1:2531487/148‑1 (MQ=255)

                                                                                       ccAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCag                                                      >  2:27189/1‑148 (MQ=255)

                                                                                        cAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCaga                                                     >  1:395449/1‑148 (MQ=255)

                                                                                          gTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAgg                                                   <  1:2085142/148‑1 (MQ=255)

                                                                                            aTGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGtt                                                 <  2:687601/148‑1 (MQ=255)

                                                                                               gATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTc                                              >  2:1588742/1‑148 (MQ=255)

                                                                                                aagagaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCa                                             <  1:232459/146‑1 (MQ=255)

                                                                                                aTGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCa                                             <  1:2479284/148‑1 (MQ=255)

                                                                                                  gagaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCAcc                                           >  2:1074824/1‑148 (MQ=255)

                                                                                                      tCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCa                                       <  2:1934269/148‑1 (MQ=255)

                                                                                                      acTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCa                                       <  2:1194218/147‑1 (MQ=255)

                                                                                                           aaaaGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCt                                  >  2:465878/1‑148 (MQ=255)

                                                                                                              aGGGAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTtgc                               >  1:146043/1‑148 (MQ=255)

                                                                                                                ggAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTtgctg                             >  2:2029861/1‑148 (MQ=255)

                                                                                                                ggAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTtgctg                             >  1:2942454/1‑148 (MQ=255)

                                                                                                                ggAGAGGGAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTtgctg                             >  1:379308/1‑148 (MQ=255)

                                                                                                                     gggAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTgcg                        >  1:2801513/1‑148 (MQ=255)

                                                                                                                      ggAAATAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTgcgc                       >  1:1829199/1‑148 (MQ=255)

                                                                                                                      ggAAATAGCCCGGTAGCCTTCACGTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTgcgc                       <  2:2400113/148‑1 (MQ=255)

                                                                                                                           tAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCAttt                  <  2:1390457/148‑1 (MQ=255)

                                                                                                                           tAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCAttt                  <  1:1808848/148‑1 (MQ=255)

                                                                                                                           tAGCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCAttt                  >  1:843142/1‑148 (MQ=255)

                                                                                                                             gCCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTg                <  2:915168/148‑1 (MQ=255)

                                                                                                                              cccGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGc               >  1:621188/1‑148 (MQ=255)

                                                                                                                              cccGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGc               >  1:595894/1‑148 (MQ=255)

                                                                                                                              cccGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGc               >  1:2360726/1‑148 (MQ=255)

                                                                                                                              cccGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGc               <  2:811551/148‑1 (MQ=255)

                                                                                                                              cccGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGc               >  1:1187156/1‑148 (MQ=255)

                                                                                                                                cGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCa             >  2:3452465/1‑148 (MQ=255)

                                                                                                                                cGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCa             <  1:3177148/148‑1 (MQ=255)

                                                                                                                                 ggTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCaa            >  2:2777270/1‑148 (MQ=255)

                                                                                                                                 ggTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCaa            <  1:713321/148‑1 (MQ=255)

                                                                                                                                    aGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCg         >  1:2036571/1‑148 (MQ=255)

                                                                                                                                     gCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGc        >  1:3314216/1‑148 (MQ=255)

                                                                                                                                      ccTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCa       >  1:384731/1‑148 (MQ=255)

                                                                                                                                       cTTCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCaa      <  2:1467262/148‑1 (MQ=255)

                                                                                                                                        ttCACTTCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCTGTATTCGGCTTTCCAGAAATTATCGTCAATGGCCTGTTGCAGGGATTCGGCAGAGGTTTTCACCGCCACGCGTTGCTGCTGCGCCATTTTGCCAACCGCtaa     <  2:1228999/148‑1 (MQ=255)

TCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGC  >  CP009273/1548074‑1548360

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