BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,913	C→G	23.6%	S213S (TCG→TCC)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	35,913	0	C	G	23.6%	623.8 / 314.1	465	S213S (TCG→TCC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (163/191); new base G (69/41); total (232/233)
Fisher's exact test for biased strand distribution p-value = 3.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

Only 100 of 471 total aligned reads displayed.

AGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTC  >  CP000730/35783‑36053

  cTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTCACTTGCTAGCCTGATGCTAATTTCGTGGCATGGCGAAAGCCCGTAg                                                                                                                                  >  2:131953/1‑141 (MQ=255)

   taaaataaaaTTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCCCTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGa                                                                                                                                 >  1:18956/1‑141 (MQ=255)

         aaaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga                                                                                                                         <  2:855521/129‑1 (MQ=11)

           aagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga                                                                                                                       <  2:20577/131‑1 (MQ=11)

            agagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                                                      <  1:359737/132‑1 (MQ=11)

              tgtgATCTAATAAAACTAGAAGGTTCTGTTGCAAAGTAAAGAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGGTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                                                      <  1:929180/141‑1 (MQ=255)

                tattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                                                                                  <  2:915771/136‑1 (MQ=11)

                 attttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGCTCTGAAGAGACCTCCg                                                                                                                 >  2:427998/5‑141 (MQ=11)

                 attttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                                                                                 >  1:415276/5‑141 (MQ=11)

                  ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                                                                <  1:810969/138‑1 (MQ=11)

                  ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATAGAACTGATCGCTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                                                                >  1:250794/4‑141 (MQ=11)

                    ttaataaAAAT‑GATGGTACTGTTGCAAAGT‑AAAAAATATAGCTAATCATTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                                                                                                              <  2:618766/140‑1 (MQ=11)

                     taataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAGGAGACCTGCGGttt                                                                                                             >  1:223600/1‑140 (MQ=11)

                      aataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGAATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt                                                                                                              >  2:405487/1‑141 (MQ=255)

                      aataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                                                                                            >  1:43771/1‑141 (MQ=11)

                        aaaaaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                                                                                            <  2:370507/140‑1 (MQ=37)

                                 aaGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGGATGGCGCAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtatata                                                                                                   <  1:727516/141‑1 (MQ=11)

                                                                 aatcACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTa                                                                   <  2:699055/138‑1 (MQ=0)

                                                                           ttgtCATGTCAGTGGTCGCTTAACGTGCTAGCATGATGCTAATTGCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATAGAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTg                                                         <  2:321688/138‑1 (MQ=1)

                                                                                    ctggtgtcgaatAACGTGCTCGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGATACCTGCGGTTCTTTGTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATAc                                                <  2:340785/130‑1 (MQ=11)

                                                                                     aGTGTTCGTTAAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt                                               >  2:717704/1‑141 (MQ=11)

                                                                                     aGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGAATTCTTTGCTGTATTAATACt                                               >  1:285578/1‑141 (MQ=9)

                                                                                      gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt                                              <  2:598443/141‑1 (MQ=17)

                                                                                       tgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACttt                                             >  1:323907/1‑141 (MQ=17)

                                                                                       tgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACttt                                             >  1:942660/1‑141 (MQ=32)

                                                                                       tgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACttt                                             >  1:129961/1‑141 (MQ=32)

                                                                                        gtgcGCTTAACTGGCTATGATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCAGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTg                                            <  2:145444/138‑1 (MQ=11)

                                                                                         ttCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGa                                           <  2:1056391/141‑1 (MQ=255)

                                                                                         ttCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGa                                           <  2:925738/141‑1 (MQ=17)

                                                                                         ttCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAAAACGTAGATCGGAAGATACCTGCGGTTCTTTTTATCCGGTGCTGAGAGACACTCAATACCTTTTAAAGTATTCTTTGCGGTATTGATACTTTGa                                           >  2:1018132/1‑141 (MQ=11)

                                                                                          tCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGGGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAt                                          <  2:380904/141‑1 (MQ=17)

                                                                                          tCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAt                                          <  2:538413/141‑1 (MQ=17)

                                                                                           cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAACTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTCAAGTATTCTTTGCTGTATTAATACTTTGAtc                                         >  2:895251/4‑140 (MQ=11)

                                                                                           cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGAGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATa                                         >  1:681934/4‑141 (MQ=11)

                                                                                              ttAACTTGGTAGCATGATGCTAATTTCGGGGCATGGCGAAAGTCCGTGGATCTGAAGAGACCTGCGGTTCTTCTTACATAGAGCGTAAATACACTCAATACCTTTTAAAGGATTCTTTGCTGTATTGATACTTTGATACCt                                      <  1:917566/141‑1 (MQ=11)

                                                                                              ttAACTTGCTGGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCt                                      <  2:1039738/141‑1 (MQ=17)

                                                                                              ttAACTTGCTAGGATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCt                                      <  1:33614/141‑1 (MQ=17)

                                                                                              ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCt                                      <  1:271953/141‑1 (MQ=17)

                                                                                              ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCGGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTACAGTATTCTTTGCTGTATTGATACTTTGATACCt                                      >  1:66975/1‑141 (MQ=12)

                                                                                              ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGAGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCt                                      <  1:492147/141‑1 (MQ=17)

                                                                                                aaCTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTg                                    >  2:471446/1‑141 (MQ=255)

                                                                                                aaCTTGCTAGCATGATGCTAATTTCGTGGCATAGCGAAAATCCGTAGATCTGAAGAGACCTGAGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCttt                                    >  1:760882/1‑140 (MQ=12)

                                                                                                aaCTTGCTAGCATGAGGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTg                                    <  1:565980/141‑1 (MQ=17)

                                                                                                 acctgcatgCATGATGCTAATTTCGGGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATAGAGAGCTTAAATACACTCAATACTTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGt                                   <  2:951918/133‑1 (MQ=11)

                                                                                                 aCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGt                                   >  1:1054000/1‑141 (MQ=17)

                                                                                                  cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtc                                  <  2:1002896/141‑1 (MQ=17)

                                                                                                     gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGAGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtcttt                               >  1:796362/1‑141 (MQ=12)

                                                                                                     gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGCTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtcttt                               <  1:632953/141‑1 (MQ=17)

                                                                                                      ttAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                              <  1:973948/140‑1 (MQ=17)

                                                                                                      cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                              <  2:676415/141‑1 (MQ=17)

                                                                                                      cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTGTTATATAGAGCGTAAATACACGCCATACCTTGTATACTATTCTTTGCTGTATTAATACTTTGATAACTTGtctttc                              >  1:128322/1‑141 (MQ=11)

                                                                                                       tAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATGGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctttct                             <  2:879468/141‑1 (MQ=11)

                                                                                                       gaGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGTGGCGTGCGGTTCTTGTTAGATGGAGCGTAAAGACACGCTATACCGTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                             <  1:740455/140‑1 (MQ=11)

                                                                                                       caGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGTTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctttct                             <  2:31074/140‑1 (MQ=11)

                                                                                                        aggatgatgCCAGTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttctt                            <  2:239836/138‑1 (MQ=255)

                                                                                                        aGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttctt                            <  1:152995/141‑1 (MQ=17)

                                                                                                        aGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttctt                            >  1:855064/1‑141 (MQ=17)

                                                                                                        aGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttctt                            <  1:475996/141‑1 (MQ=17)

                                                                                                         gCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTa                           >  1:536909/1‑141 (MQ=17)

                                                                                                         gCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATAGATACTTTGATACCTTGTCTTTCTTa                           >  1:1007847/1‑141 (MQ=17)

                                                                                                           atgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACt                         >  2:22010/1‑141 (MQ=255)

                                                                                                           atgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACt                         >  1:301526/1‑141 (MQ=21)

                                                                                                            tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACtt                        >  2:639170/1‑141 (MQ=255)

                                                                                                            tgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATGCCTTGTCTTTCTTACtt                        <  1:841948/141‑1 (MQ=25)

                                                                                                            tgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACtt                        >  2:774625/1‑141 (MQ=25)

                                                                                                             gatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACttt                       <  1:754474/141‑1 (MQ=25)

                                                                                                                 cTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAata                   >  1:955381/1‑141 (MQ=25)

                                                                                                                  tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCCATACCTCTTAAAGTATTACTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTTCTTTAatat                  >  2:444662/1‑141 (MQ=11)

                                                                                                                  tAATTTCGTGGCATGGCGAAAATCCGAAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatat                  >  2:658058/1‑141 (MQ=21)

                                                                                                                  tAAGTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTCGCAATAATGATACTTTGATACCTTGTCTTTCTTACTTTAatat                  >  2:480315/1‑141 (MQ=11)

                                                                                                                   aaTTTCGTGGTATGTGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATg                 <  1:400762/141‑1 (MQ=11)

                                                                                                                   aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGAGGTTCTTTTTTTTTATAGCGTAAATACTCTCAATACCTTTTAGGTTTTTCGGTGCTGTATGTATACTTTGATACCTTGTCTTTCTTACTTTAATATg                 <  2:1047565/141‑1 (MQ=11)

                                                                                                                     tttCGTGGCATGGCGACAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATGGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATgga               >  2:461125/1‑139 (MQ=255)

                                                                                                                      ttCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGTTACTTTAATATCTTGTCTTTCTTACTTTAATATGACg              >  1:656687/1‑141 (MQ=11)

                                                                                                                      ttCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACg              >  2:288700/1‑141 (MQ=25)

                                                                                                                      ttCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTATTTGCTGTATTGATCATTAGATACATTTTTTTTCTTACTTTAATATTACg              >  2:740689/1‑141 (MQ=11)

                                                                                                                      ttCGTGGCATGGCGAAAAGCCGTAGATCTGAAGAGACCGGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTGATCTTTCTTACTTTAATATGACg              >  2:547753/1‑141 (MQ=11)

                                                                                                                       tCGTGGCATGGCGAATATCCGTAGATGTGAAGAGTCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTATAGTATTCTTTGCTCTATTGATACTTTGATACCTTGTCCTTCTTACTTTAATATGACgg             <  1:646095/141‑1 (MQ=11)

                                                                                                                       tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg             <  2:668033/141‑1 (MQ=37)

                                                                                                                        cGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGt            <  1:613596/141‑1 (MQ=14)

                                                                                                                        cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGt            <  2:980025/141‑1 (MQ=38)

                                                                                                                        cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGt            >  1:1022064/1‑141 (MQ=37)

                                                                                                                        cGTGGCATGGCGAAAATCCGGAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGt            <  2:824567/141‑1 (MQ=255)

                                                                                                                         gTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTg           >  2:751407/1‑141 (MQ=11)

                                                                                                                         gTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTg           >  2:523339/1‑141 (MQ=11)

                                                                                                                          tGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATGTTGGCTTTCTTACTTTAATATAACGGTGa          >  1:604642/1‑141 (MQ=11)

                                                                                                                           ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAt         <  2:503491/141‑1 (MQ=39)

                                                                                                                           ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAt         >  2:18479/1‑141 (MQ=37)

                                                                                                                            gCATGGCGAAAATCCGTAGATCTGAAGAGCCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATc        <  1:949489/141‑1 (MQ=255)

                                                                                                                            gCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATc        >  2:258182/1‑141 (MQ=38)

                                                                                                                             cATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGTTACATTGTCTTTCTTACTTTAATATGACGGTGATCt       >  2:166561/1‑141 (MQ=11)

                                                                                                                             cATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCGTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTCCTTACTTTAATATGACGGTGATCt       <  1:549364/141‑1 (MQ=255)

                                                                                                                              aTGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCtt      >  2:400792/1‑141 (MQ=11)

                                                                                                                               tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATGCCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg     >  1:448449/1‑141 (MQ=255)

                                                                                                                               tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg     <  1:885406/141‑1 (MQ=38)

                                                                                                                                ggCGAAAATCCGTAGATCTCAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGTGGGATCTATACTTTGCTCCCTCTTCGTTCTTACTTTAATATGACGGTGATCTTGc    <  1:509919/141‑1 (MQ=14)

                                                                                                                                 gCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCt   <  1:580170/141‑1 (MQ=255)

                                                                                                                                 gCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCt   <  1:89529/141‑1 (MQ=255)

                                                                                                                                  cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTc  >  1:467058/1‑141 (MQ=37)

                                                                                                                                  cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAAGACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTGTTTCTTACTTTAATATGACGGTGATCTTGCta  >  1:710373/1‑140 (MQ=255)

AGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTC  >  CP000730/35783‑36053

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`