Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,993,464 Δ1 bp coding (2312/2547 nt) cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,993,4640A.89.7% 90.7 / 8.8 30coding (2312/2547 nt)cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base . (13/13);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGA  >  NC_000913/3993326‑3993589
                                                                                                                                          |                                                                                                                             
aTCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTGTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcacca                                                                                                                               >  1:370416/1‑139 (MQ=255)
                                 ttGCCTGCTAAGGGATCATCCTGTTCTTTTTCGAAGAAACGCAAGACGAGATTGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGAc                                                                                             <  1:611913/139‑1 (MQ=255)
                                  tGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACg                                                                                            >  1:399454/1‑139 (MQ=255)
                                         gAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCg                                                                                     >  2:97582/1‑139 (MQ=255)
                                         gAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCg                                                                                     <  2:399454/139‑1 (MQ=255)
                                                taatcctattCTTTCACGAAGAAACACAAGACGAGAATGGCTTTAATATCTACATTCGCGACGAAAGCAACCGGGTAGAGGTAAATCACCCTGCGAAGGCAGCAAAGACGAGCTAGTACGTGACGGCAGTCGCTTCTAc                                                                              <  1:220942/131‑1 (MQ=255)
                                                    ccAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGt                                                                          <  1:504690/139‑1 (MQ=255)
                                                      aGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGTAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCa                                                                        <  2:421122/139‑1 (MQ=255)
                                                        ttCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATc                                                                      >  1:516893/1‑139 (MQ=255)
                                                              ttCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATg                                                                 <  2:370416/139‑1 (MQ=255)
                                                               tCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGa                                                                <  2:394023/139‑1 (MQ=255)
                                                                       aCGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAgg                                                                                                                      <  1:293813/76‑1 (MQ=255)
                                                                       aCGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAgg                                                                                                                      >  2:293813/1‑76 (MQ=255)
                                                                          cAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCt                                                     <  2:208620/139‑1 (MQ=255)
                                                                          cAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTa                                                    <  1:346530/139‑1 (MQ=255)
                                                                             gACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACgg                                                 >  1:516162/1‑139 (MQ=255)
                                                                               cGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCt                                               >  1:157375/1‑139 (MQ=255)
                                                                                   aaTGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAg                                           >  1:488446/1‑139 (MQ=255)
                                                                                        cTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGa                                                                <  1:22724/113‑1 (MQ=255)
                                                                                        cTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGa                                                                >  2:22724/1‑113 (MQ=255)
                                                                                                 cTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGt                                                                                       <  2:352023/81‑1 (MQ=255)
                                                                                                 cTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGt                                                                                       >  1:352023/1‑81 (MQ=255)
                                                                                                 cTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCa                             >  2:494978/1‑139 (MQ=255)
                                                                                                 cTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCa                             >  2:575656/1‑139 (MQ=255)
                                                                                                         tcGACGAATGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCg                     >  2:400053/1‑139 (MQ=255)
                                                                                                            acgaAAGCAACCGGGTTGAGGTATATCACCCTGCGAAGGCAGCTAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAg                  <  1:400053/139‑1 (MQ=255)
                                                                                                                    aaCCGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCa          <  1:494978/139‑1 (MQ=255)
                                                                                                                      ccGGGTTGAGGTATATCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGa        <  2:516893/139‑1 (MQ=255)
                                                                                                                            tGAGGTAAATCCCCCTGCGAAGGCAGCAAAGAAGAGCTGGTACGGTACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGa  <  2:563064/139‑1 (MQ=255)
                                                                                                                                   tatCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCa                                             >  1:459745/1‑89 (MQ=255)
                                                                                                                                   tatCACCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCa                                             <  2:459745/89‑1 (MQ=255)
                                                                                                                                          |                                                                                                                             
ATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGA  >  NC_000913/3993326‑3993589

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: