BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,382,218	C→A	93.0%	noncoding (42/76 nt)	glyV →	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,382,218	0	C	A	93.0%	448.8 / 21.6	178	noncoding (42/76 nt)	glyV	tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base C (6/8); new base A (67/97); total (73/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

AAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTC  >  CP009273/4382073‑4382362

aaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAgg                                                                                                                                                >  2:907686/1‑148 (MQ=255)

aaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAgg                                                                                                                                                <  1:1522805/148‑1 (MQ=255)

 aaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGACTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAggg                                                                                                                                               >  1:3059767/1‑148 (MQ=255)

 aaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                               >  2:537314/1‑148 (MQ=255)

 aaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAggg                                                                                                                                               >  2:1413095/1‑148 (MQ=255)

 aaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAggg                                                                                                                                               >  1:1337835/1‑148 (MQ=255)

  aaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                                                                                                                                              <  1:848409/148‑1 (MQ=255)

  aaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAgggg                                                                                                                                              >  1:1850633/1‑148 (MQ=255)

  aaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAgggg                                                                                                                                              >  1:1852155/1‑148 (MQ=255)

     aCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTcg                                                                                                                                           >  1:3558063/1‑148 (MQ=255)

     aCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTcg                                                                                                                                           >  1:3078990/1‑148 (MQ=255)

      cTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                          >  1:2571828/1‑148 (MQ=255)

       ttttttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTcgcg                                                                                                                                         >  1:2981793/1‑148 (MQ=255)

        tttttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGa                                                                                                                                        <  1:2530049/148‑1 (MQ=255)

          tttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGt                                                                                                                                      <  2:1296958/148‑1 (MQ=255)

          tttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGt                                                                                                                                      <  2:1158399/148‑1 (MQ=255)

          tttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCAGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGt                                                                                                                                      <  2:1296208/148‑1 (MQ=255)

                gggTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGt                                                                                                                                >  2:2114129/1‑148 (MQ=255)

                 ggTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTAAGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGATGGTATAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGtc                                                                                                                               >  2:70865/1‑148 (MQ=255)

                  gTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGtct                                                                                                                              <  1:887006/148‑1 (MQ=255)

                    tGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCg                                                                                                                            <  2:1561886/148‑1 (MQ=255)

                     gCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGt                                                                                                                           <  1:1976010/148‑1 (MQ=255)

                     gCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGt                                                                                                                           <  1:1654977/148‑1 (MQ=255)

                       agagGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGttt                                                                                                                         <  1:2761738/148‑1 (MQ=255)

                        gagGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                                                             >  1:1508794/1‑143 (MQ=255)

                        gagGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                                                        <  2:2339339/148‑1 (MQ=255)

                          gggAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTccc                                                                                                                      >  2:1247144/1‑148 (MQ=255)

                             aaaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGGGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCt                                                                                                                   >  1:1995371/1‑148 (MQ=255)

                                 aTTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaa                                                                                                               >  1:3245372/1‑148 (MQ=255)

                                  tttCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaa                                                                                                              <  1:1178046/148‑1 (MQ=255)

                                  tttCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaa                                                                                                              <  2:1098807/148‑1 (MQ=255)

                                  tttCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaa                                                                                                              <  2:2393361/148‑1 (MQ=255)

                                       cGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTg                                                                                                         >  1:2018724/1‑148 (MQ=255)

                                         tataATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaa                                                                                                       <  2:2983436/148‑1 (MQ=255)

                                          ataATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaa                                                                                                      >  1:670310/1‑148 (MQ=255)

                                            aaTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTTTGAAAAg                                                                                                    >  2:1732399/1‑148 (MQ=255)

                                             aTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaactgt                                                                                                   >  2:2427183/1‑144 (MQ=255)

                                              tGCGCCTCCCGTAACGACGCATAAATGCGAAAATTAAGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtg                                                                                                  <  2:3235652/148‑1 (MQ=255)

                                              tGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtg                                                                                                  >  1:144863/1‑148 (MQ=255)

                                              tGCGCCTCCAGTAACGCCGCAGAATTGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCTCTCCAAAATTTGAAAAGtg                                                                                                  <  2:3819796/148‑1 (MQ=21)

                                                  ccTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtgctgc                                                                                              >  1:3139246/1‑148 (MQ=255)

                                                    tCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAACATTTGAAAAGTGCTGCaa                                                                                            >  2:2549645/1‑148 (MQ=255)

                                                    tCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCaa                                                                                            <  2:3746445/148‑1 (MQ=255)

                                                     cccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCaaa                                                                                           >  2:3407767/1‑148 (MQ=255)

                                                        gTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGca                                                                                        <  1:2114129/148‑1 (MQ=255)

                                                         tAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcac                                                                                       <  1:3069466/148‑1 (MQ=255)

                                                         tAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcac                                                                                       <  1:726033/148‑1 (MQ=255)

                                                          aaCGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcaca                                                                                      <  2:1619018/148‑1 (MQ=255)

                                                            cgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGACTAGCTCAGTTGGTAGATCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGa                                                                                    >  2:2686777/1‑148 (MQ=255)

                                                            cgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGATCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGa                                                                                    >  2:2221970/1‑148 (MQ=255)

                                                            cgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGa                                                                                    <  2:3587380/148‑1 (MQ=255)

                                                               cgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGacca                                                                                 >  1:466193/1‑148 (MQ=255)

                                                                   gAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGAACAGACCACCCa                                                                             >  1:474680/1‑148 (MQ=255)

                                                                          gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGa                                                                      <  1:2250357/148‑1 (MQ=255)

                                                                          gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGa                                                                      <  1:2250806/148‑1 (MQ=255)

                                                                          gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGa                                                                      <  1:1018556/148‑1 (MQ=255)

                                                                            aaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAGGCACAGACCACCCAAGCGGGAAt                                                                    <  2:1786353/148‑1 (MQ=255)

                                                                            aaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAAt                                                                    <  1:984862/148‑1 (MQ=255)

                                                                            aaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAAt                                                                    <  1:3573923/148‑1 (MQ=255)

                                                                               ttACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGc                                                                 <  2:2304020/148‑1 (MQ=255)

                                                                                   gAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAg                                                             <  1:3496453/148‑1 (MQ=255)

                                                                                   gAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAg                                                             <  1:3818217/148‑1 (MQ=255)

                                                                                       gCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCGTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGt                                                        <  2:1782326/148‑1 (MQ=255)

                                                                                        cAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGt                                                        <  1:3226144/148‑1 (MQ=255)

                                                                                          aaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGACTAGCTCAGTTGGTag                                                      >  2:1367389/1‑148 (MQ=255)

                                                                                          aaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGATCAGTTGGTag                                                      >  1:1311814/1‑148 (MQ=255)

                                                                                            aTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTagag                                                    <  1:3529125/148‑1 (MQ=255)

                                                                                             ttAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGtc                                                                                                                               >  1:1716165/1‑72 (MQ=255)

                                                                                             ttAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaa                                                                                                               >  2:2824842/1‑88 (MQ=255)

                                                                                             ttAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGc                                                   <  1:620776/148‑1 (MQ=255)

                                                                                              tAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaattc                                                                                                           >  2:3094682/1‑88 (MQ=255)

                                                                                                agtagtACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTTTGAACTGACGTGAAAAGCACAGACCACCCGAGCGGGAATAGCTCGGTTGGGGTAGCACg                                                >  2:2457856/1‑148 (MQ=14)

                                                                                                  tagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaa                                                                                                               <  2:1194018/83‑1 (MQ=38)

                                                                                                  tagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaa                                                                                                               >  1:1194018/1‑83 (MQ=38)

                                                                                                   agtaCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGAcc                                             <  1:901473/148‑1 (MQ=255)

                                                                                                       cgcgGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                         >  1:3095272/1‑148 (MQ=255)

                                                                                                             aaTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAATTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGt                                   >  2:3422386/1‑148 (MQ=255)

                                                                                                               gaGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCg                                 <  2:2867323/147‑1 (MQ=255)

                                                                                                                aGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                <  2:1402003/148‑1 (MQ=255)

                                                                                                                aGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                <  1:2791957/148‑1 (MQ=255)

                                                                                                                aGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                <  1:1723232/148‑1 (MQ=255)

                                                                                                                  cTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAATTGATGAAGTGCACAGACCACCCAACCGGGAATAGCTCAGTTGGTGGTGCAGGACCTTGCCAAGGtggtgg                              >  2:411135/1‑143 (MQ=11)

                                                                                                                  cTCAGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                              <  1:848249/148‑1 (MQ=255)

                                                                                                                     aGTTGGTAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcg                           <  2:2273160/148‑1 (MQ=255)

                                                                                                                           tAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGtt                     >  2:3179258/1‑148 (MQ=255)

                                                                                                                           tAGAGCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTGAAAAGTGATGAAAAGCACAGCCCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTAGCAAACGACGGGGACGCGAGtg                     >  2:3179196/1‑147 (MQ=11)

                                                                                                                              agCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTTCAAAGCACACACCACCCAAGCGGTAATAGCTCGTTTGGTATGGCACGACCTTCCCAAGGTGGGGGTCGCGGGTTCgt                  >  2:3313402/1‑147 (MQ=14)

                                                                                                                               gCACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                 <  1:2883279/148‑1 (MQ=255)

                                                                                                                                cACGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGCCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                <  2:1425527/148‑1 (MQ=255)

                                                                                                                                 aCGACCTTGCCAAGGTAGGGGTCGCGAGTTCTAGTCTCGTTTCACGCTCCAAAATTTGAAAAGTTCTGCAAAGCACAGACCTCCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc               <  1:1059504/148‑1 (MQ=14)

                                                                                                                                 aCGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc               <  1:3029658/148‑1 (MQ=255)

                                                                                                                                 aCGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc               <  2:3558063/148‑1 (MQ=255)

                                                                                                                                 aCGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc               <  1:1747769/148‑1 (MQ=255)

                                                                                                                                 aCGACCTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc               <  1:1661522/148‑1 (MQ=255)

                                                                                                                                     ccTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGTTTAAAAGCAAAAACCACCCAAGCGGGAATTGCGCCGTTGTTTTTGCACTACCTTCCCAATTTAGTGGTCGCGAGTACCAGACTCGt           >  1:1320048/1‑148 (MQ=2)

                                                                                                                                     ccTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGGTGCAAATCAAAGACCACCCAAGCGGGACTAGCTCAGGTGGGAGTGCACGCCCTGGCCGAGGTCGGGGGCGCGGGTGcgcgtcgcgt           >  2:1977858/1‑148 (MQ=14)

                                                                                                                                     ccTTGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAACGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGt           <  2:1850682/148‑1 (MQ=255)

                                                                                                                                       ttGCCAAGGTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACACACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAACGTCGGGGTCGCGAGTTCGAGTCTCGttt         >  1:2668424/1‑148 (MQ=255)

                                                                                                                                              ggTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc  <  1:1264420/148‑1 (MQ=255)

                                                                                                                                              ggTAGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACAACCCAAGCGGGAATAGCTCAGTGGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc  >  1:2416329/1‑148 (MQ=255)

AAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTC  >  CP009273/4382073‑4382362

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