BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	583,341	T→G	83.6%	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	83.6%	979.3 / 134.1	330	S33R (AGT→CGT)	ybcH	PRK09936 family protein
Reads supporting (aligned to +/- strand): ref base T (21/34); new base G (123/152); total (144/186)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 331 total aligned reads displayed.

CAGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTT  >  CP009273/583194‑583490

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

 aGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtc                                                                                                                                                    <  2:4995910/150‑1 (MQ=255)

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

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

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

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

         cccGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGtttt                                                                                                                                            <  2:4638150/150‑1 (MQ=255)

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

           cgcgTAAACACTCCCGGGCGCTGATCTGGCTGGGTAACTGCATCGCCGTAACGGGTCCATGGCAAACCACGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACGGGGTATCGGTAACCGGACTATCTCGGTTTtgt                                                                                                                                          <  2:2901884/150‑1 (MQ=255)

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

               tAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGtt                                                                                                                                      >  1:4840750/1‑150 (MQ=255)

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

                 aaCAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGc                                                                                                                                    >  2:3718690/1‑150 (MQ=255)

                     aTAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGa                                                                                                                                <  2:4013235/150‑1 (MQ=255)

                      tAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGat                                                                                                                               >  1:3454899/1‑150 (MQ=255)

                       aaCGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGata                                                                                                                              >  2:5618889/1‑150 (MQ=255)

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

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

                            gcgcTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATAcc                                                                                                                         >  2:1042846/1‑150 (MQ=255)

                             cgcTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCt                                                                                                                        >  1:1039326/1‑150 (MQ=255)

                              gcTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCtt                                                                                                                       >  2:1893999/1‑150 (MQ=255)

                                    ctggctggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCg                                                                                                                 <  2:3950503/150‑1 (MQ=255)

                                    ctggctggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCg                                                                                                                 >  2:298118/1‑150 (MQ=255)

                                        ctggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCaaa                                                                                                             <  1:4570946/150‑1 (MQ=255)

                                         tggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAg                                                                                                            >  1:4351660/1‑150 (MQ=255)

                                         tggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAg                                                                                                            <  2:2644460/150‑1 (MQ=255)

                                         tggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAg                                                                                                            <  1:4766206/150‑1 (MQ=255)

                                         tggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAg                                                                                                            >  2:2328759/1‑150 (MQ=255)

                                         tggGTAAATGCATCGCCGTAACGGGTACATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAg                                                                                                            <  1:3049238/150‑1 (MQ=255)

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

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

                                                   cATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGc                                                                                                  <  1:1090468/150‑1 (MQ=255)

                                                   cATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGc                                                                                                  >  1:2329605/1‑150 (MQ=255)

                                                    aTCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCt                                                                                                 <  2:836562/150‑1 (MQ=255)

                                                        ccGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGAcc                                                                                             <  2:2505855/150‑1 (MQ=255)

                                                        ccGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGAcc                                                                                             <  1:377142/150‑1 (MQ=255)

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

                                                           tAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGATTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGATGACCaaa                                                                                          >  2:5553968/1‑150 (MQ=255)

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

                                                             aCGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGGAAAAGGGCTGACCAAAAg                                                                                        <  2:232355/150‑1 (MQ=255)

                                                             aCGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAg                                                                                        >  2:1161806/1‑150 (MQ=255)

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

                                                                 gTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAgt                                                                                    >  1:2023038/1‑150 (MQ=255)

                                                                  tCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAgtg                                                                                   <  2:3547991/150‑1 (MQ=255)

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

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

                                                                     aTTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCa                                                                                <  1:2634500/150‑1 (MQ=255)

                                                                     aTTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCa                                                                                <  1:5358823/150‑1 (MQ=255)

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

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

                                                                           aaaCAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATa                                                                          <  1:2287220/150‑1 (MQ=255)

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

                                                                             aCAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACg                                                                        >  2:3494531/1‑150 (MQ=255)

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

                                                                                     aTCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGaa                                                                <  2:4091981/150‑1 (MQ=255)

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

                                                                                        gAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACtt                                                             <  1:2490085/150‑1 (MQ=255)

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

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

                                                                                         aaGCCTTGCAAACGTCACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTa                                                            <  2:4663202/150‑1 (MQ=255)

                                                                                          aGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTAc                                                           <  1:2376126/150‑1 (MQ=255)

                                                                                           gCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACg                                                          <  1:2664854/150‑1 (MQ=255)

                                                                                            ccTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGAATATCTAGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAAAGAACTAACGc                                                         >  1:146089/1‑150 (MQ=255)

                                                                                              ttGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAt                                                       >  2:2061097/1‑150 (MQ=255)

                                                                                                gCAAACGTAACTGACTCATCAGCCCCTGCCACTGTGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAttt                                                     <  2:295294/150‑1 (MQ=255)

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

                                                                                                 cAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTa                                                    <  2:121768/150‑1 (MQ=255)

                                                                                                  aaaCGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTAc                                                   >  1:3613411/1‑150 (MQ=255)

                                                                                                   aaCGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTAcc                                                  <  1:3524256/150‑1 (MQ=255)

                                                                                                   aaCGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTAcc                                                  <  1:229955/150‑1 (MQ=255)

                                                                                                      gTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGt                                               <  2:1248747/150‑1 (MQ=255)

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

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

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

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

                                                                                                               tcatcaGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCaa                                      >  2:1139808/1‑150 (MQ=255)

                                                                                                                  tcaGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATgg                                   >  1:156915/1‑150 (MQ=255)

                                                                                                                     gCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTga                                >  2:5191601/1‑150 (MQ=255)

                                                                                                                       cccTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTgaga                              <  1:4915121/150‑1 (MQ=255)

                                                                                                                           gCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACTAGTGCACTCCAATGGTGAGAAAcg                          <  1:2713208/150‑1 (MQ=255)

                                                                                                                               cTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgtt                      <  1:4208076/150‑1 (MQ=255)

                                                                                                                               cTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgtt                      >  2:4246748/1‑150 (MQ=255)

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

                                                                                                                                tGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgttg                     >  1:3015838/1‑150 (MQ=255)

                                                                                                                                tGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgttg                     >  2:1110501/1‑150 (MQ=255)

                                                                                                                                  ggTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgttgtt                   <  1:3877496/150‑1 (MQ=255)

                                                                                                                                   gTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGttgttg                  >  2:4588396/1‑150 (MQ=255)

                                                                                                                                    tATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTgc                 >  1:718061/1‑150 (MQ=255)

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

                                                                                                                                        ggTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGctc             <  1:2246388/150‑1 (MQ=255)

                                                                                                                                         gTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGctct            <  2:3019939/150‑1 (MQ=255)

                                                                                                                                          tAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCtctc           <  2:1990484/150‑1 (MQ=255)

                                                                                                                                           aaCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCACTACGAAACTGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTcc          <  2:5492980/150‑1 (MQ=255)

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

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

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

                                                                                                                                            aCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCg         <  1:1152924/150‑1 (MQ=255)

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

                                                                                                                                                gACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTAc     <  1:2821218/150‑1 (MQ=255)

                                                                                                                                                   tATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGtt  <  1:86697/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: `‑`