BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,922,164	C→A	47.2%	P311P (CCG→CCT)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,922,164	0	C	A	47.2%	157.6 / inf	252	P311P (CCG→CCT)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base C (62/71); new base A (47/72); total (109/143)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 255 total aligned reads displayed.

ATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAA  >  CP009273/1922017‑1922308

aTAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGc                                                                                                                                                  >  1:340260/1‑148 (MQ=255)

 tAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCg                                                                                                                                                 >  2:3651757/1‑148 (MQ=255)

   ggCCGGATCATCACAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGa                                                                                                                                               <  1:2066111/148‑1 (MQ=255)

   ggCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGa                                                                                                                                               >  2:808621/1‑148 (MQ=255)

   ggCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGa                                                                                                                                               <  2:345709/148‑1 (MQ=255)

   ggCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGa                                                                                                                                               <  1:135664/148‑1 (MQ=255)

    gCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGaa                                                                                                                                              >  1:1013004/1‑148 (MQ=255)

    gCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGaa                                                                                                                                              <  2:75392/148‑1 (MQ=255)

    gCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGaa                                                                                                                                              <  1:1323858/148‑1 (MQ=255)

     ccGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAAt                                                                                                                                             <  2:1193883/148‑1 (MQ=255)

        gATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTaa                                                                                                                                          <  1:2227924/148‑1 (MQ=255)

           catcaAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACGGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAActc                                                                                                                                       >  1:785573/1‑148 (MQ=255)

             tcaAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCtt                                                                                                                                     >  2:1720940/1‑148 (MQ=255)

              caAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTc                                                                                                                                    <  2:1797911/148‑1 (MQ=255)

               aaaGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTcc                                                                                                                                   <  1:2449293/148‑1 (MQ=255)

                aaGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTccc                                                                                                                                  >  1:389132/1‑148 (MQ=255)

                    cAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGACTTAACTCTTCCCACTg                                                                                                                              >  1:3525606/1‑148 (MQ=255)

                    cAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTg                                                                                                                              >  1:2937969/1‑148 (MQ=255)

                    cAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTg                                                                                                                              <  2:2212308/148‑1 (MQ=255)

                     aaTACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGt                                                                                                                             <  1:770050/148‑1 (MQ=255)

                     aaTACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGt                                                                                                                             <  1:3341087/148‑1 (MQ=255)

                     aaTACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGt                                                                                                                             <  2:2624254/148‑1 (MQ=255)

                          cAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCtc                                                                                                                        >  2:918550/1‑148 (MQ=255)

                           aaTGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCtca                                                                                                                       >  2:1271708/1‑148 (MQ=255)

                            aTGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGAATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCtcat                                                                                                                      >  1:1640052/1‑148 (MQ=255)

                               aCTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCAc                                                                                                                   >  1:1612813/1‑148 (MQ=255)

                                 ttCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGc                                                                                                                 >  2:2764449/1‑148 (MQ=255)

                                 ttCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGc                                                                                                                 <  1:1478851/148‑1 (MQ=255)

                                              cGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGt                                                                                                    <  2:3093661/148‑1 (MQ=255)

                                               ggTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGta                                                                                                   <  2:3259994/148‑1 (MQ=255)

                                                gTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGtat                                                                                                  >  1:1952545/1‑148 (MQ=255)

                                                 ttAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGtata                                                                                                 >  1:3003636/1‑148 (MQ=255)

                                                   aaTTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAg                                                                                               <  2:2083279/148‑1 (MQ=255)

                                                      ttGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCa                                                                                            >  2:1941525/1‑148 (MQ=255)

                                                       tGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCaa                                                                                           <  2:1239909/148‑1 (MQ=255)

                                                       tGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCaa                                                                                           <  1:1709804/148‑1 (MQ=255)

                                                       tGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTACCACCAAGATGACTTCTCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCaa                                                                                           <  1:3745644/148‑1 (MQ=255)

                                                            cAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGttt                                                                                      <  1:1282523/148‑1 (MQ=255)

                                                            cAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGttt                                                                                      <  2:2792788/148‑1 (MQ=255)

                                                            cAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGttt                                                                                      <  2:1413752/148‑1 (MQ=255)

                                                               aCTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTg                                                                                   <  2:1303800/148‑1 (MQ=255)

                                                               aCTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGACTTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTg                                                                                   >  2:2472772/1‑148 (MQ=255)

                                                                cTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGc                                                                                  >  1:1424033/1‑148 (MQ=255)

                                                                 tGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGATTTTGcc                                                                                 >  1:2305926/1‑148 (MQ=255)

                                                                     tAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTgg                                                                             >  2:281886/1‑148 (MQ=255)

                                                                        cccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCgg                                                                          <  1:1271892/148‑1 (MQ=255)

                                                                          cgcgCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGtt                                                                        <  1:1938031/148‑1 (MQ=255)

                                                                             gcTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGa                                                                     <  1:2909835/148‑1 (MQ=255)

                                                                                 aCGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGc                                                                 <  2:1612813/148‑1 (MQ=255)

                                                                                    ctctTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGa                                                              >  2:1730230/1‑148 (MQ=255)

                                                                                      ctTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATa                                                            >  2:906179/1‑148 (MQ=255)

                                                                                       ttCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaa                                                           >  1:3758808/1‑148 (MQ=255)

                                                                                        tCAACCACCCGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaa                                                          <  2:2966470/148‑1 (MQ=255)

                                                                                        tCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaa                                                          >  1:3034760/1‑148 (MQ=255)

                                                                                         cAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaa                                                         <  1:1159028/148‑1 (MQ=255)

                                                                                          aaCCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaaa                                                        <  2:51117/148‑1 (MQ=255)

                                                                                          aaCCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaaa                                                        >  2:2595872/1‑148 (MQ=255)

                                                                                              accaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGAt                                                    >  1:1003891/1‑148 (MQ=255)

                                                                                              accaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGAt                                                    >  1:2129482/1‑148 (MQ=255)

                                                                                              accaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGAt                                                    <  2:2660776/148‑1 (MQ=255)

                                                                                              accaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGAt                                                    >  1:2386980/1‑148 (MQ=255)

                                                                                                cagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGc                                                  <  2:3231875/148‑1 (MQ=255)

                                                                                                  cccagTCGGTAAACCGCTTAAACCATTCCAGCATGATGTTTTCGGGTGGAGGAATTAACTCGTCCCACTGTTGCTCATCACGGTTCGGGGTACGGTATACGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTg                                                <  1:3272031/147‑1 (MQ=255)

                                                                                                      gTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAg                                            <  2:1546596/148‑1 (MQ=255)

                                                                                                        cGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCTGATAAAAACGATGCTGGTAGTg                                          >  2:2301897/1‑148 (MQ=255)

                                                                                                        cGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTg                                          <  1:1633535/148‑1 (MQ=255)

                                                                                                               cAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAgg                                   >  1:2437771/1‑148 (MQ=255)

                                                                                                                aGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGc                                  <  1:2741269/148‑1 (MQ=255)

                                                                                                                   gTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGt                               <  1:3391290/148‑1 (MQ=255)

                                                                                                                   gTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGt                               <  2:478867/148‑1 (MQ=255)

                                                                                                                    tAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTa                              <  1:2673465/148‑1 (MQ=255)

                                                                                                                    tAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTa                              <  1:3089807/148‑1 (MQ=255)

                                                                                                                     aaaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTAt                             >  2:3059152/1‑148 (MQ=255)

                                                                                                                     aaaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTAt                             >  2:440011/1‑148 (MQ=255)

                                                                                                                     aaaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTAt                             >  2:2365000/1‑148 (MQ=255)

                                                                                                                     aaaCCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTAt                             >  2:3391139/1‑148 (MQ=255)

                                                                                                                     aaaCCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTAt                             >  2:3391241/1‑148 (MQ=255)

                                                                                                                      aaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTAtt                            >  2:3777231/1‑148 (MQ=255)

                                                                                                                      aaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTAtt                            >  2:1887782/1‑148 (MQ=255)

                                                                                                                       aCCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTc                           <  1:1871042/148‑1 (MQ=255)

                                                                                                                        cccTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCg                          >  2:2054377/1‑148 (MQ=255)

                                                                                                                         ccTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGt                         <  1:2010760/148‑1 (MQ=255)

                                                                                                                          cTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTg                        >  1:1149007/1‑148 (MQ=255)

                                                                                                                           ttCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGa                       >  2:153414/1‑148 (MQ=255)

                                                                                                                           ttCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGa                       <  2:1493662/148‑1 (MQ=255)

                                                                                                                           ttCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGa                       >  1:3704044/1‑148 (MQ=255)

                                                                                                                            tCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGAt                      <  2:1648972/148‑1 (MQ=255)

                                                                                                                            tCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGAt                      <  2:1297633/148‑1 (MQ=255)

                                                                                                                            tCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGAt                      <  1:1576724/148‑1 (MQ=255)

                                                                                                                              cAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCt                    <  2:339928/148‑1 (MQ=255)

                                                                                                                               aGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCtt                   <  2:1873171/148‑1 (MQ=255)

                                                                                                                               aGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCtt                   <  1:3794316/148‑1 (MQ=255)

                                                                                                                               aGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTAGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCtt                   <  1:505761/148‑1 (MQ=255)

                                                                                                                                   tgatgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcg               >  1:809193/1‑148 (MQ=255)

                                                                                                                                    gatgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgtgg              >  1:998426/1‑148 (MQ=255)

                                                                                                                                      tgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcg            >  2:954585/1‑148 (MQ=255)

                                                                                                                                      tgTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcg            <  1:918550/148‑1 (MQ=255)

                                                                                                                                        ttttCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcg          >  1:3381338/1‑148 (MQ=255)

                                                                                                                                                tGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGaaaa  >  2:424661/1‑148 (MQ=255)

                                                                                                                                                tGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGAGGCGCGGCAGaaaa  >  2:3531116/1‑148 (MQ=255)

ATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAA  >  CP009273/1922017‑1922308

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