BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	583,341	T→G	94.7%	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	94.7%	998.0 / 22.9	261	S33R (AGT→CGT)	ybcH	PRK09936 family protein
Reads supporting (aligned to +/- strand): ref base T (7/7); new base G (114/132); total (122/139)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 261 total aligned reads displayed.

CAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTAC  >  CP009273/583194‑583487

 aGCTGCGGCCCGCTTAGACACTAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtc                                                                                                                                                 <  1:1441689/150‑1 (MQ=255)

 aGCTGCGGCCCGCTTAAACAATACCGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtc                                                                                                                                                 <  2:2532996/150‑1 (MQ=255)

     gCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCgg                                                                                                                                             <  1:1140288/150‑1 (MQ=255)

        gCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGttt                                                                                                                                          >  1:1014968/1‑150 (MQ=255)

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

         cccGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCAGTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGtttt                                                                                                                                         <  2:3043164/150‑1 (MQ=255)

            gCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTtgtg                                                                                                                                      <  2:2427325/150‑1 (MQ=255)

              ttAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGt                                                                                                                                    >  2:605102/1‑150 (MQ=255)

                aaaCAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTg                                                                                                                                  >  2:1782329/1‑150 (MQ=255)

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

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

                  aCAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGcc                                                                                                                                >  2:4183257/1‑150 (MQ=255)

                       aaCGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGata                                                                                                                           <  2:2379144/150‑1 (MQ=255)

                       aaCGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGata                                                                                                                           <  2:915605/150‑1 (MQ=255)

                         cGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataat                                                                                                                         >  2:1081016/1‑150 (MQ=255)

                          gTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataata                                                                                                                        >  2:925452/1‑150 (MQ=255)

                          gTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataata                                                                                                                        <  2:4514128/150‑1 (MQ=255)

                            gcgcTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATAcc                                                                                                                      <  1:2370925/150‑1 (MQ=255)

                            gcgcTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATAcc                                                                                                                      <  1:4709924/150‑1 (MQ=255)

                            gcgcTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATAcc                                                                                                                      <  1:4710201/150‑1 (MQ=255)

                             cgcTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCt                                                                                                                     <  2:2165549/150‑1 (MQ=255)

                                tGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTc                                                                                                                  >  2:2747399/1‑150 (MQ=255)

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

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

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

                                   tCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATc                                                                                                               >  2:4377732/1‑150 (MQ=255)

                                   tCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATc                                                                                                               <  2:3099566/150‑1 (MQ=255)

                                        ctggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCaaa                                                                                                          >  1:885264/1‑150 (MQ=255)

                                          gggTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGAGAATACCTTTCATCGCAAAgg                                                                                                        >  2:3354625/1‑150 (MQ=255)

                                           ggTAAATGCATCGCCGTAACGGGTCCAGTGCAAAACCAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGa                                                                                                       <  1:3971429/150‑1 (MQ=255)

                                            gTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaa                                                                                                      >  2:2955164/1‑150 (MQ=255)

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

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

                                             tAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaa                                                                                                     >  2:39148/1‑150 (MQ=255)

                                             tAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaa                                                                                                     >  1:3471515/1‑150 (MQ=255)

                                             tAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaa                                                                                                     >  1:3471906/1‑150 (MQ=255)

                                                  gCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAggg                                                                                                >  1:2514175/1‑150 (MQ=255)

                                                    aTCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCt                                                                                              >  1:1154160/1‑150 (MQ=255)

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

                                                     tCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTg                                                                                             <  1:4021891/150‑1 (MQ=255)

                                                     tCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTg                                                                                             <  2:3510677/150‑1 (MQ=255)

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

                                                      cGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGa                                                                                            >  1:1056912/1‑150 (MQ=255)

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

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

                                                       gCCGTACAGGGTCCATTGCAAAACAAGGGTATCGAAGCCGTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTGTTGTGGTTTCCCGATAATACCTTTCATCGCACAGGAAACAGGTCTGAc                                                                                           >  2:2807788/1‑150 (MQ=255)

                                                       gCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGAc                                                                                           <  1:720645/150‑1 (MQ=255)

                                                          gTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCaa                                                                                        >  1:3158962/1‑150 (MQ=255)

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

                                                            aaCGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCaaaa                                                                                      <  1:1623989/150‑1 (MQ=255)

                                                              cGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGc                                                                                    >  1:2110881/1‑150 (MQ=255)

                                                              cGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGc                                                                                    <  2:1889976/150‑1 (MQ=255)

                                                                ggTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAg                                                                                  <  1:2944513/150‑1 (MQ=255)

                                                                ggTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAg                                                                                  <  2:4407477/150‑1 (MQ=255)

                                                                ggTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAg                                                                                  <  2:1267539/150‑1 (MQ=255)

                                                                ggTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAg                                                                                  <  2:3254499/150‑1 (MQ=255)

                                                                ggTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAg                                                                                  <  1:4291349/150‑1 (MQ=255)

                                                                   ccATTGCAAAACAATGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAgtgt                                                                               <  2:989474/150‑1 (MQ=255)

                                                                   ccATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAgtgt                                                                               <  1:4044450/150‑1 (MQ=255)

                                                                   ccATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAgtgt                                                                               >  1:3268247/1‑150 (MQ=255)

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

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

                                                                       tGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGc                                                                           <  1:1759844/150‑1 (MQ=255)

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

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

                                                                             accaGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACg                                                                     <  1:3637513/147‑1 (MQ=255)

                                                                              cAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGa                                                                    >  2:4228976/1‑150 (MQ=255)

                                                                              cAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGa                                                                    >  1:473713/1‑150 (MQ=255)

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

                                                                                   gTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATg                                                               >  2:2871245/1‑150 (MQ=255)

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

                                                                                    tATCGACGCCTGGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGa                                                              <  1:388383/150‑1 (MQ=255)

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

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

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

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

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

                                                                                         aaGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTa                                                         <  2:534234/150‑1 (MQ=255)

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

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

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

                                                                                                     cGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAg                                             >  1:495557/1‑150 (MQ=255)

                                                                                                     cGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAg                                             <  1:3018402/150‑1 (MQ=255)

                                                                                                         aCTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCa                                         <  2:2616625/150‑1 (MQ=255)

                                                                                                          cTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCAc                                        <  2:2365702/150‑1 (MQ=255)

                                                                                                               tcatcaGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCCTCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCaa                                   <  2:1338271/150‑1 (MQ=255)

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

                                                                                                                       cccTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTgaga                           >  1:3357010/1‑150 (MQ=255)

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

                                                                                                                              aCTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgt                    >  1:874743/1‑150 (MQ=255)

                                                                                                                                tGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgttg                  <  2:473713/150‑1 (MQ=255)

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

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

                                                                                                                                          tAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCtctc        >  1:4463245/1‑150 (MQ=255)

                                                                                                                                           aaCCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCCATGGTGAGAAACGCGTTGTTGCGCTCTcc       >  2:1598467/1‑150 (MQ=255)

                                                                                                                                           aaCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTcc       <  2:557889/150‑1 (MQ=255)

                                                                                                                                           aaCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTcc       <  1:3732196/150‑1 (MQ=255)

                                                                                                                                            aCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCg      >  1:614441/1‑150 (MQ=255)

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

                                                                                                                                                gACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTAc  >  1:2429650/1‑150 (MQ=255)

CAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTAC  >  CP009273/583194‑583487

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