BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,500,309	Δ1 bp	coding (333/363 nt)	b3372 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,500,308	0	C	.	99.2%	586.8 / ‑4.6	130	coding (332/363 nt)	b3372	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base . (65/64); minor base T (0/1); total (65/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

GCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCGCCGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTG  >  NC_000913/3500233‑3500379

gcgAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgc                                                                         <  1:3256300/76‑1 (MQ=255)

 ctaaaGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcg                                                                        <  1:2664302/74‑2 (MQ=255)

 cgAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcg                                                                        >  1:5615215/1‑75 (MQ=255)

 cgAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcg                                                                        <  2:5518655/76‑2 (MQ=255)

   aaaGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcggc                                                                      <  1:655570/76‑4 (MQ=255)

   aaaGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcggc                                                                      >  1:961430/1‑73 (MQ=255)

   aaaGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcggc                                                                      >  1:198365/1‑73 (MQ=255)

   aaaGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcggc                                                                      >  1:1021838/1‑73 (MQ=255)

     aGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGccc                                                                   >  2:5248346/1‑76 (MQ=255)

       ccTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAc                                                                 <  1:1966796/76‑1 (MQ=255)

       ccTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAc                                                                 <  2:1043080/76‑1 (MQ=255)

       ccTCGACATGAGCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAc                                                                 <  2:1884809/76‑1 (MQ=255)

          cGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAcgcg                                                              <  2:253688/76‑1 (MQ=255)

           gACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCgg                                                             <  1:3294317/76‑1 (MQ=255)

           gACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCgg                                                             <  2:2896636/76‑1 (MQ=255)

             cATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGc                                                           >  1:4166648/1‑76 (MQ=255)

              aTGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCt                                                          >  1:3154848/1‑76 (MQ=255)

              aTGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCt                                                          >  1:4733726/1‑76 (MQ=255)

               tGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTa                                                         <  2:592809/76‑1 (MQ=255)

               tGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTa                                                         >  2:319510/1‑76 (MQ=255)

               tGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTa                                                         <  2:1178005/76‑1 (MQ=255)

               tGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTa                                                         <  1:1470504/76‑1 (MQ=255)

               tGATCAAGCTGGCGAGGGATATGGCAAAAGAGATGAACGCGGGATATACGCTGATTTCCG‑CGGCCCACGCGGGCTa                                                         >  1:3983065/1‑76 (MQ=255)

                gATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAt                                                        >  1:1681178/1‑76 (MQ=255)

                 aTCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAtc                                                       <  2:2519290/76‑1 (MQ=255)

                 aTCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAtc                                                       >  1:766890/1‑76 (MQ=255)

                 aTCAAGCTGGCGATGGATATGGAAAAAGAGGTGAACGAGGGTTATAGGCTGATGTCCG‑CGGCCCACGCGGGATagc                                                       >  1:741007/1‑74 (MQ=255)

                  tCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAtct                                                      <  2:4708104/76‑1 (MQ=255)

                    aaGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCa                                                    >  1:5361365/1‑76 (MQ=255)

                     aGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCAc                                                   <  1:3896265/76‑1 (MQ=255)

                     aGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCAc                                                   <  2:3217412/76‑1 (MQ=255)

                     aGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCAc                                                   <  1:2497589/76‑1 (MQ=255)

                     aGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCAc                                                   <  2:1553548/76‑1 (MQ=255)

                     aGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCAc                                                   <  1:2028895/76‑1 (MQ=255)

                      gCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACg                                                  <  2:3386926/76‑1 (MQ=255)

                      gCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGATATCTCACg                                                  >  2:165010/1‑76 (MQ=255)

                        tGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGcc                                                >  2:4990230/1‑76 (MQ=255)

                        tGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGcc                                                >  1:1923248/1‑76 (MQ=255)

                        tGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGcc                                                <  1:3266471/76‑1 (MQ=255)

                        tGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGcc                                                <  1:4023603/76‑1 (MQ=255)

                        tGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCCATCTCAcccc                                                >  2:10944/1‑76 (MQ=255)

                          gCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCAc                                              <  2:3968465/76‑1 (MQ=255)

                            gATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCt                                            >  1:1030350/1‑76 (MQ=255)

                            gATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCt                                            >  1:414718/1‑76 (MQ=255)

                             aTGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTa                                           >  2:2324266/1‑76 (MQ=255)

                               ggATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAAt                                         >  1:2152453/1‑76 (MQ=255)

                                 atatGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGt                                       >  1:1812312/1‑76 (MQ=255)

                                 atatGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGt                                       >  2:3940123/1‑76 (MQ=255)

                                  tatGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGtt                                      >  2:654779/1‑76 (MQ=255)

                                  tatGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGtt                                      <  1:4143444/76‑1 (MQ=255)

                                  gatGGCAAAAGAGAGGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGtt                                      <  2:91364/75‑1 (MQ=255)

                                   atGGCAAAAGAGATGAACGCGGGCTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTa                                     >  1:1518106/1‑76 (MQ=255)

                                    tGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTAt                                    <  1:3909536/76‑1 (MQ=255)

                                     ggCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCGCCTAATGTTATc                                   <  1:4768384/76‑1 (MQ=255)

                                      gCAAAAGGGATGAGCGCGGGGTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCt                                  <  2:4629435/76‑1 (MQ=255)

                                      gCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTAtat                                  >  1:2991431/1‑74 (MQ=255)

                                         aaaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTggg                               >  1:2167305/1‑76 (MQ=255)

                                          aaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTgggg                              <  2:3449688/76‑1 (MQ=255)

                                          aaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTgggg                              >  1:3847970/1‑76 (MQ=255)

                                           agagCTGAACTCGGGTTATACGATGATCTCCG‑CGTCCCACGCGGGCTATCCCACGCCCCCTAAGGTTACCTGGGGa                             <  2:5139018/76‑1 (MQ=255)

                                           agagATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGa                             <  2:4865014/76‑1 (MQ=255)

                                           agagATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGa                             <  2:414718/76‑1 (MQ=255)

                                            gagaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCGCCTAATGTTATCTGGGGAc                            <  1:789131/76‑1 (MQ=255)

                                            gagaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGAc                            <  2:2918107/76‑1 (MQ=255)

                                             agaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGAc                            <  1:5576521/75‑1 (MQ=255)

                                             agaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGAc                            >  2:5576519/1‑75 (MQ=255)

                                             agaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACg                           <  2:4036447/76‑1 (MQ=255)

                                             agaTGAACGCGGGTTATACGCTGATTTCCG‑CGGACCACGCGGGCTATCTCACGCCACCGAATGTTATCTGGGGACg                           >  2:2014059/1‑76 (MQ=255)

                                              gctGAACGCGGGTTATACGCTGATTTCCG‑CGGGCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACgg                          <  2:1812012/74‑1 (MQ=255)

                                               aTGAGCGCGGGTTATGCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGc                         <  2:3282967/76‑1 (MQ=255)

                                               aTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGc                         >  2:1975944/1‑76 (MQ=255)

                                                tGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCt                        >  2:2128596/1‑76 (MQ=255)

                                                ggAACGCGGGTTATACGGTGATTTCCG‑CGGCCCGCGCGGGGTATCTCACGCCACCTAATGTTATCTGGGGACGGCt                        <  2:5538417/75‑1 (MQ=255)

                                                   aCGCGGGTTATACGCTGCTTTCCG‑CGGCCCACGCGGTCTATCTCACGCCACCTAATGTTCTCTTGGGACGGCTGGc                     >  1:4465113/1‑76 (MQ=255)

                                                   aCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGc                     >  1:3394764/1‑76 (MQ=255)

                                                   aCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGc                     >  2:951911/1‑76 (MQ=255)

                                                    cgcgGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCCCCTAATGCTATCTGGGGCCCGCCGGcc                    >  2:557167/1‑75 (MQ=255)

                                                    cgcgGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCa                    >  1:2894506/1‑76 (MQ=255)

                                                     gcgGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCCCGCCCCCTAATGTTATCTGGGGACGGCTGGccg                   >  2:5779050/1‑74 (MQ=255)

                                                     gcgGGTTATACGCTGATTTCC‑TCGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGTCTGGCag                   <  1:2014059/76‑1 (MQ=255)

                                                      cgGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCaga                  <  2:1914050/76‑1 (MQ=255)

                                                       gggTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCagag                 <  1:5445194/76‑1 (MQ=255)

                                                        ggTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAgaga                <  1:5199082/76‑1 (MQ=255)

                                                        ggTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAgaga                >  1:845003/1‑76 (MQ=255)

                                                          gtataCGCCGATTTCCG‑GGGCCCCCGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAAc              <  2:5668978/75‑1 (MQ=255)

                                                             taCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTa           >  1:510939/1‑76 (MQ=255)

                                                             taCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTa           <  1:4387046/76‑1 (MQ=255)

                                                              aCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCa                    >  1:1257346/1‑66 (MQ=255)

                                                              aCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCa                    <  2:1257346/66‑1 (MQ=255)

                                                              aCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTaa          >  1:4261874/1‑76 (MQ=255)

                                                              aCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTaa          >  1:5419962/1‑76 (MQ=255)

                                                              aCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTaa          <  2:4249686/76‑1 (MQ=255)

                                                                gCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGc        <  2:1817033/76‑1 (MQ=255)

                                                                gCTGATTTCCG‑AGGCACACGCGGGCTATCCCACGCCACCTAATGTTGTCTGGGGTCGGCTGGCAGCGAAACCAAGc        >  2:169798/1‑76 (MQ=255)

                                                                 cTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCg       <  2:4011285/76‑1 (MQ=255)

                                                                  tGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTCATGTTATCTGGGGACGGCTCGCAGAGCACCTCAGcca      >  2:5022906/1‑74 (MQ=255)

                                                                  tGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGa      >  2:811277/1‑76 (MQ=255)

                                                                     tttCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGCGACCCTAAGCGAACt   >  2:3904461/1‑76 (MQ=255)

                                                                      ttCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACtg  >  1:1765781/1‑76 (MQ=255)

                                                                      ttCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACtg  >  1:2766578/1‑76 (MQ=255)

GCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCGCCGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTG  >  NC_000913/3500233‑3500379

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 9 ≤` `ATCG/ATCG` `< 19 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`