BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	583,341	T→G	93.8%	S33R (AGT→CGT)	ybcH ←	PRK09936 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	583,341	0	T	G	93.8%	1050.9 / 31.4	275	S33R (AGT→CGT)	ybcH	PRK09936 family protein
Reads supporting (aligned to +/- strand): ref base T (7/10); new base G (124/134); total (131/144)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 276 total aligned reads displayed.

CAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTT  >  CP009273/583194‑583490

cAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAt                                                                                                                                                     <  1:2880591/150‑1 (MQ=255)

 aGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtc                                                                                                                                                    >  2:3313258/1‑150 (MQ=255)

 aGCTGCGGCCCGCTTAAACAAGAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtc                                                                                                                                                    <  2:4211154/150‑1 (MQ=255)

   cTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtctc                                                                                                                                                  <  1:923229/150‑1 (MQ=255)

    tGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCg                                                                                                                                                 <  1:3136800/150‑1 (MQ=255)

      cGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGt                                                                                                                                               >  2:4250227/1‑150 (MQ=255)

      cGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCAACAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGt                                                                                                                                               >  1:2668417/1‑150 (MQ=255)

       ggCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGtt                                                                                                                                              <  1:2161033/150‑1 (MQ=255)

         cccGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGtttt                                                                                                                                            >  2:3242193/1‑150 (MQ=255)

         cccGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGtttt                                                                                                                                            >  2:4052167/1‑150 (MQ=255)

          ccGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTtg                                                                                                                                           >  1:1769222/1‑150 (MQ=255)

          ccGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTtg                                                                                                                                           >  1:3850775/1‑150 (MQ=255)

            gCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTtgtg                                                                                                                                         >  2:252963/1‑150 (MQ=255)

             cTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTgg                                                                                                                                        <  1:1599989/150‑1 (MQ=255)

               tAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGtt                                                                                                                                      <  2:2639983/150‑1 (MQ=255)

                aaaCAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTg                                                                                                                                     >  1:2836065/1‑150 (MQ=255)

                aaaCAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTg                                                                                                                                     <  2:821367/150‑1 (MQ=255)

                   cAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCa                                                                                                                                  <  1:2052742/150‑1 (MQ=255)

                    aaTAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAg                                                                                                                                 <  2:1771486/150‑1 (MQ=255)

                         cGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataat                                                                                                                            <  2:57949/150‑1 (MQ=255)

                              gcTGTTCTGGCTGGGTAACTGCATCGCCGTAACGGGTCCATTGCAAAACAAGGTTATCGAAGCCTTTCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCtt                                                                                                                       <  1:1587477/150‑1 (MQ=255)

                                tGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTc                                                                                                                     <  1:4365743/150‑1 (MQ=255)

                                tGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTc                                                                                                                     <  2:477867/150‑1 (MQ=255)

                                 gTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCa                                                                                                                    <  1:1757159/150‑1 (MQ=255)

                                 attCTGGCTGGGTAAATGCATCGCCGTAACGGGTCAATTGAAAAACAAGGGTATCGAAGACTTGCAAACGTAACTGACTCATCAGCACCTGAAACTGGGTATCGGTAAACTGACGAACTATGTTATGTGGTTGACAGATAAAACCtataa                                                                                                                    >  2:3455873/2‑146 (MQ=255)

                                  ttCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCAt                                                                                                                   <  2:3147964/150‑1 (MQ=255)

                                          gggTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAgg                                                                                                           <  1:1461530/150‑1 (MQ=255)

                                           ggTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGa                                                                                                          <  2:2422589/150‑1 (MQ=255)

                                            gTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaa                                                                                                         <  2:1093336/150‑1 (MQ=255)

                                             tAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaa                                                                                                        >  1:208114/1‑150 (MQ=255)

                                             tAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaa                                                                                                        <  2:2947375/150‑1 (MQ=255)

                                                 tGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAgg                                                                                                    >  2:1409975/1‑150 (MQ=255)

                                                 tGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAgg                                                                                                    >  2:3911909/1‑150 (MQ=255)

                                                  gCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAggg                                                                                                   <  1:1817661/150‑1 (MQ=255)

                                                   cATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAATGGAAAAAGGGc                                                                                                  >  1:1872497/1‑150 (MQ=255)

                                                    aTCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCt                                                                                                 <  1:4346403/150‑1 (MQ=255)

                                                      cGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGa                                                                                               <  2:1271602/150‑1 (MQ=255)

                                                      cGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGa                                                                                               <  2:2667420/150‑1 (MQ=255)

                                                         cGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCa                                                                                            <  2:2976277/150‑1 (MQ=255)

                                                           tAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCaaa                                                                                          >  1:2367015/1‑150 (MQ=255)

                                                           tAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCaaa                                                                                          >  2:2182260/1‑150 (MQ=255)

                                                           tAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCaaa                                                                                          >  1:734672/1‑150 (MQ=255)

                                                               gggTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCa                                                                                      >  2:2566830/1‑150 (MQ=255)

                                                                 gTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAgt                                                                                    <  2:130868/150‑1 (MQ=255)

                                                                 gTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAgt                                                                                    <  1:3938544/150‑1 (MQ=255)

                                                                 gTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAgt                                                                                    <  2:4461372/150‑1 (MQ=255)

                                                                    cATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTc                                                                                 <  1:899930/150‑1 (MQ=255)

                                                                     aTTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCa                                                                                >  1:1257394/1‑150 (MQ=255)

                                                                     aTTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCa                                                                                >  1:377848/1‑150 (MQ=255)

                                                                        gCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCa                                                                             >  1:177920/1‑150 (MQ=255)

                                                                        gCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCa                                                                             <  2:3110611/150‑1 (MQ=255)

                                                                        gCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCa                                                                             <  2:2252873/150‑1 (MQ=255)

                                                                         cAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCaa                                                                            >  2:2105873/1‑150 (MQ=255)

                                                                          aaaaCAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAAt                                                                           <  2:88400/150‑1 (MQ=255)

                                                                           aaaCAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATa                                                                          <  1:2302262/150‑1 (MQ=255)

                                                                            aaCAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATAc                                                                         >  2:2176934/1‑150 (MQ=255)

                                                                               aaGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGaa                                                                      >  2:3546562/1‑150 (MQ=255)

                                                                                aGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGaaa                                                                     <  1:1010280/150‑1 (MQ=255)

                                                                                   gTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATg                                                                  <  2:1610155/150‑1 (MQ=255)

                                                                                    tATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGa                                                                 >  2:1956236/1‑150 (MQ=255)

                                                                                      tCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAAc                                                               <  1:4537149/150‑1 (MQ=255)

                                                                                      tCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGGTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAACAGCAGTGTCAGCAATACGAAAATGAAc                                                               <  1:2261323/150‑1 (MQ=255)

                                                                                       cGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAATGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACt                                                              >  1:3090380/1‑150 (MQ=255)

                                                                                        gAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACtt                                                             >  1:3011458/1‑150 (MQ=255)

                                                                                         aaGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTa                                                            >  2:2474245/1‑150 (MQ=255)

                                                                                          aGCCTTGCTAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTAc                                                           <  1:1324874/150‑1 (MQ=255)

                                                                                          aGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTAc                                                           <  1:584573/150‑1 (MQ=255)

                                                                                          aGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTAc                                                           <  2:2497806/150‑1 (MQ=255)

                                                                                            ccTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGc                                                         >  2:4433377/1‑150 (MQ=255)

                                                                                              ttGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAt                                                       <  1:2458487/150‑1 (MQ=255)

                                                                                               tGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAtt                                                      <  2:2769658/150‑1 (MQ=255)

                                                                                                gCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAttt                                                     >  1:4509765/1‑150 (MQ=255)

                                                                                                 cAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTa                                                    >  2:1052231/1‑150 (MQ=255)

                                                                                                    aCGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCa                                                 >  1:4456906/1‑150 (MQ=255)

                                                                                                      gTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGt                                               >  1:2410728/1‑150 (MQ=255)

                                                                                                        aaCTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGc                                             >  1:2580034/1‑150 (MQ=255)

                                                                                                        aaCTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGc                                             <  2:3645310/150‑1 (MQ=255)

                                                                                                        aaCTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGc                                             >  1:3233676/1‑150 (MQ=255)

                                                                                                         aCTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCa                                            >  2:2226082/1‑150 (MQ=255)

                                                                                                             aCTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTcc                                        >  2:1663421/1‑150 (MQ=255)

                                                                                                             aCTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTcc                                        >  2:3349587/1‑150 (MQ=255)

                                                                                                             aCTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTcc                                        >  2:3871877/1‑150 (MQ=255)

                                                                                                                catcaGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCact                                     >  2:4163308/1‑148 (MQ=255)

                                                                                                                    aGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTg                                 >  2:2364921/1‑150 (MQ=255)

                                                                                                                        ccTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGaa                             <  2:3606305/150‑1 (MQ=255)

                                                                                                                        ccTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGaa                             <  2:648178/150‑1 (MQ=255)

                                                                                                                          tGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAc                           <  2:3850969/150‑1 (MQ=255)

                                                                                                                          tGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAc                           >  1:3660459/1‑150 (MQ=255)

                                                                                                                             cACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAcgcg                        >  1:3437385/1‑150 (MQ=255)

                                                                                                                               cTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgtt                      <  2:340741/150‑1 (MQ=255)

                                                                                                                                 gggTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgttgt                    >  2:1552663/1‑150 (MQ=255)

                                                                                                                                 gggTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgttgt                    <  1:3463713/150‑1 (MQ=255)

                                                                                                                                  ggTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgttgtt                   <  1:3966231/150‑1 (MQ=255)

                                                                                                                                  ggTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgttgtt                   <  1:2979857/150‑1 (MQ=255)

                                                                                                                                     aTCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTgcg                <  2:2730063/150‑1 (MQ=255)

                                                                                                                                             ccTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGt        >  2:3842306/1‑150 (MQ=255)

                                                                                                                                              cTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGtt       >  1:2543025/1‑150 (MQ=255)

                                                                                                                                              cTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGtt       <  2:3731229/150‑1 (MQ=255)

                                                                                                                                               tGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTa      <  1:4052167/150‑1 (MQ=255)

                                                                                                                                                   tATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGtt  <  2:1646377/150‑1 (MQ=255)

CAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTT  >  CP009273/583194‑583490

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