BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,359,867	G→A	40.2%	noncoding (1228/2905 nt)	rrlA ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,359,867	0	G	A	40.2%	71.2 / 94.7	92	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (24/31); new base A (21/16); total (45/47)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTG  >  Exported/2359818‑2359931

gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                             <  1:4314411/71‑1 (MQ=34)

gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                             >  1:3216557/1‑71 (MQ=255)

gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                             >  1:2000973/1‑71 (MQ=255)

gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                             >  1:4602559/1‑71 (MQ=255)

 cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                <  1:1836757/67‑1 (MQ=31)

 cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                <  1:5153493/67‑1 (MQ=31)

 cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                            >  1:331885/1‑71 (MQ=34)

  tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                           >  1:5623932/1‑71 (MQ=34)

  tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                           >  1:4465113/1‑71 (MQ=255)

   ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                          >  1:1573809/1‑71 (MQ=255)

    tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                          >  1:4457524/1‑70 (MQ=255)

      tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                        >  1:4141641/1‑70 (MQ=35)

        gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                      <  1:1137304/70‑1 (MQ=35)

        gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     <  1:1456590/71‑1 (MQ=34)

        gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     <  1:230144/71‑1 (MQ=34)

        gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                     <  1:1034871/71‑1 (MQ=255)

         ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     >  1:2651374/1‑70 (MQ=35)

         ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     >  1:2527219/1‑70 (MQ=35)

         ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     >  1:5405943/1‑70 (MQ=35)

         ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCtt                                         <  1:4459936/66‑1 (MQ=255)

            cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 <  1:468211/71‑1 (MQ=34)

              tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                               >  1:3215094/1‑71 (MQ=255)

                tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 <  1:1179536/67‑1 (MQ=35)

                tGTCAGCATTCGCACTTCTGATACCTCCAGCATACGTCACAGTACACCTTCACAGGCTTACAGAACGCTc                              <  1:3012951/70‑1 (MQ=11)

                tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                              <  1:3318125/70‑1 (MQ=255)

                tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                              <  1:1450232/70‑1 (MQ=255)

                 gTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                            >  1:474960/1‑71 (MQ=34)

                  tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                        <  1:3607846/58‑1 (MQ=35)

                  tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                           <  1:5434838/71‑1 (MQ=255)

                   cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 >  1:4955679/1‑64 (MQ=35)

                   cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 >  1:1059741/1‑64 (MQ=35)

                   cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 >  1:5394983/1‑64 (MQ=35)

                   cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 >  1:786363/1‑64 (MQ=35)

                    aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                          <  1:1678427/70‑1 (MQ=35)

                    aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                          <  1:3078047/70‑1 (MQ=35)

                    aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                          <  1:4822968/70‑1 (MQ=255)

                    aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTAACAGGCTTACAGAACGCTCCCCTa                         <  1:1551669/71‑1 (MQ=11)

                     gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                     <  1:3139555/58‑1 (MQ=11)

                      cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                          >  1:1774057/1‑68 (MQ=255)

                        ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                        >  1:2011287/1‑68 (MQ=35)

                          cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                        >  1:2394865/1‑50 (MQ=38)

                          cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                    >  1:3728057/1‑70 (MQ=35)

                          cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                   >  1:5267017/1‑55 (MQ=11)

                          cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                   >  1:1050604/1‑55 (MQ=11)

                          cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                   >  1:5664140/1‑55 (MQ=11)

                          cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                   >  1:2862491/1‑55 (MQ=11)

                          cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGACCGCTCCCCTACCcaa                    >  1:2506401/1‑70 (MQ=11)

                          cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                    >  1:3464800/1‑70 (MQ=255)

                          cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                    >  1:5597559/1‑70 (MQ=255)

                          cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                    >  1:5040003/1‑70 (MQ=255)

                           gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaaca                  <  1:2474110/71‑1 (MQ=255)

                           gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                   <  1:1425669/70‑1 (MQ=255)

                            cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                 <  1:5349487/71‑1 (MQ=34)

                            cACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                     <  1:3668574/67‑1 (MQ=255)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     <  1:1491225/49‑1 (MQ=38)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     <  1:4719608/49‑1 (MQ=38)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     <  1:2378503/49‑1 (MQ=38)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                >  1:2604141/1‑70 (MQ=35)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               <  1:5162549/71‑1 (MQ=35)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               <  1:1655642/71‑1 (MQ=35)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               >  1:720663/1‑71 (MQ=35)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               <  1:1923196/71‑1 (MQ=35)

                              cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               <  1:271100/71‑1 (MQ=35)

                              cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg               <  1:3593738/71‑2 (MQ=255)

                               ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                       <  1:3222575/46‑1 (MQ=38)

                               ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                           <  1:4416140/58‑1 (MQ=35)

                               ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                              <  1:1190031/55‑1 (MQ=35)

                               ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               >  1:5214529/1‑70 (MQ=37)

                               ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                     <  1:3713487/48‑1 (MQ=11)

                               ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                <  1:1609437/69‑1 (MQ=255)

                                 cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                              <  1:1627755/53‑1 (MQ=35)

                                 cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca             >  1:361626/1‑70 (MQ=40)

                                    aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt            <  1:1393065/68‑1 (MQ=40)

                                    aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt            <  1:1131736/68‑1 (MQ=40)

                                    aTACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                   <  1:1378131/61‑1 (MQ=255)

                                      aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                 >  1:2301972/1‑61 (MQ=255)

                                      aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                 >  1:827228/1‑61 (MQ=255)

                                       ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc      <  1:50855/71‑1 (MQ=255)

                                        cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCg     <  1:4109123/71‑1 (MQ=255)

                                         tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                <  1:3853896/59‑1 (MQ=35)

                                         tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                <  1:165367/59‑1 (MQ=35)

                                         tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                <  1:1947825/59‑1 (MQ=35)

                                         tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                <  1:837182/59‑1 (MQ=35)

                                         tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc    >  1:170487/1‑71 (MQ=255)

                                         tCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                        >  1:5323336/1‑35 (MQ=37)

                                         tCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                        >  1:2651291/1‑35 (MQ=37)

                                         tCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                        >  1:4084914/1‑35 (MQ=37)

                                          ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca             >  1:4063043/1‑61 (MQ=40)

                                          ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca             >  1:5630974/1‑61 (MQ=40)

                                          ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt   >  1:4547525/1‑71 (MQ=255)

                                           cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg  >  1:3524338/1‑71 (MQ=255)

                                           cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg  >  1:2333166/1‑71 (MQ=255)

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTG  >  Exported/2359818‑2359931

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