Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,331,783 G→A A428T (GCG→ACG)  proP → osmolyte:H(+) symporter ProP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,331,7830GA96.3% 83.4 / ‑4.8 28A428T (GCG→ACG) proPosmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (10/16);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTG  >  NC_000913/4331539‑4332022
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
tAACAGTAACGTCATCGGCCTGATTTTCTCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACTGGCCTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTTGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACggtg                                                                                                                                                                                                                                            <  1:140599/250‑1 (MQ=255)
            cATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAAATGCTTTACGGGCGGTATGGCTTCTACCTTGCCCGCGATGTTTCCGGCGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGGGCTGGTTGCCGGGCTGACGCCAACCCTGGCGGCCTGGCGTGG‑CGA‑AAGATTCCAGGATCTGGTGATGCCCGCCTATTTCCTGCTGGGAGGGGCggggggtgggtttatc                                                                                                                                                                                                                                >  2:228644/1‑236 (MQ=255)
                           tGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATTGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTTGGTGCTGGTTTCCGGTCTGACGCCAACGCCGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAAACTGATGATGCCTGCCCATTAACTGATGGGAGGGACGGGGGGTGGTTTAATCACCGGCGTAAccctg                                                                                                                                                                                                                 >  2:74997/1‑247 (MQ=255)
                             ccGGGTTTCTGATGCTGGCGGTGAGCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCGCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGaaa                                                                                                                                                                                                              <  1:203479/251‑1 (MQ=255)
                                   tACTGATGCTGGCGGTGATCCTTTACTGCTTTACGGGCGTTATGGCTTCTAACTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGGCTGACGCCAACGCTGGCGGCCTGGG‑TGGGCGA‑AAGCCCGCAGAAACTGATTATGCCTGCCCTTTACCTGATGGGTGTGACGGGGGGTGGGTTAATCAACGGCGTAACCATGAAagagg                                                                                                                                                                                                          >  2:138686/1‑248 (MQ=255)
                                          gCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAAt                                                                                                                                                                                                  <  1:74997/250‑1 (MQ=255)
                                                   gATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTg                                                                                                                                                                                         >  1:36940/1‑250 (MQ=255)
                                                                      ggCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCGAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCggcggc                                                                                                                                                                     <  1:57738/251‑1 (MQ=255)
                                                                                       cTTGCCAGCGATGTTACCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCAACGGCGTAACCATGAAAGAGACGGCAAATAGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGa                                                                                                                                                      >  1:194261/1‑249 (MQ=255)
                                                                                         tGCCAGCGATGTTCCATACGCATATCAGCTACAGCGCGCTGGAGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAAAGCTGGCGGCCTGGC‑TGGGCGATAAG‑TCGCAGAATCTGATGATGCCTGCCTATTACCTGATTGTAGTGTCGGCGGTTGGTTTAATCACCGGCGTGACCCTGAAAGCGATGGCAAATCGTACGATGAAAGGTGCGACCCCGACGGCGTCACGTATACAGGAAga                                                                                                                                                   >  2:142844/1‑249 (MQ=255)
                                                                                            cAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGtt                                                                                                                                                                                          <  1:244755/208‑1 (MQ=255)
                                                                                            cAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGtt                                                                                                                                                                                          >  2:244755/1‑208 (MQ=255)
                                                                                             aGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAgg                                                                                                                                              <  1:138686/251‑1 (MQ=255)
                                                                                                      cccGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATtctc                                                                                                                                      <  1:140967/250‑1 (MQ=255)
                                                                                                                atCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCAtt                                                                                                                            <  1:188835/250‑1 (MQ=255)
                                                                                                                atCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTa                                                                                                                           <  1:87819/251‑1 (MQ=255)
                                                                                                                    gTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGat                                                                                                                       >  1:270659/1‑251 (MQ=255)
                                                                                                                       aCAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGataat                                                                                                                    <  2:48382/251‑1 (MQ=255)
                                                                                                                                         cATTTAATATTTCGGGGCTTGTTGCCGGTCTGACGCCAACGATGGCGGGCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCTAGCAGAAAATCga                                                                                                   <  2:199748/250‑1 (MQ=255)
                                                                                                                                             tAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCgatgat                                                                                               <  2:271478/250‑1 (MQ=255)
                                                                                                                                                          cTGGTTGCCAGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACg                                                                                     >  2:201872/1‑247 (MQ=255)
                                                                                                                                                            tgTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGAtt                                                                                <  2:36940/249‑1 (MQ=255)
                                                                                                                                                                        gACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGc                                                                      <  2:57241/248‑1 (MQ=255)
                                                                                                                                                                         aCGCCAACGCTGGCGGCCTGGC‑TGGGCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGc                                                                  <  1:263234/251‑1 (MQ=255)
                                                                                                                                                                                                        aaGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATc                                     >  1:110117/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aTTACCTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTGGAAGGGTGCGACCCCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTGGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGa         <  1:173509/250‑1 (MQ=255)
                                                                                                                                                                                                                                          cTGATGGTAGTGACGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCt   >  1:132946/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gATGGTAGTGACGGTGGTTGGTTTAATCCCCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTg  >  1:107668/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
TAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGC‑TGGTCGA‑AAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTG  >  NC_000913/4331539‑4332022

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