BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,561,616	C→T	6.2%	Q17* (CAA→TAA)	yjiS →	DUF1127 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,561,616	0	C	T	6.2%	1052.2 / 35.8	291	Q17* (CAA→TAA)	yjiS	DUF1127 family protein
Reads supporting (aligned to +/- strand): ref base C (152/121); new base T (9/9); total (161/130)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 292 total aligned reads displayed.

AACCTGTCTGTATGTTAAATAAAAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGATTTATATATAATAAAACG  >  CP009273/4561471‑4561763

aaCCTGTCTGTATGTTAAATAAAAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAAc                                                                                                                                                  <  1:3909723/149‑1 (MQ=255)

aaCCTGTCTGTATGTTAAATAAAAGTAATATTGGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAAc                                                                                                                                                  <  1:1865931/149‑1 (MQ=255)

 aCCTGTCTGTATGTTAAATAAAAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAAct                                                                                                                                                 <  2:1991262/149‑1 (MQ=255)

  ccTGTCTGTATGTTAAATAAAAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAActc                                                                                                                                                >  1:3488549/1‑149 (MQ=255)

  ccTGTCTGTATGTTAAATAAAAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAActc                                                                                                                                                >  1:4648891/1‑149 (MQ=255)

      tctgtATGTTAAATAAAAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGc                                                                                                                                            >  1:1055625/1‑149 (MQ=255)

      tctgtATGTTAAATAAAAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGc                                                                                                                                            >  1:3254471/1‑149 (MQ=255)

       ctgtATTTTAAATAAAATTAATATTTGAATCTTTATTGTTAACTGCCAGAAGGTGGATGATGAAGCATCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCa                                                                                                                                           <  1:1634986/149‑1 (MQ=255)

         gtATGTTAAATAAAAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAgg                                                                                                                                         >  1:1534106/1‑149 (MQ=255)

         gtATGTTAAATAAAAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAgg                                                                                                                                         >  2:3390187/1‑149 (MQ=255)

               taaataaaAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGa                                                                                                                                   >  1:3980426/1‑149 (MQ=255)

               taaataaaAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGa                                                                                                                                   >  2:238535/1‑149 (MQ=255)

                 aataaaAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAgg                                                                                                                                 >  2:3673159/1‑149 (MQ=255)

                 aataaaAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAgg                                                                                                                                 >  1:305770/1‑149 (MQ=255)

                  ataaaAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGc                                                                                                                                >  1:192290/1‑149 (MQ=255)

                   taaaAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGTAACTCTGGCAGGCGGTGAGGCg                                                                                                                               >  2:3035434/1‑149 (MQ=255)

                   taaaAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGTAACTCTGGCAGGCGGTGAGGCg                                                                                                                               >  2:644320/1‑149 (MQ=255)

                   taaaAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCg                                                                                                                               >  2:2134116/1‑149 (MQ=255)

                     aaaGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGtt                                                                                                                             <  2:3680701/149‑1 (MQ=255)

                     aaaGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGtt                                                                                                                             >  1:2359254/1‑149 (MQ=255)

                     aaaCTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGtt                                                                                                                             >  2:2100265/1‑149 (MQ=255)

                       aGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTTGATGATGAAGCCTCTGAAAGATGAGGATGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCTCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGAGTTgg                                                                                                                           >  2:1808550/1‑149 (MQ=255)

                         tAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTggcg                                                                                                                         >  1:1113883/1‑149 (MQ=255)

                          aaTATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTggcgg                                                                                                                        >  1:1825312/1‑149 (MQ=255)

                           atatTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTggcggc                                                                                                                       >  1:573256/1‑149 (MQ=255)

                             atTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTggcggcgg                                                                                                                     >  2:1932403/1‑149 (MQ=255)

                                  aaTCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAAt                                                                                                                >  2:1347538/1‑149 (MQ=255)

                                   aTCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATg                                                                                                               >  1:223172/1‑149 (MQ=255)

                                   aTCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATg                                                                                                               >  1:4283537/1‑149 (MQ=255)

                                    tCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGc                                                                                                              >  2:575435/1‑149 (MQ=255)

                                      tGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAg                                                                                                            >  1:2951583/1‑149 (MQ=255)

                                         aTTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGAcc                                                                                                         >  1:83278/1‑149 (MQ=255)

                                         aTTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGAcc                                                                                                         <  2:4403963/149‑1 (MQ=255)

                                          ttGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGAccc                                                                                                        >  2:2708765/1‑149 (MQ=255)

                                          ttGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGAccc                                                                                                        >  2:4353429/1‑149 (MQ=255)

                                          ttGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGAccc                                                                                                        >  2:4353402/1‑149 (MQ=255)

                                           tgtgAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCg                                                                                                       >  1:2298740/1‑149 (MQ=255)

                                                  tGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCgtg                                                                                                <  1:2527597/149‑1 (MQ=255)

                                                     cAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTa                                                                                             >  2:3117094/1‑149 (MQ=255)

                                                      aGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGTAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTAc                                                                                            <  1:2772728/149‑1 (MQ=255)

                                                       gAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTAca                                                                                           <  1:261250/149‑1 (MQ=255)

                                                            tGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGa                                                                                      <  2:1899809/149‑1 (MQ=255)

                                                            tGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGa                                                                                      <  1:660107/149‑1 (MQ=255)

                                                            tGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGa                                                                                      <  1:2868559/149‑1 (MQ=255)

                                                             ggATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGAt                                                                                     <  2:667657/149‑1 (MQ=255)

                                                              tatgatgAAGACTATTAAAGATGAGGAGTTAGAAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATg                                                                                    <  1:4753924/148‑1 (MQ=255)

                                                              gatgatGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATg                                                                                    <  1:3518955/149‑1 (MQ=255)

                                                               atgatgAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGa                                                                                   <  1:1556806/149‑1 (MQ=255)

                                                                 gatgaAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGt                                                                                 >  1:3651575/1‑149 (MQ=255)

                                                                  atgaAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGtg                                                                                <  2:1108590/149‑1 (MQ=255)

                                                                    gaAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGtgat                                                                              >  1:2285181/1‑149 (MQ=255)

                                                                     aaGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGtgatg                                                                             >  1:1813541/1‑149 (MQ=255)

                                                                      aGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGtgatga                                                                            <  1:4195126/149‑1 (MQ=255)

                                                                       gCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAg                                                                           >  2:2806352/1‑149 (MQ=255)

                                                                        ccTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGc                                                                          >  1:1278284/1‑149 (MQ=255)

                                                                        ccTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGc                                                                          <  2:2435127/149‑1 (MQ=255)

                                                                             gAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGTAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTg                                                                     <  2:3095532/149‑1 (MQ=255)

                                                                               aaGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGaa                                                                   >  2:3554224/1‑149 (MQ=255)

                                                                                 gATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAgg                                                                 <  1:871867/149‑1 (MQ=255)

                                                                                   tGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGat                                                               <  1:3863385/149‑1 (MQ=255)

                                                                                    gaggagGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGata                                                              >  2:1184000/1‑149 (MQ=255)

                                                                                    gaggagGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGata                                                              >  1:1596209/1‑149 (MQ=255)

                                                                                     aggaggTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGatat                                                             >  2:3295404/1‑149 (MQ=255)

                                                                                      ggaggTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATc                                                            >  2:3353511/1‑149 (MQ=255)

                                                                                      ggaggTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATc                                                            <  2:4507682/149‑1 (MQ=255)

                                                                                       gaggTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCg                                                           <  1:4035723/149‑1 (MQ=255)

                                                                                        aggTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCgg                                                          >  1:3110204/1‑149 (MQ=255)

                                                                                         gggaGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCggg                                                         <  2:2616789/146‑1 (MQ=255)

                                                                                           tAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGtt                                                       <  1:1487629/149‑1 (MQ=255)

                                                                                            aGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTa                                                      >  2:2312593/1‑149 (MQ=255)

                                                                                             gCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTAc                                                     >  2:761115/1‑149 (MQ=255)

                                                                                                 aTGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAg                                                 >  1:3051268/1‑149 (MQ=255)

                                                                                                  tGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAgg                                                >  1:1720511/1‑149 (MQ=255)

                                                                                                             gAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGa                                     <  1:679649/149‑1 (MQ=255)

                                                                                                              aacaCAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAg                                    <  2:3937843/146‑1 (MQ=255)

                                                                                                              aaaaCAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCATCTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAg                                    <  1:2067459/149‑1 (MQ=255)

                                                                                                               aaaCAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGt                                   <  2:1655845/149‑1 (MQ=255)

                                                                                                               aaaCAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGt                                   >  1:4349650/1‑149 (MQ=255)

                                                                                                                 aCAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGa                                 >  2:4229636/1‑149 (MQ=255)

                                                                                                                    gagCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAggg                              <  1:2683828/149‑1 (MQ=255)

                                                                                                                     agCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAgggg                             <  2:4584044/149‑1 (MQ=255)

                                                                                                                      gCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAggggg                            >  2:511250/1‑149 (MQ=255)

                                                                                                                        tAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGa                          <  2:4699722/149‑1 (MQ=255)

                                                                                                                         aaaGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGat                         <  2:1047898/149‑1 (MQ=255)

                                                                                                                         aaaGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGat                         <  2:3771516/149‑1 (MQ=255)

                                                                                                                           aGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGatat                       >  2:2224278/1‑149 (MQ=255)

                                                                                                                            gcgcCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGAtata                      <  1:4665179/149‑1 (MQ=255)

                                                                                                                            gcgcCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGAtata                      <  1:1932403/149‑1 (MQ=255)

                                                                                                                            gcgcCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGAtata                      <  2:804357/149‑1 (MQ=255)

                                                                                                                             cgcCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAg                     <  2:4617625/149‑1 (MQ=255)

                                                                                                                               ccGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGAt                   <  2:2926957/149‑1 (MQ=255)

                                                                                                                                cGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGAtt                  >  2:1599925/1‑149 (MQ=255)

                                                                                                                                cGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGAtt                  >  1:4057885/1‑149 (MQ=255)

                                                                                                                                 gTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGAttt                 <  1:2121381/149‑1 (MQ=255)

                                                                                                                                  tttATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGATTta                >  2:3877224/1‑149 (MQ=255)

                                                                                                                                    tATCGGCCTGGTGTAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGATTtata              >  1:1217203/1‑149 (MQ=255)

                                                                                                                                     aTCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGATTtatat             >  2:2852972/1‑149 (MQ=255)

                                                                                                                                              gTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGATTTATATATAATaaaa    >  1:4092384/1‑149 (MQ=255)

                                                                                                                                              gTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGATTTATATATAATaaaa    >  1:4043347/1‑149 (MQ=255)

                                                                                                                                               tGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGATTTATATATAATAAAAc   <  2:1261375/149‑1 (MQ=255)

AACCTGTCTGTATGTTAAATAAAAGTAATATTTGAATCTGTATTGTGAACTGCCAGAAGGTGGATGATGAAGCCTCTGAAAGATGAGGAGGTAGCAGATGGAATTTCACGAAAACAGAGCTAAAGCGCCGTTTATCGGCCTGGTGCAACTCTGGCAGGCGGTGAGGCGTTGGCGGCGGCAAATGCAGACCCGACGCGTGTTACAGCAGATGAGTGATGAGCGGTTGAAGGATATCGGGTTACGCAGGGAGGATGTGGAGTGAGGGGGGATATAGATTTATATATAATAAAACG  >  CP009273/4561471‑4561763

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