BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	109,050	C→T	6.4%	E13K (GAA→AAA)	zapD ←	FtsZ stabilizer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	109,050	0	C	T	6.4%	1030.9 / 40.4	297	E13K (GAA→AAA)	zapD	FtsZ stabilizer
Reads supporting (aligned to +/- strand): ref base C (141/137); new base T (14/5); total (155/142)
Fisher's exact test for biased strand distribution p-value = 5.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 300 total aligned reads displayed.

CTCAGTGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCATC  >  CP009273/108901‑109199

 tCAGTGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCa                                                                                                                                                      >  2:866952/1‑150 (MQ=255)

  cAGTGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCAt                                                                                                                                                     >  1:1870779/1‑150 (MQ=255)

   agagCGGACTTCGCCGCGCTCGAAAACATACAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATATTTTCAtt                                                                                                                                                    >  2:3660821/4‑150 (MQ=255)

    gTGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCAttt                                                                                                                                                   <  2:4186248/150‑1 (MQ=255)

     tGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTa                                                                                                                                                  <  2:1882921/150‑1 (MQ=255)

      gCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAg                                                                                                                                                 <  2:3916165/150‑1 (MQ=255)

        ggACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTg                                                                                                                                               <  2:899194/150‑1 (MQ=255)

        ggACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTg                                                                                                                                               >  2:1787673/1‑150 (MQ=255)

         taCTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTgg                                                                                                                                              <  2:3385829/149‑1 (MQ=255)

            ttCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATg                                                                                                                                           >  1:4366254/1‑150 (MQ=255)

            ttCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATg                                                                                                                                           <  2:2173296/150‑1 (MQ=255)

               gccgcGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTc                                                                                                                                        <  2:796911/150‑1 (MQ=255)

                   cgcTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCaaaa                                                                                                                                    >  1:3242293/1‑150 (MQ=255)

                   cgcTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCaaaa                                                                                                                                    <  1:2240205/150‑1 (MQ=255)

                   cgcTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCaaaa                                                                                                                                    <  1:1379538/150‑1 (MQ=255)

                     cTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAg                                                                                                                                  >  2:2969627/1‑150 (MQ=255)

                      tCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAgg                                                                                                                                 <  2:3578327/150‑1 (MQ=255)

                      tCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAgg                                                                                                                                 >  2:2332135/1‑150 (MQ=255)

                        gAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGAc                                                                                                                               >  2:3727103/1‑150 (MQ=255)

                          aaaCATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCt                                                                                                                             >  1:2889333/1‑150 (MQ=255)

                           aaCATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTg                                                                                                                            >  2:3601762/1‑150 (MQ=255)

                           aaCATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTg                                                                                                                            <  2:3063642/150‑1 (MQ=255)

                               tCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTc                                                                                                                        <  2:4268365/150‑1 (MQ=255)

                               tCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTc                                                                                                                        >  1:1255065/1‑150 (MQ=255)

                                  aGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGc                                                                                                                     >  2:3465853/1‑150 (MQ=255)

                                      aTTCACTGACATTTCGGAAGACATGCAGCGCGCCAGCGTGGTCAACGATGGGTAATTTATCGGTGAGTTGCTGAATCAAAATCTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTTGATGTTCAAAAAGGACCTGGGTCTGCATTa                                                                                                                 <  1:1200337/150‑1 (MQ=255)

                                        tCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGGTCAAAAAGGACCTGGGTCTGCATTACg                                                                                                               >  2:4111735/1‑150 (MQ=255)

                                         cACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGGGCAGCCATGTACGCGTTTTTTCATTTAGTGGATGTTCAAAAAGGTCCTGGGTCTGCATTACgg                                                                                                              <  2:2908365/150‑1 (MQ=255)

                                            tGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGttt                                                                                                           <  1:2720240/150‑1 (MQ=255)

                                            tGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGttt                                                                                                           >  2:2642339/1‑150 (MQ=255)

                                               cATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTc                                                                                                        >  1:3989034/1‑150 (MQ=255)

                                                 ttACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCt                                                                                                      >  1:3735783/1‑150 (MQ=255)

                                                   aCGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCtgt                                                                                                    >  1:4530818/1‑150 (MQ=255)

                                                   aCGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCtgt                                                                                                    <  1:1835345/150‑1 (MQ=255)

                                                     ggAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTga                                                                                                  <  1:419035/150‑1 (MQ=255)

                                                        agaaATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGAc                                                                                               <  1:3967847/150‑1 (MQ=255)

                                                         taaaTGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACa                                                                                              <  2:3973964/149‑1 (MQ=255)

                                                            aTGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAAc                                                                                           <  2:718279/150‑1 (MQ=255)

                                                            aTGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAAc                                                                                           >  1:4007192/1‑150 (MQ=255)

                                                               cAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGc                                                                                        >  1:3930450/1‑150 (MQ=255)

                                                                aGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCg                                                                                       <  1:2332135/150‑1 (MQ=255)

                                                                 gcgcgcCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGa                                                                                      <  1:2255517/150‑1 (MQ=255)

                                                                  cgcgcCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGAc                                                                                     <  2:3609570/150‑1 (MQ=255)

                                                                   gcgcCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACg                                                                                    <  1:3479222/150‑1 (MQ=255)

                                                                     gcCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCa                                                                                  <  1:381570/150‑1 (MQ=255)

                                                                       cAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAg                                                                                >  1:284779/1‑150 (MQ=255)

                                                                        aGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGc                                                                               >  2:174547/1‑150 (MQ=255)

                                                                        aGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGc                                                                               >  2:2215801/1‑150 (MQ=255)

                                                                         gCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCt                                                                              <  1:2969627/150‑1 (MQ=255)

                                                                          cGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTg                                                                             >  2:1027048/1‑150 (MQ=255)

                                                                          cGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTg                                                                             >  2:3305473/1‑150 (MQ=255)

                                                                          cGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTg                                                                             <  1:1482674/150‑1 (MQ=255)

                                                                           gTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGc                                                                            >  1:1046757/1‑150 (MQ=255)

                                                                           gTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGc                                                                            <  2:3704389/150‑1 (MQ=255)

                                                                              gTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGaaaa                                                                         >  1:3339479/1‑150 (MQ=255)

                                                                                 aaCGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAg                                                                      >  2:1933925/1‑150 (MQ=255)

                                                                                  aCGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAgt                                                                     >  1:3958193/1‑150 (MQ=255)

                                                                                    gATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAgtgt                                                                   <  2:1638742/150‑1 (MQ=255)

                                                                                        ggTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTgtgcgt                                                               >  1:165051/1‑150 (MQ=255)

                                                                                            aaTTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAg                                                           <  1:2197401/150‑1 (MQ=255)

                                                                                                aaCGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCgg                                                       >  2:3018146/1‑150 (MQ=255)

                                                                                                  cGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGc                                                     >  2:510526/1‑150 (MQ=255)

                                                                                                       aGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAGTGGCTGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACAt                                                <  1:2601792/150‑1 (MQ=255)

                                                                                                        gTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATc                                               >  1:1736488/1‑150 (MQ=255)

                                                                                                        gTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATc                                               <  1:592727/150‑1 (MQ=255)

                                                                                                         ttGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATcc                                              >  2:2408964/1‑150 (MQ=255)

                                                                                                          tGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCg                                             >  2:1215438/1‑150 (MQ=255)

                                                                                                             tGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACAACCGATg                                          >  2:4089262/1‑150 (MQ=255)

                                                                                                               aaTCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCg                                        >  1:3090298/1‑150 (MQ=255)

                                                                                                               aaTCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCg                                        >  1:3485614/1‑150 (MQ=255)

                                                                                                                aTCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCTACATCCGATGCGa                                       <  2:1870779/150‑1 (MQ=255)

                                                                                                                 tCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGAt                                      <  2:1000381/150‑1 (MQ=255)

                                                                                                                 tCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGAt                                      >  2:3719150/1‑150 (MQ=255)

                                                                                                                  cAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATa                                     <  1:2652318/150‑1 (MQ=255)

                                                                                                                  cAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATa                                     <  1:3438348/150‑1 (MQ=255)

                                                                                                                  cAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATa                                     <  1:354813/150‑1 (MQ=255)

                                                                                                                    aaaaCTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGc                                   >  1:363969/1‑150 (MQ=255)

                                                                                                                      aaCTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCAt                                 <  2:3920105/150‑1 (MQ=255)

                                                                                                                        cTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATcc                               >  1:3071332/1‑150 (MQ=255)

                                                                                                                         tCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATTTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCATGTGAGACAAACTGCGACGCAAGCTGAAAATAGTGTGCGTGAAGGCGGGCAACATCAGATGCGATAGCATcag                              >  2:2758567/1‑148 (MQ=255)

                                                                                                                         tCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCg                              >  2:2758654/1‑150 (MQ=255)

                                                                                                                           aaTGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGt                            <  2:4481525/150‑1 (MQ=255)

                                                                                                                           aaTGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGt                            >  2:3162429/1‑150 (MQ=255)

                                                                                                                               cgcAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTgcgc                        >  2:4437344/1‑150 (MQ=255)

                                                                                                                                  aGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGt                     >  1:2043419/1‑150 (MQ=255)

                                                                                                                                     cATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGttat                  >  2:4218679/1‑150 (MQ=255)

                                                                                                                                     cATGTACGCATTTTTTTATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGttat                  >  2:2577856/1‑150 (MQ=255)

                                                                                                                                         tACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATc              <  2:287732/150‑1 (MQ=255)

                                                                                                                                             cATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATg          >  1:679627/1‑150 (MQ=255)

                                                                                                                                             cATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATGACCGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATg          >  1:3703008/1‑150 (MQ=255)

                                                                                                                                              aTTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGa         >  1:3763441/1‑150 (MQ=255)

                                                                                                                                               ttttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGAc        <  1:3727102/150‑1 (MQ=255)

                                                                                                                                               ttttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGAc        <  1:2240537/150‑1 (MQ=255)

                                                                                                                                               ttttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGAc        >  1:2932837/1‑150 (MQ=255)

                                                                                                                                                tttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACg       >  1:3859531/1‑150 (MQ=255)

                                                                                                                                                 ttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGt      >  1:2211028/1‑150 (MQ=255)

                                                                                                                                                 ttttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGt      >  2:2018549/1‑150 (MQ=255)

                                                                                                                                                  tttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTc     <  1:244193/150‑1 (MQ=255)

                                                                                                                                                  tttCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTc     >  2:2257176/1‑150 (MQ=255)

                                                                                                                                                     cATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCATc  >  1:932506/1‑150 (MQ=255)

CTCAGTGCGGACTTCGCCGCGCTCGAAAACATCCAGTAATTCACTGACATTACGGAAGAAATGCAGCGCGCCAGCGTGGTCAACGATGGGTAAATTAACGGTGAGTTGCTGAATCAAAAACTCAATGCGCAGCCATGTACGCATTTTTTCATTTAGTGGATGTTCAAAAAGGACCTGGGTCTGCATTACGGTTTTTCCTGTGAGACAAACTGCGACGCAAGCTGCAAATAGTGTGCGTGCAGGCGGGCAACATCCGATGCGATAGCATCCGGTGCGCCGTTATTATCAATGACGTCATC  >  CP009273/108901‑109199

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