D66 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,136,229	A→G	A116A (GCT→GCC)	SAUSA300_RS05595 ←	ribonuclease HIII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,136,229	0	A	G	98.4%	158.5 / NA	62	A116A (GCT→GCC)	SAUSA300_RS05595	ribonuclease HIII
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (1/60); total (1/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TTCTACAATTTTGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCAGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGTTAAATTGATCATACATTAAAGTTTGTTCAAG  >  NC_007793/1136142‑1136321

ttCTACAATTTTGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTc                                                                                  <  1:1446777/100‑1 (MQ=255)

ttCTACAATTTTGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTc                                                                                  <  1:5425082/100‑1 (MQ=255)

 tCTACAATTTTGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCa                                                                                 <  1:1544721/100‑1 (MQ=255)

 tCTACAATTTTGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCa                                                                                 <  1:3653814/100‑1 (MQ=255)

 tCTACAATTTTGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCa                                                                                 <  1:1795831/100‑1 (MQ=255)

 tCTACAATTTTGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCGTTGGTAACGAAAGCGGCACACACAGTCa                                                                                 <  1:4591272/100‑1 (MQ=255)

     cAATTTTGGTATCGGTTAATTTTTTTGAATCATCTACTCCCGAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCaagga                                                                             <  1:3699029/100‑4 (MQ=255)

        ttttGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGAcc                                                                          <  1:3155/100‑1 (MQ=255)

         tttGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCa                                                                         <  1:5687459/100‑1 (MQ=255)

           tGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACggggctgg                                                                       <  1:2988639/100‑7 (MQ=255)

           tGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACAGatcgga                                                                       >  1:1520051/1‑95 (MQ=255)

                      aaTTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCag                                                            <  1:1836025/100‑2 (MQ=255)

                          tttttGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTcc                                                        <  1:4422553/100‑1 (MQ=255)

                             ttGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAga                                                     <  1:844078/100‑2 (MQ=255)

                              cttccgatcTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGca                                                    <  1:2198519/94‑2 (MQ=255)

                               ttccgatcTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCt                                                   <  1:3264513/95‑1 (MQ=255)

                               ttccgatcTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCt                                                   <  1:1286780/95‑1 (MQ=255)

                               gAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCt                                                   <  1:2057955/100‑1 (MQ=255)

                               gAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCt                                                   <  1:4725377/100‑1 (MQ=255)

                                ttccgatctcacCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCtc                                                  <  1:5306527/91‑1 (MQ=255)

                                tccgatcTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCtc                                                  <  1:4224017/96‑1 (MQ=255)

                                tccgatcTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCtc                                                  <  1:3108166/96‑1 (MQ=255)

                                tccgatcTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCtc                                                  <  1:360681/96‑1 (MQ=255)

                                tccgatcTACACCCAAAGTTTTTAAAATAAGGACATGCTATTTGGTAACGAAAGCAGCACACACAGTAAACGGACCAAAGTAATCGCCACTTCCAGCCtt                                                  <  1:613537/96‑2 (MQ=255)

                                 ccgatcTGCACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCtca                                                 <  1:4069190/97‑1 (MQ=255)

                                 ccgatcTACATCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCtcg                                                 <  1:5042426/97‑2 (MQ=255)

                                 accatcTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGGCCAAAGTAATCGCCACTTCCAGCCtca                                                 <  1:2267439/98‑1 (MQ=255)

                                  cgatcTACGCCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCtcaa                                                <  1:3131826/98‑2 (MQ=255)

                                    atcTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCtcatca                                              <  1:1996255/100‑1 (MQ=255)

                                      cTACGCCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACt                                            <  1:5752458/100‑1 (MQ=255)

                                      cTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACt                                            <  1:5573775/100‑1 (MQ=255)

                                      cTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACt                                            <  1:4059047/100‑1 (MQ=255)

                                       tACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACtg                                           <  1:1995743/100‑2 (MQ=255)

                                       tACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTa                                           <  1:3962602/100‑1 (MQ=255)

                                         atcccAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTAcc                                         <  1:4708677/98‑1 (MQ=255)

                                          tcctaaaGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCa                                        <  1:5248479/96‑1 (MQ=255)

                                          tcccAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCa                                        <  1:4536428/99‑1 (MQ=255)

                                          gcccAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCAACTTCCAGCCTCATCACTACCa                                        <  1:4802617/99‑1 (MQ=255)

                                          actcAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCa                                        <  1:3261100/97‑1 (MQ=255)

                                          acCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCa                                        <  1:4147125/100‑1 (MQ=255)

                                          acCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCa                                        <  1:1888688/100‑1 (MQ=255)

                                           cccAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGGCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCag                                       <  1:2554762/100‑2 (MQ=255)

                                              aaaGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTGCCAATAc                                    <  1:1558294/100‑1 (MQ=255)

                                              aaaGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATAc                                    <  1:4721170/100‑1 (MQ=255)

                                              aaaGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATAc                                    <  1:5376623/100‑1 (MQ=255)

                                              aaaGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATAc                                    <  1:4441235/100‑1 (MQ=255)

                                              aaaGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATAc                                    <  1:109247/100‑1 (MQ=255)

                                              aaaGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATAc                                    <  1:1150374/100‑1 (MQ=255)

                                              aaaGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATAc                                    <  1:2818433/100‑1 (MQ=255)

                                               aaGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATAcg                                   <  1:4854196/100‑2 (MQ=255)

                                                aGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAg                                  <  1:826553/100‑1 (MQ=255)

                                                 gTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAgg                                 <  1:2305764/100‑2 (MQ=255)

                                                 gTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGt                                 <  1:1655104/100‑1 (MQ=255)

                                                  tttttAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGGCCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGtt                                <  1:2340410/100‑1 (MQ=255)

                                                  tttttAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGtt                                <  1:2693841/100‑1 (MQ=255)

                                                  tttttAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGtt                                <  1:2670367/100‑1 (MQ=255)

                                                   ttttAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGttc                               <  1:1576032/100‑2 (MQ=255)

                                                      tAAAATAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGTTaagg                            <  1:4579154/100‑3 (MQ=255)

                                                           tAGGGACATGCTCTTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGTTAAATTGATc                       <  1:3377892/100‑1 (MQ=255)

                                                                     ctctTTGGTAACGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGTTAAATTGATCATACacagag             <  1:5123935/100‑6 (MQ=255)

                                                                              aaCGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGTTAAATTGATCATACATTAAAGTTTGTTCa    <  1:3460589/100‑1 (MQ=255)

                                                                                cGAAAGCGGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGTTAAATTGATCATACATTAAAGTTTGTTCaag  <  1:4679672/100‑1 (MQ=255)

TTCTACAATTTTGGTATCGGTTAATTTTTTTGAATCATCTACACCCAAAGTTTTTAAAATAGGGACATGCTCTTTGGTAACGAAAGCAGCACACACAGTCAACGGACCAAAGTAATCGCCACTTCCAGCCTCATCACTACCAATACAGTTAAATTGATCATACATTAAAGTTTGTTCAAG  >  NC_007793/1136142‑1136321

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