Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,518,856 G→A 7.2% D73D (GAC→GAT 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,8560GA 7.2% 272.9 / 8.7 83D73D (GAC→GATspoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (46/31);  new base A (4/2);  total (50/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTA  >  Exported/2518787‑2518926
                                                                     |                                                                      
cAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                                       >  1:2977191/1‑71 (MQ=255)
 aGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                                       >  1:2964796/1‑70 (MQ=255)
  gCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCa                                                                     <  1:1435286/71‑1 (MQ=255)
  gCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCa                                                                     <  1:2189596/71‑1 (MQ=255)
    tCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATg                                                                   >  1:3332155/1‑71 (MQ=255)
     cGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATg                                                                   >  1:2715722/1‑70 (MQ=255)
      ggCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATgca                                                                 <  1:605792/71‑1 (MQ=255)
        cGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATgcag                                                                >  1:2694040/1‑70 (MQ=255)
           cGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATgca                                                                 >  1:3796510/1‑66 (MQ=255)
           cGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCagcagc                                                            >  1:828202/1‑71 (MQ=255)
           cGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCagcagc                                                            >  1:52649/1‑71 (MQ=255)
            gCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAgcg                                                           >  1:940297/1‑71 (MQ=255)
                   ccAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCAt                                                    >  1:1009818/1‑71 (MQ=255)
                   ccAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCAt                                                    >  1:4001146/1‑71 (MQ=255)
                   ccAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCAt                                                    >  1:3928624/1‑71 (MQ=255)
                               tCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTc                                           >  1:1291767/1‑68 (MQ=255)
                                ccATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAg                                       >  1:355024/1‑71 (MQ=255)
                                 cATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACg                                                                        >  1:4266451/1‑37 (MQ=255)
                                  atatCCTGGTATGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTc                                     <  1:2124787/71‑1 (MQ=255)
                                  atatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                                       >  1:3118645/1‑37 (MQ=255)
                                  atatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                                       >  1:436797/1‑37 (MQ=255)
                                  atatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                                       >  1:1637602/1‑37 (MQ=255)
                                  atatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                                       >  1:4058352/1‑37 (MQ=255)
                                  atatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGt                                                                       >  1:2769226/1‑37 (MQ=255)
                                   tatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTc                                     >  1:3861253/1‑70 (MQ=255)
                                   tatCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTc                                     >  1:3261600/1‑70 (MQ=255)
                                     tCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATgcagca                                                              >  1:728838/1‑43 (MQ=255)
                                     tCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATgca                                                                 >  1:2172637/1‑40 (MQ=255)
                                     tCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGcagcag                                                             >  1:2469247/1‑44 (MQ=255)
                                     tCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGcagcag                                                             >  1:1308013/1‑44 (MQ=255)
                                      ccTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGt                                 >  1:3383147/1‑71 (MQ=255)
                                      ccTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGt                                 >  1:471984/1‑71 (MQ=255)
                                      ccTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGt                                 >  1:633174/1‑71 (MQ=255)
                                      ccTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGt                                 >  1:492457/1‑71 (MQ=255)
                                      ccTGGTAGGTGGAGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGt                                 >  1:2738902/1‑71 (MQ=255)
                                       cTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGtt                                <  1:945360/71‑1 (MQ=255)
                                       cTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGtt                                <  1:3573452/71‑1 (MQ=255)
                                        tGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCg                                               <  1:1956127/55‑1 (MQ=255)
                                        tGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCg                                               <  1:3272235/55‑1 (MQ=255)
                                        tGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCg                                               <  1:2022786/55‑1 (MQ=255)
                                        tGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGtt                                >  1:3466994/1‑70 (MQ=255)
                                          gTAGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCg                                               <  1:3813527/53‑1 (MQ=39)
                                          gTAGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCg                                               <  1:2014347/53‑1 (MQ=39)
                                          gTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGttt                               <  1:4276813/69‑1 (MQ=255)
                                          gTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGttt                               <  1:3314495/69‑1 (MQ=255)
                                          gTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGttt                               <  1:2595588/69‑1 (MQ=255)
                                          gTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTc                              <  1:2848010/70‑1 (MQ=255)
                                            aGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCGtt                                             >  1:3643209/1‑53 (MQ=39)
                                            aGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCGtt                                             >  1:985662/1‑53 (MQ=39)
                                            aGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCGtt                                             >  1:3331001/1‑53 (MQ=39)
                                            aGGTGGCGGGTGTATCTTCAATCACATCATGCAGCAGCGCCGCCATCAGCGtt                                             >  1:3182362/1‑53 (MQ=39)
                                            aGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCAtc                           >  1:3493867/1‑71 (MQ=255)
                                             ggTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCAtct                          >  1:2268266/1‑71 (MQ=255)
                                               tGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGttt                               <  1:904977/64‑1 (MQ=255)
                                               tGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGttt                               <  1:3811393/64‑1 (MQ=255)
                                                ggCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCgg                       >  1:3122182/1‑71 (MQ=255)
                                                 gCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCa                             >  1:716077/1‑64 (MQ=255)
                                                 gCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCa                             >  1:2512466/1‑64 (MQ=255)
                                                 gCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTAGAGTTTCa                             >  1:1914530/1‑64 (MQ=255)
                                                    ggAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAg                   <  1:2001533/71‑1 (MQ=255)
                                                       gTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCa                    <  1:239353/67‑1 (MQ=255)
                                                       gTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCa                    <  1:892558/67‑1 (MQ=255)
                                                       gTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCa                    <  1:3061575/67‑1 (MQ=255)
                                                       gTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCa                    <  1:3176079/67‑1 (MQ=255)
                                                       gTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAAt                <  1:1700880/71‑1 (MQ=255)
                                                        tATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGt                                 <  1:1748317/53‑1 (MQ=255)
                                                        tATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATg               >  1:353537/1‑71 (MQ=255)
                                                        tATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATg               <  1:698745/71‑1 (MQ=255)
                                                          tCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGttt                               >  1:1365264/1‑53 (MQ=255)
                                                          tCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTc                              <  1:2796675/54‑1 (MQ=255)
                                                          tCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTc                              <  1:154886/54‑1 (MQ=255)
                                                          tCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTc                              <  1:2425740/54‑1 (MQ=255)
                                                          tCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATg               >  1:60134/1‑69 (MQ=255)
                                                          tCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGc              <  1:174514/70‑1 (MQ=255)
                                                          tCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGc              <  1:382476/70‑1 (MQ=255)
                                                           cTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAAt                >  1:1109435/1‑67 (MQ=255)
                                                              cAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAg            <  1:2034555/68‑1 (MQ=255)
                                                               aaTCACGTCATGCAGCAGCGCCGCCATCAGCGttt                                            >  1:3697887/1‑35 (MQ=255)
                                                               aaTCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCa         <  1:3772184/70‑1 (MQ=255)
                                                               aaTCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCa         <  1:3429017/70‑1 (MQ=255)
                                                                   aCGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCg                                    >  1:1330608/1‑39 (MQ=255)
                                                                   aCGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCa                             >  1:50352/1‑46 (MQ=255)
                                                                   aCGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCa                             >  1:316596/1‑46 (MQ=255)
                                                                   aCGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGc    >  1:3623092/1‑71 (MQ=255)
                                                                     gTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTa  >  1:753324/1‑71 (MQ=255)
                                                                     gTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTa  >  1:3086410/1‑71 (MQ=255)
                                                                     gTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTa  >  1:2774090/1‑71 (MQ=255)
                                                                     |                                                                      
CAGCTCGGCGACGCTTTTACCAAAAAGCTGTTCCATATCCTGGTAGGTGGCGGGAGTATCTTCAATCACGTCATGCAGCAGCGCCGCCATCAGCGTTTCATAGTCGAGTTTCATCTCGGCCAGAATGCAGGCAACCGCTA  >  Exported/2518787‑2518926

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: