BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,287,829	C→A	94.7%	S185S (TCG→TCT)	ccmF ←	heme lyase, CcmF subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,287,829	0	C	A	94.7%	520.0 / 18.8	179	S185S (TCG→TCT)	ccmF	heme lyase, CcmF subunit
Reads supporting (aligned to +/- strand): ref base C (5/5); new base A (83/86); total (88/91)
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 180 total aligned reads displayed.

AGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAG  >  CP009273/2287683‑2287972

aGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAga                                                                                                                                               <  2:914272/149‑1 (MQ=255)

 gTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGATTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACagag                                                                                                                                              <  2:1159250/149‑1 (MQ=255)

   tCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGaa                                                                                                                                            >  1:800973/1‑149 (MQ=255)

    cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGaaa                                                                                                                                           <  2:3554631/149‑1 (MQ=255)

      tAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAAcc                                                                                                                                         >  1:1143796/1‑149 (MQ=255)

      tAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAAcc                                                                                                                                         >  1:1553214/1‑149 (MQ=255)

       aaTAGGCCCATGCGGAACCGAGCACGATGCCGAGCTTCAGTAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAACCGCGAATAAGTGGTGTCCAGACGCCCGCTCAGAAAAGAAGAAATGGCAAAAGCAAACGCCACAGAGAAAccc                                                                                                                                        <  1:835579/149‑1 (MQ=255)

        aTAGGCCCATGCGGAACCGAGCACGTTGCCGAGCGTCAGTAAGATCCACGCTGCCTGCGTCCACGGACGAGTAAAACGCGAATAAGTGCTGTCCAGACGCCCGATCAGCAAAGAAGCAATGGAAAAAGCAAACGCCACAGAGAAACCCa                                                                                                                                       <  1:3171106/149‑1 (MQ=255)

        aTAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCa                                                                                                                                       >  2:248989/1‑149 (MQ=255)

         tAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCTTCCACTGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGAAAACGCCACAGAGAAACCCAc                                                                                                                                      <  2:3483747/149‑1 (MQ=255)

           ggCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGt                                                                                                                                    >  1:1033323/1‑149 (MQ=255)

           ggCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGt                                                                                                                                    >  1:1033318/1‑149 (MQ=255)

             cccATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATATGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTAc                                                                                                                                  <  1:1516119/149‑1 (MQ=255)

               cATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGAAAATAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTAccc                                                                                                                                <  1:526963/149‑1 (MQ=255)

                 tGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCa                                                                                                                              >  2:773320/1‑149 (MQ=255)

                     gAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCAtata                                                                                                                          <  2:791321/149‑1 (MQ=255)

                     gAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCAtata                                                                                                                          >  1:2706033/1‑149 (MQ=255)

                     gAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCAtata                                                                                                                          >  1:1915616/1‑149 (MQ=255)

                      aaCCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATaa                                                                                                                         >  1:2483100/1‑149 (MQ=255)

                       aCCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATaaa                                                                                                                        >  2:3418299/1‑149 (MQ=255)

                        ccGAGCACGATGCCGAGCGTCAGGTAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAg                                                                                                                       <  2:290786/149‑1 (MQ=255)

                          gAGCACTATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGATTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGAAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCa                                                                                                                     <  1:2616578/149‑1 (MQ=255)

                          gAGCACGATGCCGCGAGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCAAGACGCCCGCTCAGCAAAGAAGAAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCa                                                                                                                     <  1:2616576/149‑1 (MQ=255)

                            ncACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCaga                                                                                                                   >  1:3239732/2‑149 (MQ=255)

                            gCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCTGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGGTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCaga                                                                                                                   <  2:700208/149‑1 (MQ=255)

                            gCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCaga                                                                                                                   >  1:1228155/1‑149 (MQ=255)

                            gCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCaga                                                                                                                   >  1:2654936/1‑149 (MQ=255)

                            gCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCaga                                                                                                                   >  1:3239544/1‑149 (MQ=255)

                            gCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCaga                                                                                                                   <  2:2411729/149‑1 (MQ=255)

                             cACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCagag                                                                                                                  <  2:40200/149‑1 (MQ=255)

                             cACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCagag                                                                                                                  <  2:2118977/149‑1 (MQ=255)

                              aCGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAgg                                                                                                                 >  1:3178267/1‑149 (MQ=255)

                               cGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAGGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGc                                                                                                                <  2:1360239/149‑1 (MQ=255)

                                 aTGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCATCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCgg                                                                                                              >  1:1708662/1‑149 (MQ=255)

                                 aTGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCgg                                                                                                              >  2:1320359/1‑149 (MQ=255)

                                       aGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGaa                                                                                                        >  2:876666/1‑149 (MQ=255)

                                        gCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAg                                                                                                       >  1:1190127/1‑149 (MQ=255)

                                           tCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATc                                                                                                    >  1:1300759/1‑149 (MQ=255)

                                           tCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATc                                                                                                    >  1:3351321/1‑149 (MQ=255)

                                           tCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATc                                                                                                    <  1:2805235/149‑1 (MQ=255)

                                               gAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCAGATAAAGCAGAGGCGGAGGGAAGATCAGcc                                                                                                >  1:3523148/1‑149 (MQ=255)

                                                aaGATCCTCGCTGACAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGccc                                                                                               <  2:3981426/149‑1 (MQ=255)

                                                aaGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGccc                                                                                               <  1:3754086/149‑1 (MQ=255)

                                                aaGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGccc                                                                                               <  1:1709764/149‑1 (MQ=255)

                                                 aGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGcccc                                                                                              >  1:3349566/1‑149 (MQ=255)

                                                 aGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGcccc                                                                                              <  2:3324631/149‑1 (MQ=255)

                                                       aCGCTGCCAGCGTCCACGGACGAGTACAACGCGCATACGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATcc                                                                                        <  2:160854/149‑1 (MQ=255)

                                                       aCGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATcc                                                                                        <  2:2236672/149‑1 (MQ=255)

                                                       aCGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATcc                                                                                        >  1:1062502/1‑149 (MQ=255)

                                                       aCGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATcc                                                                                        <  1:1245532/149‑1 (MQ=255)

                                                         gCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTg                                                                                      >  1:827328/1‑149 (MQ=255)

                                                                gcttCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGc                                                                               <  2:1263691/146‑1 (MQ=255)

                                                                 cGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCg                                                                              <  1:1347811/149‑1 (MQ=255)

                                                                    ccACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGt                                                                           >  2:2850989/1‑149 (MQ=255)

                                                                    ccACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGt                                                                           >  1:45311/1‑149 (MQ=255)

                                                                      aCGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTa                                                                         >  1:142353/1‑149 (MQ=255)

                                                                      aCGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTa                                                                         >  2:3818951/1‑149 (MQ=255)

                                                                        ggACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAg                                                                       <  1:3078538/149‑1 (MQ=255)

                                                                         gACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGa                                                                      >  1:3860693/1‑149 (MQ=255)

                                                                          aCGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGCAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGAt                                                                     <  2:331441/149‑1 (MQ=255)

                                                                          aCGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGAt                                                                     >  1:963485/1‑149 (MQ=255)

                                                                          aCGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGAt                                                                     <  2:2082699/149‑1 (MQ=255)

                                                                          aCGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGAt                                                                     >  1:3502944/1‑149 (MQ=255)

                                                                              gTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATcgcg                                                                 >  2:542350/1‑149 (MQ=255)

                                                                              gTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATcgcg                                                                 >  2:540912/1‑149 (MQ=255)

                                                                               tAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGa                                                                >  1:3819462/1‑149 (MQ=255)

                                                                                aaaaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGACGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAc                                                               >  2:2008420/1‑149 (MQ=255)

                                                                                 aaaCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAcc                                                              <  1:1997561/149‑1 (MQ=255)

                                                                                    cgcgCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTc                                                           <  1:3560282/149‑1 (MQ=255)

                                                                                          tAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCgg                                                     <  2:142353/149‑1 (MQ=255)

                                                                                                cTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGtt                                               <  1:1970321/149‑1 (MQ=255)

                                                                                                   tCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCgg                                            <  2:3356802/149‑1 (MQ=255)

                                                                                                          gCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTg                                     <  2:2706033/149‑1 (MQ=255)

                                                                                                            ccGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCg                                   >  2:3063512/1‑149 (MQ=255)

                                                                                                              gCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCACCGTGCgag                                 >  2:3146481/1‑149 (MQ=255)

                                                                                                              gCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgag                                 <  2:1268757/149‑1 (MQ=255)

                                                                                                                tCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgagag                               >  1:631968/1‑149 (MQ=255)

                                                                                                                 cAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgagaga                              <  2:446893/149‑1 (MQ=255)

                                                                                                                 cAGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgagaga                              >  2:3970775/1‑149 (MQ=255)

                                                                                                                  aGCAAAGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGaa                             <  1:3175037/149‑1 (MQ=255)

                                                                                                                   gCAAAGAAGCAATGGCAAACGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAAc                            <  1:1935270/149‑1 (MQ=255)

                                                                                                                   gCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAAc                            >  1:1414728/1‑149 (MQ=255)

                                                                                                                   gCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAAc                            >  1:1415145/1‑149 (MQ=255)

                                                                                                                   gCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAAc                            >  1:1414559/1‑149 (MQ=255)

                                                                                                                     aaaGAAGCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACgg                          <  2:1632530/149‑1 (MQ=255)

                                                                                                                       agaagCAATGGCAAAAGCAAACGCCACAGAGAAACCCCCGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGt                        <  2:1224412/149‑1 (MQ=255)

                                                                                                                       agaagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGt                        <  2:3830434/149‑1 (MQ=255)

                                                                                                                       agaagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGt                        <  2:3937940/149‑1 (MQ=255)

                                                                                                                         aagCAATGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTa                      <  1:1607257/149‑1 (MQ=255)

                                                                                                                               tGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTa                <  2:1792548/149‑1 (MQ=255)

                                                                                                                               tGGCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTa                <  1:470076/149‑1 (MQ=255)

                                                                                                                                 gCAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTaaa              >  1:3451795/1‑149 (MQ=255)

                                                                                                                                  cAAAAGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAag             >  2:2886194/1‑149 (MQ=255)

                                                                                                                                      aGCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGaa         <  2:3239544/149‑1 (MQ=255)

                                                                                                                                       gCAAACGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAgagact        >  2:3787088/1‑147 (MQ=255)

                                                                                                                                          aaCGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGaa     >  2:519758/1‑149 (MQ=255)

                                                                                                                                          aaCGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGaa     >  2:1650012/1‑149 (MQ=255)

                                                                                                                                            cGCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACa   <  2:817028/149‑1 (MQ=255)

                                                                                                                                             gCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCTGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAAGGGGTTAGAGGTAAAGAGAATGAACAg  <  1:2248149/149‑1 (MQ=255)

                                                                                                                                             gCCACAGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAg  <  1:1599199/149‑1 (MQ=255)

AGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAG  >  CP009273/2287683‑2287972

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