BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	256,882	C→T	7.2%	L262F (CTT→TTT)	proB →	gamma‑glutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	256,882	0	C	T	7.2%	971.1 / 38.4	277	L262F (CTT→TTT)	proB	gamma‑glutamate kinase
Reads supporting (aligned to +/- strand): ref base C (125/132); new base T (13/7); total (138/139)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

GCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAA  >  CP009273/256733‑257031

gcggcATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGc                                                                                                                                                       <  2:3866004/150‑1 (MQ=255)

gcggcATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGc                                                                                                                                                       <  1:2071523/150‑1 (MQ=255)

 cggcATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGGTCCATGCCCAGGCGACTCCGCt                                                                                                                                                      <  2:169640/150‑1 (MQ=255)

 cggcATGAGTACCAAAGTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCt                                                                                                                                                      <  2:2484976/150‑1 (MQ=255)

   gcATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTg                                                                                                                                                    >  2:3482651/1‑150 (MQ=255)

   gcATGAGTACCAAAGTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTg                                                                                                                                                    <  2:3446943/150‑1 (MQ=255)

    cATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGa                                                                                                                                                   <  2:2146712/150‑1 (MQ=255)

    cATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGa                                                                                                                                                   <  1:482156/150‑1 (MQ=255)

      tGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACAACATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTACGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGaaa                                                                                                                                                 >  2:1253595/1‑150 (MQ=255)

        aGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAAc                                                                                                                                               <  1:2826906/150‑1 (MQ=255)

           aCCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGt                                                                                                                                            <  1:224286/150‑1 (MQ=255)

            ccAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTa                                                                                                                                           <  1:2149436/150‑1 (MQ=255)

             cAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTaa                                                                                                                                          >  2:1329737/1‑150 (MQ=255)

             cAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTaa                                                                                                                                          <  1:2043610/150‑1 (MQ=255)

               aaTTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAAc                                                                                                                                        >  2:3597583/1‑150 (MQ=255)

               aaTTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAAc                                                                                                                                        >  1:1098905/1‑150 (MQ=255)

                aTTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACg                                                                                                                                       <  2:4463488/150‑1 (MQ=255)

                   gCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTg                                                                                                                                    >  2:3036305/1‑150 (MQ=255)

                     aGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGa                                                                                                                                  >  2:3916057/1‑150 (MQ=255)

                     aGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGa                                                                                                                                  >  1:3754451/1‑150 (MQ=255)

                       gCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGTTTGAAAACCGTAAACGCTGGAtt                                                                                                                                <  1:2752214/150‑1 (MQ=255)

                         cGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGTTTGAAAACCGTAAACGCTGGAtttt                                                                                                                              <  2:2788180/150‑1 (MQ=255)

                         cGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGAtttt                                                                                                                              >  1:3189443/1‑150 (MQ=255)

                         cGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGAtttt                                                                                                                              <  1:3980804/150‑1 (MQ=255)

                         cGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGAtttt                                                                                                                              <  1:3872731/150‑1 (MQ=255)

                           cTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCg                                                                                                                            <  2:1478771/150‑1 (MQ=255)

                               catGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTgc                                                                                                                        <  2:1151449/148‑1 (MQ=255)

                                  ggCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGcgcc                                                                                                                     >  2:3385956/1‑150 (MQ=255)

                                  ggCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGcgcc                                                                                                                     >  2:1006368/1‑150 (MQ=255)

                                     ttGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCgccgcc                                                                                                                  >  2:536983/1‑150 (MQ=255)

                                      tGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCGTTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGccgccg                                                                                                                 >  1:1388883/1‑150 (MQ=255)

                                      tGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGccgccg                                                                                                                 <  1:3361564/150‑1 (MQ=255)

                                        ccGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCcggc                                                                                                               <  2:4171670/150‑1 (MQ=255)

                                         cGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCcggcg                                                                                                              <  2:4058315/150‑1 (MQ=255)

                                           tGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCggg                                                                                                            <  1:2474020/150‑1 (MQ=255)

                                             cGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTg                                                                                                          <  2:2669158/150‑1 (MQ=255)

                                             cGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTg                                                                                                          <  2:336065/150‑1 (MQ=255)

                                               ggTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGaa                                                                                                        <  1:712322/150‑1 (MQ=255)

                                                gTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGaaa                                                                                                       >  2:144420/1‑150 (MQ=255)

                                                  aTCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGTTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATc                                                                                                     >  2:534838/1‑150 (MQ=255)

                                                  aTCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATc                                                                                                     <  2:2725473/150‑1 (MQ=255)

                                                  aTCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGGTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATc                                                                                                     >  2:925732/1‑150 (MQ=255)

                                                   tCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGTTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCa                                                                                                    <  2:4419022/150‑1 (MQ=255)

                                                       cacCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGt                                                                                                >  2:4376923/1‑150 (MQ=255)

                                                        acCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGTAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTa                                                                                               <  1:3482651/150‑1 (MQ=255)

                                                           attattGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGAt                                                                                            >  1:3713628/1‑150 (MQ=255)

                                                             tattGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGa                                                                                          <  1:3572480/150‑1 (MQ=255)

                                                             tattGCCGCGGGCAGCAAGCCGGGCGGTATTGGTTATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGa                                                                                          >  1:1140260/1‑150 (MQ=255)

                                                                tGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCATGGTAGATGAAgg                                                                                       <  2:2963480/150‑1 (MQ=255)

                                                                tGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAgg                                                                                       <  1:1525176/150‑1 (MQ=255)

                                                                 gCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGTAGGCATTTACGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAATGCGTTAACGCTGGATTTTCGGTGCGCCGCCGGAGGGTGAATTCACTGTAGATGTAggg                                                                                      <  2:821696/150‑1 (MQ=255)

                                                                 gCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAggg                                                                                      <  1:2289027/150‑1 (MQ=255)

                                                                     cgGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCaa                                                                                  >  2:3006979/1‑150 (MQ=255)

                                                                     cgGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCaa                                                                                  >  2:686694/1‑150 (MQ=255)

                                                                     cgGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCaa                                                                                  <  2:3544613/150‑1 (MQ=255)

                                                                       ggCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGTTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACt                                                                                >  1:463909/1‑150 (MQ=255)

                                                                       ggCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCCACt                                                                                >  1:3222003/1‑150 (MQ=255)

                                                                          agcaAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTgcc                                                                             >  1:4008804/1‑150 (MQ=255)

                                                                            caAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTgccgc                                                                           >  2:1781297/1‑150 (MQ=255)

                                                                             aaGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGTTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTgccgcc                                                                          >  2:1085921/1‑150 (MQ=255)

                                                                             aaGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTgccgcc                                                                          >  2:2353162/1‑150 (MQ=255)

                                                                               gCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCAt                                                                        <  2:1969486/150‑1 (MQ=255)

                                                                                ccGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCAtt                                                                       <  1:4104146/150‑1 (MQ=255)

                                                                                 cGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTc                                                                      >  1:591037/1‑150 (MQ=255)

                                                                                  gggCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCt                                                                     <  2:562430/150‑1 (MQ=255)

                                                                                      gTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGaa                                                                 >  1:1211455/1‑150 (MQ=255)

                                                                                      gTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGaa                                                                 >  1:181075/1‑150 (MQ=255)

                                                                                       ttATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGTTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAAc                                                                >  2:327981/1‑150 (MQ=255)

                                                                                         aTTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAAcgc                                                              >  2:3202205/1‑150 (MQ=255)

                                                                                             gtgatgtgatGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCa                                                          >  1:1456671/1‑150 (MQ=255)

                                                                                             gtgatgtgatGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCa                                                          <  1:3386672/150‑1 (MQ=255)

                                                                                              tgatgtgatgGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAg                                                         >  2:3022206/1‑150 (MQ=255)

                                                                                               gatgtgatgGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGc                                                        >  1:1250642/1‑150 (MQ=255)

                                                                                                   tgatgGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTccc                                                    >  2:1675161/1‑150 (MQ=255)

                                                                                                   tgatgGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTccc                                                    >  1:750234/1‑150 (MQ=255)

                                                                                                       ggAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGtt                                                <  2:2904455/150‑1 (MQ=255)

                                                                                                       ggAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGtt                                                >  2:4334070/1‑150 (MQ=255)

                                                                                                       ggAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGtt                                                <  1:1604440/150‑1 (MQ=255)

                                                                                                           ggCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCg                                            >  2:2310539/1‑150 (MQ=255)

                                                                                                             cATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGaa                                          >  1:83541/1‑150 (MQ=255)

                                                                                                               tttCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAg                                        >  1:1518955/1‑150 (MQ=255)

                                                                                                                ttCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAgg                                       >  2:239872/1‑150 (MQ=255)

                                                                                                                ttCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAgg                                       <  2:293284/150‑1 (MQ=255)

                                                                                                                   cGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCAt                                    <  2:1742875/150‑1 (MQ=255)

                                                                                                                     tCGGTACGCTGTTCCATGCCCAGGCGACTCCGTTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTa                                  <  2:282699/150‑1 (MQ=255)

                                                                                                                       ggTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTaaa                                >  1:3945300/1‑150 (MQ=255)

                                                                                                                        gTACGCTGTTCCATGCCCAGGCGACTCCGTTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTaaaa                               <  2:2113393/150‑1 (MQ=255)

                                                                                                                            gCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGt                           <  1:4238383/150‑1 (MQ=255)

                                                                                                                               gTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGAc                        <  2:2018786/150‑1 (MQ=255)

                                                                                                                                ttCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACt                       <  1:432582/150‑1 (MQ=255)

                                                                                                                                ttCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACt                       <  1:3312102/150‑1 (MQ=255)

                                                                                                                                  ccATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTgg                     >  2:4611586/1‑150 (MQ=255)

                                                                                                                                       cccAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAAtt                >  2:3996052/1‑150 (MQ=255)

                                                                                                                                          aGGCTACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTtct             >  1:3652887/1‑150 (MQ=255)

                                                                                                                                           ggCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTtctc            <  2:4099536/150‑1 (MQ=255)

                                                                                                                                           ggCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTtctc            >  2:2575625/1‑150 (MQ=255)

                                                                                                                                           ggCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTtctc            >  1:3253187/1‑150 (MQ=255)

                                                                                                                                            gCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTcg           >  1:3017718/1‑150 (MQ=255)

                                                                                                                                             cGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTcgc          >  2:593849/1‑150 (MQ=255)

                                                                                                                                                cTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCgtg       >  1:2604341/1‑150 (MQ=255)

GCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGTATCGACACCATTATTGCCGCGGGCAGCAAGCCGGGCGTTATTGGTGATGTGATGGAAGGCATTTCCGTCGGTACGCTGTTCCATGCCCAGGCGACTCCGCTTGAAAACCGTAAACGCTGGATTTTCGGTGCGCCGCCGGCGGGTGAAATCACGGTAGATGAAGGGGCAACTGCCGCCATTCTGGAACGCGGCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAA  >  CP009273/256733‑257031

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