BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,118,190	C→T	G27D (GGC→GAC)	b1057 ←	putative cytochrome

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,118,190	0	C	T	97.2%	237.9 / ‑3.4	71	G27D (GGC→GAC)	b1057	putative cytochrome
Reads supporting (aligned to +/- strand): ref base C (0/2); new base T (19/50); total (19/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

ATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACATAAGTGCA  >  NC_000913/1118100‑1118276

aTACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCg                                                                                      >  1:3748808/1‑93 (MQ=255)

 tACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGAACTAAAGCAAACATGTCGt                                                                                     >  2:4385338/1‑93 (MQ=255)

  aCCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTa                                                                                    >  1:319578/1‑93 (MQ=255)

     gATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGAc                                                                                 >  2:2357817/1‑93 (MQ=255)

       tGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACaa                                                                               >  1:1233321/1‑93 (MQ=255)

       tGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACaa                                                                               >  2:1691572/1‑93 (MQ=255)

        gCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACaat                                                                              <  2:4130011/93‑1 (MQ=255)

        gCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACaat                                                                              <  2:2039736/93‑1 (MQ=255)

            ttATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAAtt                                                                          <  1:1776305/93‑1 (MQ=255)

             tATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTg                                                                         <  1:4104823/93‑1 (MQ=255)

              aTGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTgc                                                                        <  1:3509981/93‑1 (MQ=255)

              aTGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTgc                                                                        <  1:3683729/93‑1 (MQ=255)

                  aGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCtt                                                                    >  2:1079371/1‑93 (MQ=255)

                  aGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCtt                                                                    >  2:162096/1‑93 (MQ=255)

                   gTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTa                                                                   >  2:532584/1‑93 (MQ=255)

                   gTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTa                                                                   >  1:2010150/1‑93 (MQ=255)

                   gTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTa                                                                   >  2:4129423/1‑93 (MQ=255)

                    ttCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAg                                                                  <  1:2290092/93‑1 (MQ=255)

                      cGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGcc                                                                <  1:211058/93‑1 (MQ=255)

                        ggTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCaa                                                              >  2:176997/1‑93 (MQ=255)

                        ggTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCaa                                                              <  2:3145698/93‑1 (MQ=255)

                        ggTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCaa                                                              <  2:4011445/93‑1 (MQ=255)

                         gTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAAt                                                             >  1:549076/1‑93 (MQ=255)

                         gTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAAt                                                             <  2:463476/93‑1 (MQ=255)

                          tGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATg                                                            <  2:654946/93‑1 (MQ=255)

                          tGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATg                                                            <  1:2110232/93‑1 (MQ=255)

                           gCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATgg                                                           <  1:1025631/93‑1 (MQ=255)

                           gCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATgg                                                           <  1:733533/93‑1 (MQ=255)

                             tttGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGaa                                                         <  2:2726588/93‑1 (MQ=255)

                              ttGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAg                                                        >  1:3348974/1‑93 (MQ=255)

                               tgtgATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAgg                                                       <  1:1010619/93‑1 (MQ=255)

                                ttgATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGc                                                      <  1:918360/92‑1 (MQ=255)

                                gtgATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGc                                                      <  1:1876124/93‑1 (MQ=255)

                                gtgATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGc                                                      <  2:4388253/93‑1 (MQ=255)

                                gtgATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGc                                                      <  1:917226/93‑1 (MQ=255)

                                 tgATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGcc                                                     <  2:4662526/93‑1 (MQ=255)

                                 tgATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGcc                                                     <  1:1249877/93‑1 (MQ=255)

                                      ccAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGa                                                <  1:3113922/93‑1 (MQ=255)

                                      ccAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGa                                                <  1:3667028/93‑1 (MQ=255)

                                         gCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGatat                                             <  1:2632406/93‑1 (MQ=255)

                                          ccATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGAtata                                            <  2:2160286/93‑1 (MQ=255)

                                           cATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATaa                                           <  1:1596684/93‑1 (MQ=255)

                                            aTCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAac                                          <  2:3834503/93‑1 (MQ=255)

                                            aTCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAac                                          >  1:1854268/1‑93 (MQ=255)

                                            aTCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAac                                          <  1:3275235/93‑1 (MQ=255)

                                                taataaCTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCa                                      <  1:3969149/93‑1 (MQ=255)

                                                taataaCTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCa                                      >  2:3069536/1‑93 (MQ=255)

                                                 aataaCTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCAt                                     <  1:301364/93‑1 (MQ=255)

                                                      cTGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGt                                <  1:2366692/93‑1 (MQ=255)

                                                       tGAGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGtt                               >  1:2165539/1‑93 (MQ=255)

                                                         aGCGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCa                             <  1:1563669/93‑1 (MQ=255)

                                                           cttGACCATCCACAGACATAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAgg                           <  2:1358095/91‑1 (MQ=255)

                                                           cGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAgg                           <  2:4270331/93‑1 (MQ=255)

                                                           cGTGACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAgg                           <  1:566317/93‑1 (MQ=255)

                                                              gACCATCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGt                        >  1:1407353/1‑93 (MQ=255)

                                                                   tCCACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTtg                   >  2:3824325/1‑93 (MQ=255)

                                                                    ccACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTtgt                  <  2:1011928/93‑1 (MQ=255)

                                                                    ccACAGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTtgt                  <  2:2470362/93‑1 (MQ=255)

                                                                       caGACCTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGaa               <  1:3916880/93‑1 (MQ=255)

                                                                        aGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAAt              <  2:2546021/93‑1 (MQ=255)

                                                                            cTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACa          <  2:2567670/93‑1 (MQ=255)

                                                                            cTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACa          >  1:2071741/1‑93 (MQ=255)

                                                                            cTAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACa          <  1:1009290/93‑1 (MQ=255)

                                                                             tAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACAt         <  1:2151956/93‑1 (MQ=255)

                                                                             tAAAGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACAt         <  1:1746756/93‑1 (MQ=255)

                                                                              aaaGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACATa        <  1:4601798/93‑1 (MQ=255)

                                                                              aaaGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACATa        <  1:1792690/93‑1 (MQ=255)

                                                                               aaGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACATaa       <  2:134881/93‑1 (MQ=255)

                                                                               aaGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACATaa       <  2:1368116/93‑1 (MQ=255)

                                                                                aGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACATAAg      <  1:1294436/93‑1 (MQ=255)

                                                                                aGCAAACATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACATAAg      <  2:3798432/93‑1 (MQ=255)

                                                                                    aaCATGTCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACATAAGTGCa  <  1:595331/93‑1 (MQ=255)

ATACCGATGCTTTTATGCAGTTCGGGTGCTTTGTGATACCAGCCATCGTAATAACTGAGCGTGACCATCCACAGACCTAAAGCAAACATGCCGTAGACAATAATTGCGCTTAGCCAATGGAAGGCCGCAGATATAACACCATAGCGTTCAGGGGTATTTGTGAATGACATAAGTGCA  >  NC_000913/1118100‑1118276

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