BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,512,321	A→T	11.4%	S35T (TCC→ACC)	aroK ←	shikimate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,512,321	0	A	T	11.4%	586.8 / 55.7	202	S35T (TCC→ACC)	aroK	shikimate kinase I
Reads supporting (aligned to +/- strand): ref base A (105/74); new base T (7/16); total (112/90)
Fisher's exact test for biased strand distribution p-value = 1.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 205 total aligned reads displayed.

AGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTC  >  CP009273/3512174‑3512468

agCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGa                                                                                                                                                     <  1:1749880/148‑1 (MQ=255)

  ccgccgccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTAt                                                                                                                                                   <  1:2500515/148‑1 (MQ=255)

  ccgccgccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTAt                                                                                                                                                   <  1:1930100/148‑1 (MQ=255)

  ccgccgccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAAt                                                                                                                                                   >  2:680964/1‑148 (MQ=255)

     ccgccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGt                                                                                                                                                >  2:359248/1‑148 (MQ=255)

         cAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTaaaa                                                                                                                                            <  2:1195953/148‑1 (MQ=255)

            tAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTc                                                                                                                                         <  2:1746601/148‑1 (MQ=255)

            tAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTc                                                                                                                                         <  1:1388552/148‑1 (MQ=255)

              gCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCa                                                                                                                                       <  1:557776/148‑1 (MQ=255)

                   cacaATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTg                                                                                                                                  <  1:3490708/148‑1 (MQ=255)

                     caATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAg                                                                                                                                <  2:2717806/148‑1 (MQ=255)

                       aTACCCTGTTTCTCGGTCAACTCAGTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCTGTTCGTTTCTCAATCACTTGATCGGAATCGTAAAATTCCATATTGAgtt                                                                                                                              <  1:1911739/148‑1 (MQ=255)

                        tACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAgttg                                                                                                                             <  1:1316182/148‑1 (MQ=255)

                        tACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAgttg                                                                                                                             <  2:2044837/148‑1 (MQ=255)

                        tACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATAGGAATCGTAAAATTCCATATTGAgttg                                                                                                                             >  1:3165575/1‑148 (MQ=255)

                          accggTTTCGCGGTCAACTCATTGTTGTCCATTTCTTGGCGTTCGGGGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAgttgt                                                                                                                            <  1:1511813/144‑1 (MQ=255)

                           ccTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGttgttg                                                                                                                          >  1:1046571/1‑148 (MQ=255)

                            cTGGTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGa                                                                                                                         <  2:781000/148‑1 (MQ=255)

                             tGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAg                                                                                                                        <  1:1268424/148‑1 (MQ=255)

                                 tctcGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTaa                                                                                                                    <  2:307541/148‑1 (MQ=255)

                                 tctcGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTaa                                                                                                                    <  2:2105193/148‑1 (MQ=255)

                                 tctcGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTaa                                                                                                                    <  2:1972530/148‑1 (MQ=255)

                                   tcGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACt                                                                                                                  >  2:3310297/1‑148 (MQ=255)

                                   tcGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACt                                                                                                                  <  1:1001445/148‑1 (MQ=255)

                                   tcGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACt                                                                                                                  >  2:1508137/1‑148 (MQ=255)

                                     ggTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTgg                                                                                                                >  1:2613604/1‑148 (MQ=255)

                                        cAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGgcgc                                                                                                             >  1:490681/1‑148 (MQ=255)

                                        cAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGgcgc                                                                                                             >  2:1383624/1‑148 (MQ=255)

                                          aCTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGccc                                                                                                           >  1:2995175/1‑148 (MQ=255)

                                             cATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCTAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTGGATGAGCTCACTGGCGCCAAAt                                                                                                        >  2:3011077/1‑148 (MQ=255)

                                             cATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCGCAAt                                                                                                        >  2:2123294/1‑148 (MQ=255)

                                             cATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAAt                                                                                                        >  2:1901880/1‑148 (MQ=255)

                                             cATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAAt                                                                                                        >  1:659262/1‑148 (MQ=255)

                                               ttGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAg                                                                                                      <  2:462943/148‑1 (MQ=255)

                                               ttGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAg                                                                                                      >  2:1488562/1‑148 (MQ=255)

                                                tgatgaCCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATCTTGAGTTGTTGAGCTAACTGGCGCCCAATAGt                                                                                                     >  2:2832099/1‑148 (MQ=255)

                                                tgatgaCCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGt                                                                                                     >  1:1032825/1‑148 (MQ=255)

                                                tgatgaCCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGt                                                                                                     >  1:2033252/1‑148 (MQ=255)

                                                  atgaCCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGc                                                                                                   >  2:1834936/1‑148 (MQ=255)

                                                  atgaCCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGc                                                                                                   >  2:3351448/1‑148 (MQ=255)

                                                    gaCCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCtt                                                                                                 >  1:180243/1‑148 (MQ=255)

                                                       cTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCttttt                                                                                              >  1:2683963/1‑148 (MQ=255)

                                                       cTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCttttt                                                                                              >  2:1041373/1‑148 (MQ=255)

                                                        ttttCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTc                                                                                             <  1:826567/148‑1 (MQ=255)

                                                        ttttCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTc                                                                                             <  2:3249980/148‑1 (MQ=255)

                                                         tttCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTcc                                                                                            >  1:2924818/1‑148 (MQ=255)

                                                                 cgATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCa                                                                                    >  2:1948123/1‑148 (MQ=255)

                                                                  gATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCAt                                                                                   <  1:2534082/148‑1 (MQ=255)

                                                                  gATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCAt                                                                                   <  2:381898/148‑1 (MQ=255)

                                                                   aTC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATa                                                                                  >  1:3242867/1‑148 (MQ=255)

                                                                   aTC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATa                                                                                  <  1:3091999/148‑1 (MQ=255)

                                                                     c‑gcgGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATCTTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAgg                                                                                >  2:2980800/1‑148 (MQ=255)

                                                                                ttcttcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGa                                                                      >  1:1130951/1‑148 (MQ=255)

                                                                                ttcttcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGa                                                                      >  2:2178336/1‑148 (MQ=255)

                                                                                  cttcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGaaa                                                                    >  2:2634378/1‑148 (MQ=255)

                                                                                   ttcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAg                                                                   >  2:91160/1‑148 (MQ=255)

                                                                                    tcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGa                                                                  >  2:614987/1‑148 (MQ=255)

                                                                                      gCCTTCTAAATCGAAAACCCAGCCCACCTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGata                                                                <  2:3108631/148‑1 (MQ=255)

                                                                                         ttCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTg                                                             <  1:88429/148‑1 (MQ=255)

                                                                                          tCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAg                                                                                                      <  2:1337767/106‑1 (MQ=255)

                                                                                          tCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAg                                                                                                      >  1:1337767/1‑106 (MQ=255)

                                                                                           cTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCg                                                           <  1:3540331/148‑1 (MQ=255)

                                                                                            tAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGt                                                          <  2:2392922/148‑1 (MQ=255)

                                                                                            tAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGt                                                          >  2:493543/1‑148 (MQ=255)

                                                                                                tCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTtct                                                      >  2:1033466/1‑148 (MQ=255)

                                                                                                tCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTtct                                                      >  1:2798077/1‑148 (MQ=255)

                                                                                                        ccAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAtt                                              >  1:2177371/1‑148 (MQ=255)

                                                                                                          aGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAtttt                                            <  1:97159/148‑1 (MQ=255)

                                                                                                           gCCCACCTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAttttt                                           <  2:2813278/148‑1 (MQ=255)

                                                                                                             ccACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTc                                         <  2:2194850/148‑1 (MQ=255)

                                                                                                                caTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGt                                      <  2:2033252/148‑1 (MQ=255)

                                                                                                                caTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGt                                      >  2:3121419/1‑148 (MQ=255)

                                                                                                                 ctCAGCTCCGGTTCGCTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGta                                     <  2:3475976/147‑1 (MQ=255)

                                                                                                                 aTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGta                                     <  1:3242560/148‑1 (MQ=255)

                                                                                                                   cAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGtact                                   >  1:94483/1‑148 (MQ=255)

                                                                                                                     gCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTAGGtactac                                 >  1:1240056/1‑148 (MQ=255)

                                                                                                                      cTCCGGTTCGTTACTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTactact                                <  2:1203464/148‑1 (MQ=255)

                                                                                                                       tCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTActacta                               <  2:2585678/148‑1 (MQ=255)

                                                                                                                           gTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGAc                           >  2:969425/1‑148 (MQ=255)

                                                                                                                            ttCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACt                          >  2:508113/1‑148 (MQ=255)

                                                                                                                               gTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTAtt                       <  2:2953545/148‑1 (MQ=255)

                                                                                                                                tttCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTc                      >  1:3156082/1‑148 (MQ=255)

                                                                                                                                 ttCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCg                     >  2:1831734/1‑148 (MQ=255)

                                                                                                                                  tctcAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGt                    >  1:2262845/1‑148 (MQ=255)

                                                                                                                                     cAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTaa                 <  1:1850243/148‑1 (MQ=255)

                                                                                                                                     cAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTaa                 <  1:1851688/148‑1 (MQ=255)

                                                                                                                                     cAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTaa                 <  1:2482127/148‑1 (MQ=255)

                                                                                                                                        tctcTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGa              <  2:1444887/148‑1 (MQ=255)

                                                                                                                                         ctctTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAATACTATTCGTTAATGAt             >  2:1001999/1‑148 (MQ=255)

                                                                                                                                          tctTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATa            <  2:500266/148‑1 (MQ=255)

                                                                                                                                          tctTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATa            <  1:359248/148‑1 (MQ=255)

                                                                                                                                           ctTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATaa           >  1:3436191/1‑148 (MQ=255)

                                                                                                                                             tGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAAc         >  2:2468028/1‑148 (MQ=255)

                                                                                                                                              gATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAAcc        <  2:1640048/148‑1 (MQ=255)

                                                                                                                                               aTCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAAccc       >  2:831581/1‑148 (MQ=255)

                                                                                                                                                tCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCg      >  1:3210402/1‑148 (MQ=255)

                                                                                                                                                tCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCg      >  1:1611311/1‑148 (MQ=255)

                                                                                                                                                 cGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGc     >  1:108143/1‑148 (MQ=255)

                                                                                                                                                 cGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGc     <  2:1294055/148‑1 (MQ=255)

                                                                                                                                                  ggAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCt    >  1:1692744/1‑148 (MQ=255)

AGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATC‑GCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTC  >  CP009273/3512174‑3512468

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