Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,321,803 T→C H495H (CAT→CAC proP → glycine betaine/L‑proline transporter ProP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,321,8030TC96.4% 85.3 / ‑3.3 28H495H (CAT→CACproPglycine betaine/L‑proline transporter ProP
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (10/17);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-01

AGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTG  >  NZ_CP009273/4321688‑4321930
                                                                                                                   |                                                                                                                               
aGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTg                                                                                                           >  2:81401/1‑138 (MQ=255)
     gCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACgg                                                                                                     <  2:40758/139‑1 (MQ=255)
     gCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACg                                                                                                      >  2:249597/1‑138 (MQ=255)
                    gtcgAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGAc                                                                                      <  1:249597/139‑1 (MQ=255)
                    gtcgAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGAc                                                                                      <  2:150556/139‑1 (MQ=255)
                     tcgAGCATTACGATAATATCGAGCAGAAAACCGATGATATTGACCACGAGGTTGCGGATTTGCAGGCGAAAGGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACg                                                                                     <  2:326316/139‑1 (MQ=255)
                            ttACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCt                                                                               >  1:329746/1‑138 (MQ=255)
                                    atatCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGc                                                                       <  2:329746/138‑1 (MQ=255)
                                     tatCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACg                                                                                     >  1:204084/1‑123 (MQ=255)
                                     tatCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACg                                                                                     <  2:204084/123‑1 (MQ=255)
                                         gAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTg                                                                 <  1:81401/139‑1 (MQ=255)
                                          cgCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTg                                                                 >  1:157636/2‑138 (MQ=255)
                                          aGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTg                                                                 >  1:431374/1‑138 (MQ=255)
                                           gCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGaa                                                               <  2:327547/139‑1 (MQ=255)
                                                 aaTCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCtt                                                          >  2:383841/1‑138 (MQ=255)
                                                           atTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACt                                                >  1:142848/1‑138 (MQ=255)
                                                                   cGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGGGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATctct                                        >  2:249626/1‑138 (MQ=255)
                                                                    gagaTTGCCGATTTGCAGGCGAACCGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGACCTCTTAAGGTTCACTTAATCTCTGa                                      <  1:249626/139‑1 (MQ=255)
                                                                      gaTTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGAcg                                    <  2:282898/139‑1 (MQ=255)
                                                                       aTTGCCGATGTGCAGCCGAAACGTGCCCGCCTGGTGCAGCATCACCCGCGAATTGATGAATAATGTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAAGTGAACCCTTAAGGTTCACTTAAGCTCTTAcgc                                   <  2:232245/139‑1 (MQ=255)
                                                                                        gAAACGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAgg                                                      >  2:417857/1‑103 (MQ=255)
                                                                                        gAAACGTACCAGCCTGGTGCAGCAACACCCTCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAgg                                                      <  1:417857/103‑1 (MQ=255)
                                                                                         aaaCGTACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGt                 <  2:125151/139‑1 (MQ=255)
                                                                                              tACCCGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACg            <  2:7688/139‑1 (MQ=255)
                                                                                                cccGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACg            >  2:211102/1‑137 (MQ=255)
                                                                                                  cGCCTGGTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACggagg        <  1:211102/139‑1 (MQ=255)
                                                                                                        gTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTg  <  2:142848/139‑1 (MQ=255)
                                                                                                        gTGCAGCAACACCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTg  <  2:74075/139‑1 (MQ=255)
                                                                                                                   |                                                                                                                               
AGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTG  >  NZ_CP009273/4321688‑4321930

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