BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,747,389	Δ1 bp	coding (572/591 nt)	b4540 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,747,388	0	C	.	98.0%	442.5 / 1.6	102	coding (571/591 nt)	b4540	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/1); new base . (51/48); total (53/49)
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.53e-01

Only 100 of 110 total aligned reads displayed.

TTGCGCACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATCCTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGA  >  NC_000913/2747313‑2747459

ttGCGCACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATc                                                                         >  1:4660118/1‑76 (MQ=255)

ttGCGCACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATc                                                                         <  1:1944195/76‑1 (MQ=255)

ttGCGCACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATc                                                                         <  1:1508894/76‑1 (MQ=255)

 tGCGCACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATct                                                                        >  1:3568425/1‑75 (MQ=255)

 tGCGCACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATct                                                                        >  2:1108778/1‑75 (MQ=255)

  gcgcACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATctg                                                                       >  2:81095/1‑74 (MQ=255)

  gcgcACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATctg                                                                       <  2:1192461/76‑3 (MQ=255)

     cACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGaa                                                                   >  1:1977127/1‑76 (MQ=255)

       cTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGaaaa                                                                 <  2:4282148/76‑1 (MQ=255)

        tATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGCCGGTGCTCGAT‑CTGGaaaaa                                                                >  1:3617363/1‑76 (MQ=255)

         ctCGTGCGCGAATCACGCTCACAAATTTCCCTTCGCAATCAGTATATGCTGCTGTCGGTGCTCGAT‑CTGGaaaaac                                                               <  1:1533985/75‑2 (MQ=255)

         aTCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGaaaaaa                                                               >  2:1895505/1‑76 (MQ=255)

          tCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAAGCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAAt                                                              <  2:2492313/76‑1 (MQ=255)

           cGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATg                                                             >  2:3128130/1‑76 (MQ=255)

           cGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATg                                                             <  1:2213269/76‑1 (MQ=255)

           cGTGCACGAATCACGCTAACAAAGTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATg                                                             <  1:2222359/76‑1 (MQ=255)

            gTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATgc                                                            >  2:4636778/1‑75 (MQ=255)

            gTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGa                                                            <  2:1690279/76‑1 (MQ=255)

            gTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGa                                                            >  1:848997/1‑76 (MQ=255)

            gTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGa                                                            <  2:1492956/76‑1 (MQ=255)

            gTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGa                                                            <  1:3422571/76‑1 (MQ=255)

              gCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGAcc                                                          >  2:122915/1‑76 (MQ=255)

              gCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGAcc                                                          <  1:3789736/76‑1 (MQ=255)

               cACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGCT‑CTGGAAAAAATGcccc                                                         >  2:355613/1‑75 (MQ=255)

               cACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCg                                                         <  1:4869529/76‑1 (MQ=255)

                 cGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGtt                                                       <  2:3107419/76‑1 (MQ=255)

                  gAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTtg                                                      >  1:4587518/1‑76 (MQ=255)

                   aaTCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTtga                                                     >  1:897374/1‑76 (MQ=255)

                    aTCACGCTCACAAATTTCCCGTCTCAATCAGGATATGGTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTtgat                                                    <  2:3526072/76‑1 (MQ=255)

                     tCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTtgatg                                                   >  2:1802403/1‑76 (MQ=255)

                      cACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTtgatga                                                  <  2:1784748/76‑1 (MQ=255)

                        cGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAca                                                <  2:2344404/76‑1 (MQ=255)

                        cGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGGGCTCGAT‑CTGGAAAAAATGACCGTTGGTGTca                                                >  1:1885645/1‑76 (MQ=255)

                          cTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAcatc                                              <  1:4039118/76‑1 (MQ=255)

                          cTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAcatc                                              >  1:1929800/1‑76 (MQ=255)

                          cTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAcata                                              <  2:2048753/76‑2 (MQ=255)

                           tCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAcatca                                             <  2:1208879/76‑1 (MQ=255)

                           tCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCCAT‑CTGGAAAAAATGACCGTTGATGAcatca                                             >  1:1640446/1‑76 (MQ=255)

                            cacaAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAcatcat                                            <  2:3278934/76‑1 (MQ=255)

                            cacaAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGAcatcat                                            <  2:3927038/76‑1 (MQ=255)

                            cacaAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATCAcatcat                                            >  1:3537209/1‑76 (MQ=255)

                             acaAATTTCCCGTCGCAATCAGGATATGGTGCTGTGGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATg                                           <  2:1888122/76‑1 (MQ=255)

                             acaAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATg                                           >  2:3033421/1‑76 (MQ=255)

                               aaaTTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGCCCGTTGATGACATCATGGt                                         <  1:4154069/76‑1 (MQ=255)

                                aaTTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGA‑GCTGGAAAAAATGACCGATAATGACATCATGGTg                                        >  2:4171503/1‑76 (MQ=255)

                                 agttCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGAGc                                       >  1:125964/3‑76 (MQ=255)

                                   ttCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCcg                                     <  1:2312234/76‑1 (MQ=255)

                                   ttCCCGTAGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGGCATCATGGTGccn                                     <  1:3531210/76‑2 (MQ=255)

                                      ccGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCgcgc                                  <  1:1607676/76‑1 (MQ=255)

                                      ccGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCgcgc                                  <  1:1272783/76‑1 (MQ=255)

                                      ccGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCgcgc                                  >  2:154299/1‑76 (MQ=255)

                                      ccGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCgcgc                                  >  2:3504673/1‑76 (MQ=255)

                                       cGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCa                                 <  1:2211804/76‑1 (MQ=255)

                                         tCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGt                               >  2:2592144/1‑76 (MQ=255)

                                         tCGCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGt                               <  1:2281706/76‑1 (MQ=255)

                                           gCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGa                             >  1:2625219/1‑76 (MQ=255)

                                           gCAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGa                             >  2:2022569/1‑76 (MQ=255)

                                            cAATCAGGATGTGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGaa                            <  2:4295879/76‑1 (MQ=255)

                                            cAATCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGaa                            >  2:589508/1‑76 (MQ=255)

                                            cAATCAGGATATGCTGCGGTCGGTGCTCGAT‑CTGGAAAAAATGGCCTTTGATTACATCATGGTGCCGCGCGGTgca                            <  2:1640111/76‑3 (MQ=255)

                                            cAATCAGGAGATGATGCTGCCGGTACATGAT‑CGGGTAAAAATGACCGTTGCAGACATCATTGCGCCCGGCAGTgta                            >  1:778632/1‑74 (MQ=255)

                                             aaTCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAAGGACCGTTGATGACATCATGGTGCCGCGCAGTGaaa                           <  2:3121150/76‑1 (MQ=255)

                                              aTCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAAt                          >  2:2802106/1‑76 (MQ=255)

                                              aTCAGGATATGCTGCTGTCGGGGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGTAGTGAAAt                          >  2:4626605/1‑76 (MQ=255)

                                               tCAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAAtt                         >  1:3176889/1‑76 (MQ=255)

                                                cAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTa                        >  1:2999880/1‑76 (MQ=255)

                                                cAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCCCCGGGAAAtca                        >  2:4798364/1‑74 (MQ=255)

                                                cAGGATATGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGCTGATGACATCATGGGGCCGCCCAGCGAAAttc                        >  2:4823188/1‑75 (MQ=255)

                                                    ggatGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCgg                    <  2:169843/74‑1 (MQ=255)

                                                    atatGCTGCTGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCga                    <  1:3128130/76‑2 (MQ=255)

                                                      aggctgctGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTa                  <  1:4466202/74‑1 (MQ=255)

                                                       tgctgcTGTCGGTGGTCGATCGGGGAAAAAA‑AACGGCTGACGGAATCACGGGGACGCGAAGTGAAATTATTGGTAt                 >  1:1731075/1‑76 (MQ=255)

                                                        gctgctGTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTAtt                >  2:3633507/1‑76 (MQ=255)

                                                         ctgctgTCGGTGATCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTg               >  1:1116010/1‑76 (MQ=255)

                                                          tgctgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGa              >  2:4107523/1‑76 (MQ=255)

                                                          tgctgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGa              >  1:372398/1‑76 (MQ=255)

                                                          tgctgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGa              <  1:1677523/76‑1 (MQ=255)

                                                          tgctgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGa              >  2:3275090/1‑76 (MQ=255)

                                                           gctgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGat             <  2:1776264/76‑1 (MQ=255)

                                                           gctgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGat             >  2:1083781/1‑76 (MQ=255)

                                                            gtgTCGGGGCTGGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGata            <  2:271948/75‑1 (MQ=255)

                                                            ctgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGata            >  2:3275462/1‑76 (MQ=255)

                                                             tgTCGTCGCTCGA‑CCTGGAAAAAATGACCGTTGATTACATCATGGTGCCGCGCAGTGAAATTATCGGTATTTatat           <  2:1063030/76‑1 (MQ=255)

                                                             tgTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGatat           <  1:2875547/76‑1 (MQ=255)

                                                             ggTCGGTGCTCGAT‑CTGGAAAAAATGACCGGTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGatat           <  2:4781957/75‑1 (MQ=255)

                                                              gTCGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGAATTGATATc          >  2:1755446/1‑76 (MQ=255)

                                                               tCGGTGCTCGAT‑CTGGAAAAAATGACCGGTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCa         >  2:294002/1‑76 (MQ=255)

                                                                cGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCaa        >  1:526437/1‑76 (MQ=255)

                                                                cGGTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCaa        >  1:675304/1‑76 (MQ=255)

                                                                 ggTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAAc       <  2:2396106/76‑1 (MQ=255)

                                                                 ggTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAAc       <  1:2134480/76‑1 (MQ=255)

                                                                 ggTGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAAc       >  1:2819901/1‑76 (MQ=255)

                                                                  gTGCTCGA‑GCTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACg      <  2:4803192/76‑1 (MQ=255)

                                                                   tGCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGa     >  1:1397258/1‑76 (MQ=255)

                                                                    gCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGAt    >  2:2041791/1‑76 (MQ=255)

                                                                    gCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGAt    <  2:3227778/76‑1 (MQ=255)

                                                                    gCTCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGAt    >  2:4794055/1‑76 (MQ=255)

                                                                      tCGAT‑CTGGAAAAAATGAGTGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGa  <  2:4070831/76‑1 (MQ=255)

                                                                      tCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGa  >  2:4454302/1‑76 (MQ=255)

                                                                      tCGAT‑CTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGa  <  2:1637889/76‑1 (MQ=255)

TTGCGCACTATCGTGCACGAATCACGCTCACAAATTTCCCGTCGCAATCAGGATATGCTGCTGTCGGTGCTCGATCCTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGA  >  NC_000913/2747313‑2747459

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