BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,311,496	+C	coding (301/408 nt)	phnQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,311,495	1	.	C	99.3%	419.5 / ‑3.8	149	coding (302/408 nt)	phnQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (91/57); total (91/58)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 164 total aligned reads displayed.

GATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTAT  >  NC_000913/4311420‑4311565

 aTAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑c                                                                       >  1:4212683/1‑76 (MQ=255)

 aTAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑c                                                                       >  2:5497467/1‑76 (MQ=255)

 aTAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑c                                                                       >  1:4948842/1‑76 (MQ=255)

 aTAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATGCACCCTGCTCACGCTGTCCA‑c                                                                       <  2:1090453/76‑1 (MQ=255)

  tAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑cc                                                                      >  1:2343653/1‑75 (MQ=255)

  tAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑cc                                                                      <  1:1461569/76‑2 (MQ=255)

   agagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑ccg                                                                     <  1:2484652/76‑3 (MQ=255)

   agagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑ccg                                                                     <  1:1778502/76‑3 (MQ=255)

    gagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑ccgt                                                                    >  2:2923158/1‑73 (MQ=255)

       tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTcc                                                                  <  1:1738447/76‑1 (MQ=255)

       tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTcc                                                                  <  1:60477/76‑1 (MQ=255)

        ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCg                                                                 <  1:168383/76‑1 (MQ=255)

         tCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCgg                                                                <  1:1659787/76‑1 (MQ=255)

         tCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCgg                                                                >  1:4333315/1‑76 (MQ=255)

         tCTTGTTGAGATTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCgg                                                                <  2:4936421/76‑1 (MQ=255)

          cTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGc                                                               <  2:3891772/76‑1 (MQ=255)

           tttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCCCCCTGCTCACGCTGTCCACCGTTCCGGcc                                                              <  1:3441817/76‑1 (MQ=255)

           tttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGcc                                                              >  1:2060186/1‑76 (MQ=255)

           tttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGcc                                                              >  1:308769/1‑76 (MQ=255)

             tttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCa                                                            <  1:3295615/76‑1 (MQ=255)

              ttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAg                                                           <  2:1630341/76‑1 (MQ=255)

               tGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAgg                                                          >  2:3818526/1‑76 (MQ=255)

                gAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGc                                                         >  1:3907886/1‑76 (MQ=255)

                 aGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCGCCGTTCCGGCCCAGGcc                                                        <  2:3248430/76‑1 (MQ=255)

                 aGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGcc                                                        >  2:279743/1‑76 (MQ=255)

                 aGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGcc                                                        <  1:2681785/76‑1 (MQ=255)

                  ggTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCg                                                       >  1:1970509/1‑76 (MQ=255)

                  gcttGGCTAATGACGCGCTTTCCTCTATTTCTTTTATTCTCCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCg                                                       <  1:468198/74‑1 (MQ=255)

                   gTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGa                                                      >  1:1080690/1‑76 (MQ=255)

                   gTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGa                                                      >  2:5748925/1‑76 (MQ=255)

                    ttGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGat                                                     >  1:2892991/1‑76 (MQ=255)

                     tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGata                                                    >  2:111782/1‑76 (MQ=255)

                      ggCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGatat                                                   <  1:4320346/76‑1 (MQ=255)

                        gtAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtatta                                                 <  2:5288779/75‑1 (MQ=255)

                        ctcaTGACACGCTTTCCGCTATTTCTTTTCTTCACCCTGCTCACGCTCCCCACCGTTCCGGCCCAGGCCGAtatta                                                 <  1:451217/73‑1 (MQ=255)

                        cTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtatta                                                 >  2:860786/1‑76 (MQ=255)

                        cTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtatta                                                 >  1:2547051/1‑76 (MQ=255)

                         tAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtattat                                                <  2:830244/76‑1 (MQ=255)

                         tAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtattat                                                >  1:97284/1‑76 (MQ=255)

                          aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc                                               >  2:2659068/1‑76 (MQ=255)

                          aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc                                               <  2:443925/76‑1 (MQ=255)

                          aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc                                               >  2:4531290/1‑76 (MQ=255)

                           aTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcg                                              <  2:1316321/76‑1 (MQ=255)

                           aTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcg                                              >  2:1203626/1‑76 (MQ=255)

                           aTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcg                                              <  1:3182608/76‑1 (MQ=255)

                            tGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcga                                             >  2:686301/1‑76 (MQ=255)

                            tGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcga                                             <  1:5718883/76‑1 (MQ=255)

                              aaaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgacg                                           >  1:918400/1‑76 (MQ=255)

                               aaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgacga                                          <  2:366741/76‑1 (MQ=255)

                                  gCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGacac                                       >  1:1311110/1‑76 (MQ=255)

                                  gCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGacac                                       >  2:1679428/1‑76 (MQ=255)

                                  gCTTTCCGCTAGTTCTTTTATTCACCCTGCTCACGCTCGCCACCGTTCCGCCCCAGGCCGATATTATCGACCacac                                       >  2:5330763/1‑76 (MQ=255)

                                   cTTTCCGCTATTTCTTTTATACACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACAcc                                      <  2:3486459/76‑1 (MQ=255)

                                    tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCa                                     >  2:2028538/1‑76 (MQ=255)

                                    tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCa                                     <  2:2171373/76‑1 (MQ=255)

                                    tttCCGCTATTTATTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCa                                     >  1:1941353/1‑76 (MQ=255)

                                     ttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCgacgac                                         >  1:4081847/1‑71 (MQ=255)

                                     ttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCgacgac                                         <  2:4081846/71‑1 (MQ=255)

                                     ttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAt                                    >  1:3692955/1‑76 (MQ=255)

                                     gcccGCTATTTCTTTCATTCACCCTGCTCACGCTGCCCACCGCTCCGGCCCAGGCCGATATTATCGACCACACCac                                    >  2:5137856/3‑75 (MQ=255)

                                       ccGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGccaccagtc                                  >  2:2367128/1‑73 (MQ=255)

                                        cGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                 >  1:1361907/1‑76 (MQ=255)

                                        cGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                 >  1:1176749/1‑76 (MQ=255)

                                        cGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                 >  1:2603475/1‑76 (MQ=255)

                                         gCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGc                                >  2:2138761/1‑76 (MQ=255)

                                         gCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGc                                >  2:4996010/1‑76 (MQ=255)

                                         gCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGc                                >  2:4321030/1‑76 (MQ=255)

                                          cTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCa                               <  1:2938490/76‑1 (MQ=255)

                                           tATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                              >  1:4332963/1‑76 (MQ=255)

                                            aTTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                              >  2:1934593/1‑75 (MQ=255)

                                            aTTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTCCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                              <  1:1934593/75‑1 (MQ=255)

                                              ttCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAatat                           >  1:2177257/1‑76 (MQ=255)

                                               tCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATAtt                          <  2:2085344/76‑1 (MQ=255)

                                               tCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATAtt                          >  2:2823161/1‑76 (MQ=255)

                                                  tttATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACAATTGGCAATATTcac                       >  2:5695142/1‑75 (MQ=255)

                                                   ttATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                      <  2:3897268/76‑1 (MQ=255)

                                                   ttATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                      <  1:3845456/76‑1 (MQ=255)

                                                     aTTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                    <  2:2383729/76‑1 (MQ=255)

                                                     aTTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                    >  1:1564817/1‑76 (MQ=255)

                                                      ttCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAgg                   >  2:3323589/1‑76 (MQ=255)

                                                      ttCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCC‑ATATTATCCACGACCACAATTGCAATAATcaacaagg                  >  1:5136595/1‑73 (MQ=255)

                                                      ttCACCCTGCTCACGCTGTCCACCGTTCCGGCCCACGCCGATATTATCGACGACACCATTGGCAACATTCAGCAgg                   >  2:788272/1‑76 (MQ=255)

                                                       tCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGc                  >  2:3256125/1‑76 (MQ=255)

                                                        cACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCg                 >  2:1848195/1‑76 (MQ=255)

                                                        cACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCg                 <  1:1596536/76‑1 (MQ=255)

                                                         ccccTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAACAGGCGa                <  2:5024751/75‑1 (MQ=255)

                                                         aCCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGa                >  2:3143772/1‑76 (MQ=255)

                                                         aCCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGa                >  1:3734479/1‑76 (MQ=255)

                                                         aCCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGa                >  2:4911574/1‑76 (MQ=255)

                                                         aCCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGa                >  2:1788099/1‑76 (MQ=255)

                                                          cccTGCTCACGCTGTCCCCCGTCCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt               <  1:2923158/76‑1 (MQ=255)

                                                           ccTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAtt              <  1:5064932/76‑1 (MQ=255)

                                                              gCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAAt           >  1:5012231/1‑76 (MQ=255)

                                                              gCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAAt           >  2:3606409/1‑76 (MQ=255)

                                                                tCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGa         <  2:2343653/76‑1 (MQ=255)

                                                                 cACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCCGCAGGCGATTAATGAc        >  2:3706339/1‑76 (MQ=255)

                                                                    gCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt     >  1:1965931/1‑76 (MQ=255)

                                                                    gCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt     >  1:978145/1‑76 (MQ=255)

                                                                     cTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCtt    >  2:2764039/1‑76 (MQ=255)

                                                                       gTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtat  <  2:308769/76‑1 (MQ=255)

GATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTAT  >  NC_000913/4311420‑4311565

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