Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,518,526 T→C 16.0% E183E (GAA→GAG spoT ← bifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,518,5260TC16.0% 191.4 / 26.8 81E183E (GAA→GAGspoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (25/43);  new base C (2/11);  total (27/54)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

CGCGCGGCTTTCACCACTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAATTTCGAGAG  >  Exported/2518458‑2518590
                                                                    |                                                                
cgcgcgGCTTTCACCACTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTc                                                                <  1:3239522/71‑1 (MQ=255)
 gcgcgGCTTTCACCACTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTc                                                                <  1:676921/70‑1 (MQ=255)
       cTTTCACCACTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCg                                                         <  1:4621371/71‑1 (MQ=255)
       cTTTCACCACTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCg                                                         >  1:2165083/1‑71 (MQ=255)
       cTTTCACCACTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCg                                                         >  1:2403398/1‑71 (MQ=255)
           caccacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGttt                                                     >  1:4139813/1‑71 (MQ=255)
           caccacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGttt                                                     >  1:3024704/1‑71 (MQ=255)
            accacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGa                                                              >  1:410222/1‑61 (MQ=255)
            accacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGa                                                              >  1:1365010/1‑61 (MQ=255)
            accacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCGGtttt                                                    <  1:2318656/71‑1 (MQ=255)
            accacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCGGtttt                                                    <  1:4778043/71‑1 (MQ=255)
            accacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCGGtttt                                                    <  1:1742169/71‑1 (MQ=255)
            accacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCGGtttt                                                    <  1:2774278/71‑1 (MQ=255)
            accacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCGGtttt                                                    <  1:2832962/71‑1 (MQ=255)
            accacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCGGtttt                                                    <  1:5187496/71‑1 (MQ=255)
            accacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCGGtttt                                                    <  1:3290205/71‑1 (MQ=255)
            accacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCGGtttt                                                    <  1:141482/71‑1 (MQ=255)
            accacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCGGtttt                                                    <  1:3554772/71‑1 (MQ=255)
             ccacTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCGGtttt                                                    <  1:1110540/70‑1 (MQ=255)
                cttcttTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAAtg                                                >  1:4675920/1‑71 (MQ=255)
                cttcttTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAAtg                                                >  1:4658690/1‑71 (MQ=255)
                     ttGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTaa                                                  >  1:4683500/1‑64 (MQ=255)
                     ttGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTaa                                                  >  1:4678700/1‑64 (MQ=255)
                     ttGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTaa                                                  >  1:4929056/1‑64 (MQ=255)
                     ttGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTaa                                                  >  1:4279353/1‑64 (MQ=255)
                     ttGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATgtggtg                                           >  1:351391/1‑71 (MQ=255)
                     ttGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATgtggtg                                           >  1:677965/1‑71 (MQ=255)
                     ttGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATgtggtg                                           >  1:4441799/1‑71 (MQ=255)
                            cgcgATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAAtgtg                                              >  1:5074782/1‑61 (MQ=255)
                            cgcgATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATgtgg                                             >  1:3852888/1‑62 (MQ=255)
                            cgcgATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATgtgg                                             >  1:2330823/1‑62 (MQ=255)
                             gcgATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTa                                   <  1:4162005/71‑1 (MQ=255)
                             gcgATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTa                                   <  1:3507849/71‑1 (MQ=255)
                             gcgATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTa                                   <  1:4088610/71‑1 (MQ=255)
                             gcgATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTa                                   <  1:3384871/71‑1 (MQ=255)
                             gcgATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTa                                   <  1:3240786/71‑1 (MQ=255)
                             gcgATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTa                                   <  1:2164445/71‑1 (MQ=255)
                              cgATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTa                                   >  1:4171455/1‑70 (MQ=255)
                              cgATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTa                                   >  1:3502028/1‑70 (MQ=255)
                                aTAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGt                                                            <  1:5167706/43‑1 (MQ=255)
                                aTAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGt                                                            <  1:583651/43‑1 (MQ=255)
                                aTAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGt                                                            <  1:5135791/43‑1 (MQ=255)
                                aTAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGt                                                            <  1:4291802/43‑1 (MQ=255)
                                aTAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGt                                                            <  1:5014607/43‑1 (MQ=255)
                                aTAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGt                                                            <  1:4937090/43‑1 (MQ=255)
                                aTAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGt                                                            <  1:3781126/43‑1 (MQ=255)
                                aTAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGt                                                            <  1:368524/43‑1 (MQ=255)
                                aTAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGt                                                            <  1:225921/43‑1 (MQ=255)
                                aTAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGt                                                            <  1:2256882/43‑1 (MQ=255)
                                  aaCGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACgg                              >  1:708977/1‑71 (MQ=255)
                                   aCGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACgg                              <  1:1349125/70‑1 (MQ=255)
                                   aCGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACgg                              <  1:3032270/70‑1 (MQ=255)
                                         gggATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACgg                              <  1:954308/64‑1 (MQ=255)
                                         gggATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACgg                              <  1:2692623/64‑1 (MQ=255)
                                         gggATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACgg                              <  1:4433652/64‑1 (MQ=255)
                                         gggATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACgg                              <  1:5116314/64‑1 (MQ=255)
                                         gggATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACgg                              <  1:4787572/64‑1 (MQ=255)
                                          ggATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACgg                              <  1:1394430/63‑1 (MQ=255)
                                          ggATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      <  1:3928325/71‑1 (MQ=255)
                                          ggATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      <  1:4813196/71‑1 (MQ=255)
                                          ggATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGGTTTAATGTGGTGGATACCTAAACgg                              <  1:1595611/63‑1 (MQ=255)
                                           gATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACgg                              <  1:3318573/62‑1 (MQ=255)
                                           gATACAGCGCCTCAAAACCCAGCTCCTCGAGTTCGGtttt                                                    <  1:5104129/40‑1 (MQ=255)
                                               cAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGc                 >  1:4498343/1‑71 (MQ=255)
                                                   gcCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCtata             <  1:1420732/71‑1 (MQ=255)
                                                   gcCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCtata             <  1:3037910/71‑1 (MQ=255)
                                                   gcCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCtata             <  1:4362886/71‑1 (MQ=255)
                                                      tCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTa                                   <  1:3219926/46‑1 (MQ=255)
                                                      tCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTa                                   <  1:2841326/46‑1 (MQ=255)
                                                      tCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTa                                   <  1:2702336/46‑1 (MQ=255)
                                                        aaaaCCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAAttt        >  1:1705523/1‑71 (MQ=255)
                                                        aaaaCCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAAttt        >  1:3952897/1‑71 (MQ=255)
                                                        aaaaCCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAAtt         >  1:2313538/1‑70 (MQ=255)
                                                        aaaaCCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAAt          <  1:4629320/69‑1 (MQ=255)
                                                        aaaaCCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAAt          >  1:4732018/1‑69 (MQ=255)
                                                        aaaaCCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAAt          <  1:1665749/69‑1 (MQ=255)
                                                        aaaaCCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAAt          <  1:4924957/69‑1 (MQ=255)
                                                          aaCCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAATTTCg      <  1:3948752/71‑1 (MQ=255)
                                                          aaCCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAATTTCg      <  1:70544/71‑1 (MQ=255)
                                                           aCCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAATTTCga     <  1:4618131/71‑1 (MQ=255)
                                                            cccAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAATTTc       >  1:4412315/1‑68 (MQ=255)
                                                                   ctTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCATCGGGCTATAAATTTc       <  1:4593598/61‑1 (MQ=255)
                                                                   ctTCGAGTTCGGTTTTAATGTGGGGGATACCTAAACGGTGCGCCAGCGGGCTATAAATTTc       <  1:245229/61‑1 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:1830498/1‑45 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:4979263/1‑45 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:169150/1‑45 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:2217517/1‑45 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:3675197/1‑45 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:346978/1‑45 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:635020/1‑45 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:4064690/1‑45 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:3790893/1‑45 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:1038846/1‑45 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:3763780/1‑45 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAATTTCGagag  >  1:3078951/1‑65 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAATTTCGagag  >  1:1651608/1‑65 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAATTTCGagag  >  1:1637542/1‑65 (MQ=255)
                                                                    ttCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAATTTCGagag  >  1:4084540/1‑65 (MQ=255)
                                                                    ttCGAGTTCGGTATTAATGTGGTGGATACCTAAACGGTGCGCCAg                      >  1:363782/1‑45 (MQ=255)
                                                                    |                                                                
CGCGCGGCTTTCACCACTTCTTTGATTACGCGATAACGGTTGGGATACAGCGCCTCAAAACCCAGCTCTTCGAGTTCGGTTTTAATGTGGTGGATACCTAAACGGTGCGCCAGCGGGCTATAAATTTCGAGAG  >  Exported/2518458‑2518590

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: