BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,720,597	Δ1 bp	coding (4/153 nt)	hokA ←	small toxic polypeptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,720,596	0	G	.	100.0%	202.4 / NA	42	coding (5/153 nt)	hokA	small toxic polypeptide
Reads supporting (aligned to +/- strand): ref base G (0/0); new base . (14/28); total (14/28)

AAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTG  >  NC_000913/3720549‑3720639

aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                           >  1:2180301/1‑48 (MQ=255)

aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgc                                            <  1:1946593/49‑2 (MQ=255)

   gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                       <  1:1907158/50‑1 (MQ=255)

   gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                       <  1:614710/50‑1 (MQ=255)

   gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                       <  1:346598/50‑1 (MQ=255)

     gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                      >  1:1096412/1‑49 (MQ=39)

      tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                    <  1:2068742/50‑1 (MQ=255)

      tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                    <  1:441513/50‑1 (MQ=255)

        aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                   <  1:2200185/49‑1 (MQ=255)

        aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                   <  1:349460/49‑1 (MQ=255)

         aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                 >  1:339217/1‑50 (MQ=255)

             aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                             >  1:1893311/1‑50 (MQ=255)

              aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                            >  1:2227819/1‑50 (MQ=255)

                 aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                         <  1:2091893/50‑1 (MQ=255)

                 aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                         <  1:1152617/50‑1 (MQ=255)

                 aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                         >  1:8742/1‑50 (MQ=255)

                  cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        <  1:672938/50‑1 (MQ=255)

                  cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        <  1:1130791/50‑1 (MQ=255)

                    aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                      >  1:1621554/1‑50 (MQ=255)

                    aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                      >  1:1458581/1‑50 (MQ=255)

                     ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                      <  1:1918847/49‑1 (MQ=255)

                     ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                      <  1:1144175/49‑1 (MQ=255)

                      tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                    >  1:1304171/1‑50 (MQ=255)

                       aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                    <  1:774781/49‑1 (MQ=255)

                       aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                    <  1:9254/49‑1 (MQ=255)

                       aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                    <  1:2040848/49‑1 (MQ=255)

                        aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                   <  1:2204146/49‑1 (MQ=255)

                        aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                   <  1:600365/49‑1 (MQ=255)

                        aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                   <  1:1420432/49‑1 (MQ=255)

                          gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                <  1:109944/50‑1 (MQ=255)

                          gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                <  1:1161514/50‑1 (MQ=255)

                           aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               >  1:2235997/1‑50 (MQ=255)

                             agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACg             <  1:1222590/50‑1 (MQ=255)

                             agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACg             <  1:597867/50‑1 (MQ=255)

                              gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGa            <  1:1598267/50‑1 (MQ=255)

                                   cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt        >  1:139160/1‑49 (MQ=255)

                                   cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt        >  1:229986/1‑49 (MQ=255)

                                   cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt        >  1:1832980/1‑49 (MQ=255)

                                      tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTGTaact    >  1:2102394‑M1/1‑46 (MQ=255)

                                      tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt        <  1:1641853/46‑1 (MQ=255)

                                      tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact    >  1:1046644‑M1/1‑46 (MQ=255)

                                      tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa       <  1:1827738‑M1/47‑2 (MQ=255)

                                        tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt  <  1:2230046‑M1/50‑7 (MQ=255)

                                        tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt  <  1:1202222‑M1/50‑7 (MQ=255)

AAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTG  >  NC_000913/3720549‑3720639

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