Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,331,890 C→T L463L (CTC→CTT proP → osmolyte:H(+) symporter ProP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,331,8900CT95.0% 44.1 / ‑3.7 20L463L (CTC→CTTproPosmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (5/14);  total (6/14)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

CGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATG  >  NC_000913/4331643‑4332134
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
cGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCTTCAGATATACAGGAAGCGAAGGAAATTCttgtc                                                                                                                                                                                                                                                   <  1:253721/251‑1 (MQ=255)
               acggcgcgcTGGTGGCGGCATTTAATATTTCGGTGCTGGTTGCCTGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCTGCGTCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGTGCATTACGataa                                                                                                                                                                                                                                     <  1:294150/247‑1 (MQ=255)
                                  aTTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCgatg                                                                                                                                                                                                                 <  1:287829/251‑1 (MQ=255)
                                    ttAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCgatga                                                                                                                                                                                                                <  1:276728/250‑1 (MQ=255)
                                                                                                     gagaaggtgagtatgCCTGCCTATTACCGGAGGGTAGTGGGGTGGGTTGGTTTAATCACCGGTGTATCCATGAAAGAGAGGGCAAATCGTCCGTTGAAAGGGGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAAAAGCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAAt                                                                                                                                              <  1:173280/239‑1 (MQ=255)
                                                                                                                          tATTACCGGAGGGTAGTGGCGGGGGTTGGTTTAATCACCGGCGTAACCATGAAAGTGACGGCAAACCGTCCGTGGAAAGGTGCGACACCGGCGGCTTCAGATAACCAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAACCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACg                                                                                                                           <  2:36308/249‑1 (MQ=255)
                                                                                                                                      gTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGAtgtg                                                                                                             <  1:93175/251‑1 (MQ=255)
                                                                                                                                            gCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTg                                                                                                       <  1:31517/251‑1 (MQ=255)
                                                                                                                                                       gttTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGCTATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGGTTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAgg                                                                                             >  2:139667/1‑251 (MQ=255)
                                                                                                                                                                           aTGAAAGAGCCGGAAAACCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATAAAGGAAGCGAAGGAAATTCTTGTCGAGAATTACGATAATATCGAGCAGAAAATCGATGATATTGACTACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACTGATGGCCTGATGTGACGCTGTCGTATCAGGCCAATTGAACTCTTAAGGtt                                                                         <  1:234071/251‑1 (MQ=255)
                                                                                                                                                                            tGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGGAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTc                                                                        >  1:147366/1‑251 (MQ=255)
                                                                                                                                                                                              gTCCGTTGAAAGGTGCGACACCGGCGGAGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGACACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACgcgc                                                      >  2:216612/1‑251 (MQ=255)
                                                                                                                                                                                                            gcgaacccGGCGGGGTCAGAATTAAAGGAAGAGAAGGAAATTCTTGTGGAGCTTTACGATAATATCGAGCAGAAAAGCGATGATATTGACCTCGAGATTGCCGATGTGCAGGCGAAACGTACCCGCCGGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACGTAATCTCTGACGCGCATACTCTCCTCCa                                         <  2:254499/244‑1 (MQ=255)
                                                                                                                                                                                                                cacCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGtt                                     >  1:29951/1‑250 (MQ=255)
                                                                                                                                                                                                                   cggcggCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACg                                 <  1:256189/251‑1 (MQ=255)
                                                                                                                                                                                                                       gggcgcagttataCAGGAAGCGAAGAAAATTCTTGTCGAGCAGTACGATAATATCTGGCAGAAAATCGATGATATGGACCACGAGATTTCCGTTTTGCAGGCTAAACTTACCCGCCTGGGGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTACCAGCCCAATTGAACTCTTAAGGTGTACTTAATCTCAGACGCGCATACTCTCCTCCAGGTGAACggagg                             <  2:296075/242‑1 (MQ=255)
                                                                                                                                                                                                                       ggcgTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGCCCACGAGATTGCCGATTTGCAGGCGAAACTTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACggagg                             <  1:239187/251‑1 (MQ=255)
                                                                                                                                                                                                                             aGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTg                       >  2:255516/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gatgCGAAGGAAATTCTTGTCGAGCATTACGTTAATTTCGAGAAGAAAATCGTTGATATTGCCCCCGAGATTGCCGATTTGCAGGCGAAACGTCCCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGTAGAGTGCAATGaaaaa             <  1:242959/248‑1 (MQ=255)
                                                                                                                                                                                                                                              aagAAATTCTTGTCGAGCATTACGAAAATATCGAGCAGAAAATCGATGATATGGACCACGAGATGGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCGCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGtt      <  2:124861/249‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaTTCTTGTCGAGCATTACGATAATATCGTGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACTGAGGAGAGTGCAATGAAAAACCGTGTTTATg  >  2:48979/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
CGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATG  >  NC_000913/4331643‑4332134

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