BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,182,778	A→G	12.2%	N12N (AAT→AAC)	sprL ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,182,778	0	A	G	12.2%	92.6 / 6.7	41	N12N (AAT→AAC)	sprL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (18/18); new base G (3/2); total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

CTCCATAATGCTCATATTGTTTTGGATTAATTTCAATATCATGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAATACTGTAAATCCAAACGTCAACGATATCACCAACACTGACAATATCCATTGGATT  >  CP000730/2182642‑2182913

ctcCATAATGCTCATATTGTTTTGGATTAATTTCAATATCATGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATtct                                                                                                                                     >  2:91864/1‑141 (MQ=255)

     tAATGCTCATATTGTTTTGGATTAATTTCAATATCATGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAAcc                                                                                                                                <  1:264567/141‑1 (MQ=255)

          cTCATATTGTTTTGGATTAATTTCAATATCATGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCg                                                                                                                           <  1:382337/141‑1 (MQ=255)

            cATATTGTTTTGGATTAATTTCAATATCATGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCt                                                                                                                         >  1:428393/1‑141 (MQ=255)

               atTGTTTTGGATTAATTTCAATATCATGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGGTTCTCAACCATTCGCTGaa                                                                                                                      >  2:29888/1‑141 (MQ=255)

                ttGTTTTGGATTAATTTCAATATCATGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAg                                                                                                                     <  2:82986/141‑1 (MQ=255)

                 tGTTTTGGATTAATTTCAATATCATGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGt                                                                                                                    <  2:407549/141‑1 (MQ=255)

                    tttGGATTAATTTCAATATCATGTGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTg                                                                                                                 <  1:374887/141‑1 (MQ=255)

                        gATTAATTTCAATATCATGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCat                                                                                                             <  2:126437/141‑1 (MQ=255)

                                 cAATATCATGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCAtgt                                                                                                    >  2:163582/1‑141 (MQ=255)

                                        aTGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGatcaa                                                                                             >  1:396344/1‑141 (MQ=255)

                                             aCTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGGTTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAatcatc                                                                                        >  1:43393/1‑141 (MQ=255)

                                             aCTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGCAGTTTGTCATTATTCATGTGGATCAatcatc                                                                                        >  1:100683/1‑141 (MQ=255)

                                                ttAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGtt                                                                                     <  1:288785/141‑1 (MQ=255)

                                                 tAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTa                                                                                    >  1:240347/1‑141 (MQ=255)

                                                  aaaaGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAatcat                                                                                         >  2:90142/1‑135 (MQ=255)

                                                  aaaaGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAatcat                                                                                         <  1:90142/135‑1 (MQ=255)

                                                     aGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGGTTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGa                                                                                >  2:102220/1‑141 (MQ=255)

                                                      gATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAACTATGCACAATGTCGAAACATACGTCCAAATTTTTCTTCCGACAGATTCCCAACCATTCGCTGAAGTTTGTCATTATTCATGTCGATCAATCATCGTTAATGAt                                                                               <  2:208370/141‑1 (MQ=255)

                                                                gCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCtt                                                                     <  1:337882/141‑1 (MQ=255)

                                                                   tgttgtACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGGTTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTAt                                                                  <  2:43393/141‑1 (MQ=255)

                                                                        tACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTg                                                             >  1:403907/1‑141 (MQ=255)

                                                                         aCGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGt                                                            >  1:34804/1‑141 (MQ=255)

                                                                           gTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCa                                                          <  1:23054/141‑1 (MQ=255)

                                                                            tAACCTTTTATTAAAATATGCACAATGTCTAACCGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTATGTcta                                                         <  2:394616/141‑3 (MQ=255)

                                                                             aaCCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGGTTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAAt                                                        <  1:102220/141‑1 (MQ=255)

                                                                                   ttattaAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAATACTGTa                                                  >  1:434419/1‑141 (MQ=255)

                                                                                     attaAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCACTACTGTaaa                                                >  2:80254/1‑141 (MQ=255)

                                                                                     attaAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAATACTGTaaa                                                >  2:66709/1‑141 (MQ=255)

                                                                                             atGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAATACTGTAAATCCAAACg                                        >  1:126374/1‑141 (MQ=255)

                                                                                                cacaATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCTATACTGTAAATCCAAACGTCa                                     >  2:163347/1‑141 (MQ=255)

                                                                                                 acaATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCACTACTGTAAATACAAACGTCaa                                    >  1:282225/1‑141 (MQ=255)

                                                                                                    aTGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGatcaatc                                                                                           <  1:189422/83‑1 (MQ=255)

                                                                                                    aTGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGatcaatc                                                                                           >  2:189422/1‑83 (MQ=255)

                                                                                                                   tCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAATACTGTAAATCCAAACGTCAACGATATCACCAACACTGa                  <  2:253215/141‑1 (MQ=255)

                                                                                                                          ttttCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAATACTGTAAATCCAAACGTCAACGATATCACCAACACTGACAATATc           <  2:257297/141‑1 (MQ=255)

                                                                                                                            ttcttcCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAATACTGTAAATCCAAACGTCAACGATATCACCAACACTGACAATATCCa         >  1:74805/1‑141 (MQ=255)

                                                                                                                               ttcCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTTTCAATACTGTAAATCCAAACGTCAACGATATCACCAACACTGACAATATCCATTg      >  1:88242/1‑141 (MQ=255)

                                                                                                                                tcCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAATACTGTAAATCCAAACGTCAACGATATCACCAACACTGACAATATCCATTgg     <  2:9256/141‑1 (MQ=255)

                                                                                                                                  cGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAATACTGTAAATCCAAACGTCAACGATATCACCAACACTGACAATATCCATTGGAt   <  1:49414/141‑1 (MQ=255)

                                                                                                                                   gAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAATACTGTAAATCCAAACGTCAACGATATCACCAACACTGACAATATCCATTGGAtt  <  1:91197/141‑1 (MQ=255)

CTCCATAATGCTCATATTGTTTTGGATTAATTTCAATATCATGGGACTTTAAAAGATAACGTCCGCCTGTTGTACGTAACCTTTTATTAAAATATGCACAATGTCGAAACGTACGTCCAAATTTTTCTTCCGAAAGATTCTCAACCATTCGCTGAAGTTTGTCATTATTCATGTGGATCAATCATCGTTAATGATACTTTGTCTTTATTTTTGTCAATACTGTAAATCCAAACGTCAACGATATCACCAACACTGACAATATCCATTGGATT  >  CP000730/2182642‑2182913

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