BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	372,581	G→T	85.3%	A118A (GCG→GCT)	yaiL →	DUF2058 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	372,581	0	G	T	85.3%	410.9 / 76.3	184	A118A (GCG→GCT)	yaiL	DUF2058 family protein
Reads supporting (aligned to +/- strand): ref base G (16/11); new base T (77/80); total (93/91)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 184 total aligned reads displayed.

AAAGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACG  >  CP009273/372435‑372723

 acgcccAGGTGAAGCAGCTTATTGACCTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTAt                                                                                                                                              <  2:3300244/143‑1 (MQ=255)

 aaGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTAt                                                                                                                                              >  2:2710552/1‑148 (MQ=255)

 aaGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTAt                                                                                                                                              >  2:1195000/1‑148 (MQ=255)

 aaGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGAt                                                                                                                                              <  2:3236571/148‑1 (MQ=255)

   gCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTg                                                                                                                                            <  2:1891499/148‑1 (MQ=255)

   gCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTg                                                                                                                                            <  2:221977/148‑1 (MQ=255)

       aGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCg                                                                                                                                        <  1:2284075/148‑1 (MQ=255)

        ggTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGc                                                                                                                                       <  1:2365465/148‑1 (MQ=255)

            aaGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCttgt                                                                                                                                   >  1:2421526/1‑148 (MQ=255)

            aaGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCttgt                                                                                                                                   >  2:3114209/1‑148 (MQ=255)

               cagcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTgg                                                                                                                                >  1:3063585/1‑148 (MQ=255)

                agcTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGt                                                                                                                               <  1:460158/148‑1 (MQ=255)

                  cTTATTGAAATGAACCGAATCACCCTTGCCAATGGCGATATTGTTTTTAACTTCACTGACGGGAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCg                                                                                                                             <  2:2559975/148‑1 (MQ=255)

                  cTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCg                                                                                                                             <  2:1515164/148‑1 (MQ=255)

                  cTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCg                                                                                                                             <  2:2287658/148‑1 (MQ=255)

                     aTTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATa                                                                                                                          <  1:1258762/148‑1 (MQ=255)

                       tGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAAc                                                                                                                        >  1:136666/1‑148 (MQ=255)

                       tGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAAc                                                                                                                        >  1:136522/1‑148 (MQ=255)

                        gAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAGCa                                                                                                                       <  2:1178197/148‑1 (MQ=255)

                        gAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACa                                                                                                                       <  2:109603/148‑1 (MQ=255)

                         aaaTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACaa                                                                                                                      >  1:632096/1‑148 (MQ=255)

                           aTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATa                                                                                                                    >  2:382511/1‑148 (MQ=255)

                           aTGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATa                                                                                                                    >  2:3149581/1‑148 (MQ=255)

                              aaCCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCg                                                                                                                 >  1:927182/1‑148 (MQ=255)

                               aCCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGa                                                                                                                >  1:1380516/1‑148 (MQ=255)

                               aCCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGa                                                                                                                <  2:1736678/148‑1 (MQ=255)

                               aCCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGa                                                                                                                <  2:1736272/148‑1 (MQ=255)

                               aCCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGa                                                                                                                <  2:1465409/148‑1 (MQ=255)

                                 cGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAg                                                                                                              <  2:205406/148‑1 (MQ=255)

                                   aaTCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGt                                                                                                            <  2:3126183/148‑1 (MQ=255)

                                   aaTCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGt                                                                                                            >  1:1039615/1‑148 (MQ=255)

                                     tCACCATTGCCAATGGCGATATTGTTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTTATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGa                                                                                                          <  1:1767177/148‑1 (MQ=255)

                                     tCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGa                                                                                                          >  1:2892071/1‑148 (MQ=255)

                                     tCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGa                                                                                                          >  2:2413685/1‑148 (MQ=255)

                                     tCACCATTGCCAAGGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGa                                                                                                          <  1:1767226/148‑1 (MQ=255)

                                      cACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGaa                                                                                                         >  2:1687603/1‑148 (MQ=255)

                                      cACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGaa                                                                                                         >  1:2084608/1‑148 (MQ=255)

                                      cACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGaa                                                                                                         >  2:66057/1‑148 (MQ=255)

                                       aCCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAAt                                                                                                        <  2:3141931/148‑1 (MQ=255)

                                        ccATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATa                                                                                                       >  1:663709/1‑148 (MQ=255)

                                        ccATTGCCAATGGCGATATTGGTTTTAACTTCACGGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATa                                                                                                       <  2:98877/148‑1 (MQ=255)

                                         cATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATAc                                                                                                      >  2:367360/1‑148 (MQ=255)

                                           ttGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGc                                                                                                    >  2:777514/1‑148 (MQ=255)

                                           ttGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGc                                                                                                    >  2:2802565/1‑148 (MQ=255)

                                           ttGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGc                                                                                                    >  1:1890845/1‑148 (MQ=255)

                                            tGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGcc                                                                                                   >  1:2450472/1‑148 (MQ=255)

                                             gCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCa                                                                                                  <  2:2385650/148‑1 (MQ=255)

                                               cAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCatt                                                                                                >  2:2991128/1‑148 (MQ=255)

                                               cAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCatt                                                                                                >  2:2170306/1‑148 (MQ=255)

                                                aaTGGCGATATTGTTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCatta                                                                                               <  1:698465/148‑1 (MQ=255)

                                                aaTGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCatta                                                                                               <  2:3097695/148‑1 (MQ=255)

                                                aaTGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCatta                                                                                               >  2:1290181/1‑148 (MQ=255)

                                                  tGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCattatt                                                                                             <  1:1766810/148‑1 (MQ=255)

                                                   ggCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTc                                                                                            <  1:1267055/148‑1 (MQ=255)

                                                    gCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTcc                                                                                           <  1:258776/148‑1 (MQ=255)

                                                        tatTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGcc                                                                                       >  1:1042520/1‑148 (MQ=255)

                                                          ttGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAg                                                                                     >  1:2844259/1‑148 (MQ=255)

                                                              ttttAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGtt                                                                                 >  1:2703343/1‑148 (MQ=255)

                                                               tttAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTg                                                                                <  2:148150/148‑1 (MQ=255)

                                                               tttAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTg                                                                                >  1:2595616/1‑148 (MQ=255)

                                                                   aCTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGa                                                                            >  2:1930921/1‑148 (MQ=255)

                                                                    cTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGAt                                                                           >  1:848361/1‑148 (MQ=255)

                                                                           gACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTg                                                                    >  1:943797/1‑148 (MQ=255)

                                                                           gACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTg                                                                    <  2:3071450/148‑1 (MQ=255)

                                                                           gACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTg                                                                    <  2:2421526/148‑1 (MQ=255)

                                                                            aCGGCAATCTTATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGATATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGc                                                                   <  1:646398/148‑1 (MQ=255)

                                                                                cAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATATAATTGCTCAg                                                               >  1:2140918/1‑148 (MQ=255)

                                                                                cAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAg                                                               >  2:2740789/1‑148 (MQ=255)

                                                                                    cTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGcgcg                                                           <  2:1712118/148‑1 (MQ=255)

                                                                                    cTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGcgcg                                                           <  1:2468586/148‑1 (MQ=255)

                                                                                       aTTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATg                                                        >  2:583872/1‑148 (MQ=255)

                                                                                       aTTAAGAAGATTTTTGTAGATAAGCTCACACAGTCGCAGTTGAATAATGTTCGTCTTGCTATTTCCCGCTTGTTGGTCGATAAAAATATCGAAGGTGAAAACTCCATTATTCCCTCAAGCGATGACGATAAAAATGCTCAGATCGATg                                                        >  2:2508803/1‑148 (MQ=255)

                                                                                           agaagaTTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCcag                                                    >  2:366277/1‑148 (MQ=255)

                                                                                               gaTTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGt                                                >  2:1314437/1‑148 (MQ=255)

                                                                                               gaTTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGt                                                >  1:1399087/1‑148 (MQ=255)

                                                                                                  ttttGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGtata                                             >  1:728294/1‑148 (MQ=255)

                                                                                                  ttttGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGtata                                             >  1:728291/1‑148 (MQ=255)

                                                                                                  ttttGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGtata                                             <  2:259094/148‑1 (MQ=255)

                                                                                                     tGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAgtg                                          <  2:3042146/148‑1 (MQ=255)

                                                                                                         gATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGc                                      <  2:420025/148‑1 (MQ=255)

                                                                                                           tAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGcac                                    <  2:2509684/148‑1 (MQ=255)

                                                                                                             aGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAg                                  <  2:2226204/148‑1 (MQ=255)

                                                                                                               cTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAgcg                                <  1:233909/148‑1 (MQ=255)

                                                                                                                 cacacaGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAgcgcg                              >  1:863117/1‑148 (MQ=255)

                                                                                                                 cacacaGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAgcgcg                              <  1:1587900/148‑1 (MQ=255)

                                                                                                                  acacaGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAgcgcgc                             <  2:3272220/148‑1 (MQ=255)

                                                                                                                   cacaGGCGCAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCt                            >  2:167732/1‑148 (MQ=255)

                                                                                                                        gcgcAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGAGTTGCACAGCGCGCTCAgcg                       >  1:2134734/1‑148 (MQ=255)

                                                                                                                           cAGTTGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCg                    >  1:1113046/1‑148 (MQ=255)

                                                                                                                              ttGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAag                 >  1:1292059/1‑148 (MQ=255)

                                                                                                                              ttGATTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAag                 >  1:2183160/1‑148 (MQ=255)

                                                                                                                               tGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAaga                >  1:2649145/1‑148 (MQ=255)

                                                                                                                                 aTTAATGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGc              <  1:3281179/148‑1 (MQ=255)

                                                                                                                                    aaTGGTCTTGTGGCTATTGCCCGCTTGTTGGTCTATAAAATTAGCGAAGGTGAATACTCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCATCTGTATAGTGTTGCACAGCGAGCTCAGCGCCGAAGAGCAgg           <  1:2508803/148‑1 (MQ=255)

                                                                                                                                    aaTGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAgg           <  1:583872/148‑1 (MQ=255)

                                                                                                                                    aaTGGTCGTCTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAgg           <  1:2246676/148‑1 (MQ=255)

                                                                                                                                        gtcgtcTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGa       <  2:1441699/148‑1 (MQ=255)

                                                                                                                                        gtcgtcTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGa       <  2:1428771/148‑1 (MQ=255)

                                                                                                                                        gtcgtcTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGa       <  1:3149581/148‑1 (MQ=255)

                                                                                                                                           gtcTGGCTATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAga    <  1:66057/148‑1 (MQ=255)

AAAGCTCAGGTGAAGCAGCTTATTGAAATGAACCGAATCACCATTGCCAATGGCGATATTGGTTTTAACTTCACTGACGGCAATCTGATTAAGAAGATTTTTGTCGATAAGCTCACACAGGCGCAGTTGATTAATGGTCGTCTGGCGATTGCCCGCTTGTTGGTCGATAACAATAGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACG  >  CP009273/372435‑372723

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: `‑`