Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,234,667 Δ1 bp intergenic (‑1434/‑499) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,234,6660A.96.8% 125.6 / ‑3.2 31intergenic (‑1433/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (13/17);  minor base C (1/0);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTG  >  NC_000913_3_pae_pgi/4234520‑4234810
                                                                                                                                                   |                                                                                                                                                
tCTATATACGCCGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                                >  2:154041/1‑147 (MQ=255)
    tataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcg                                                                                                                                           <  1:154041/149‑1 (MQ=255)
      taCGTAGTAGAGATGGACCCCAGGGGCGTTTTATGCTTTCTACGGGATGTGGGCATCGCCTCTGTAAAACTGGCAGCGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcgtc                                                                                                                                         <  2:116863/149‑1 (MQ=255)
         gtaGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTGACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGat                                                                                                                                      <  2:106318/147‑1 (MQ=255)
         gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGa                                                                                                                                                  >  1:328674/1‑138 (MQ=255)
         gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGa                                                                                                                                                  <  2:328674/138‑1 (MQ=255)
         gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGat                                                                                                                                      >  2:289400/1‑149 (MQ=255)
             tAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                                >  1:237021/1‑134 (MQ=255)
             tAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                                <  2:237021/136‑3 (MQ=255)
               gagaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAACCTGGCAGAGATGA‑GGGGCGTAAATCCCCTGGCCACCCCTAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGGTATTTAg                                                                                                                                >  2:333146/1‑149 (MQ=255)
                  aTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGcc                                                                                                                             >  1:281608/1‑149 (MQ=255)
                       ccccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATa                                                                                                                        >  2:380249/1‑149 (MQ=255)
                           aaGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                    >  1:9804/1‑149 (MQ=255)
                                cGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                                                                               <  2:143741/149‑1 (MQ=255)
                                    ttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCTAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCCcgcc                                                                                                           <  1:102114/149‑1 (MQ=255)
                                     tATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACCCACCTA‑ATGAGGATGCCTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTGGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGGTAAcgccg                                                                                                          >  1:98106/1‑149 (MQ=255)
                                           ttttACGAGATGTGGGCATCGCTCCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                    <  1:12798/149‑1 (MQ=255)
                                               aCGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAAACCCCTCGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                                <  1:95418/149‑1 (MQ=255)
                                                cGAGATGTGGGCATCGTCTTTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAGCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                               <  1:352172/149‑1 (MQ=255)
                                                cGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGACCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                               >  2:68190/1‑149 (MQ=255)
                                                               gCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGa                                                                                <  1:225529/149‑1 (MQ=255)
                                                                                agggaTTA‑GCGGCGTAAATCCCCTGGCCAACACGAGAATAGATGCATTGAAGAGGCGCTTGCACTG‑ATCATCGTCGAGATTTAGGCCGGATAAGTCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                               <  2:98106/146‑1 (MQ=255)
                                                                                            gCGTAAATACCCGGGCCACCCCGAGAATAGATGCAGTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGCATAAGGCGTTCACGCCGCATGTGGAAACCGATGCCTTATAAGACGCGCCCGCGTCGCATCAGGcatcggtt                                                    <  2:69090/149‑8 (MQ=255)
                                                                                                       cTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                         <  1:380249/149‑1 (MQ=255)
                                                                                                        tgcccACACCAAGAATAGAGCGATTTAAGAGGCGGTAGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                        <  1:131442/146‑1 (MQ=21)
                                                                                                        tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCCTCGGCGAAATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCGGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                        >  1:388804/1‑149 (MQ=25)
                                                                                                        tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                        <  2:9804/149‑1 (MQ=255)
                                                                                                         ggCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGAGAAGCCGTTCACGCCGCATACGGCAACCGTTGCCTGATGCACCGCGGTCGCGTCTTGATAGGCCTACCGGTTCATGCCGatat                                       >  2:122041/1‑149 (MQ=255)
                                                                                                          gCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGT‑CATCATCGTCGATATATAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCTTGATGCGACGCGGTCGCGTCTTATCAGTCTTACAGGTCTATGCCGATATg                                      >  2:391557/1‑149 (MQ=37)
                                                                                                                      aTAGATGCAGTGAAGAGGCGCTTGCAGTG‑AGCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATTCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                          <  1:46503/149‑1 (MQ=255)
                                                                                                                          aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                      >  1:205904/1‑149 (MQ=255)
                                                                                                                          aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                      >  2:401088/1‑149 (MQ=255)
                                                                                                                              aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACTGATGCCTGATGCGACGCGGTCGCGTGTTGTCAGGCGTACAGGTCGATGCAGATATGTGTATCGGATGCGGCAAGTa                  >  1:121846/1‑149 (MQ=25)
                                                                                                                               ttGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa                 <  1:68190/149‑1 (MQ=255)
                                                                                                                                 gAAGAGGCGCTTGCAGTG‑ATCATCGTCGAGATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCTACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAGCGTATTCGGCAATTattc               <  2:246030/149‑4 (MQ=25)
                                                                                                                                              cAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTAGATCGTTTTCGGCACTTAATACATAGCACgattg  >  1:75889/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTA‑GAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTG  >  NC_000913_3_pae_pgi/4234520‑4234810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑