Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,125,913 G→T 17.9% G218W (GGG→TGG)  PP_5374 → choline/carnitine/betaine transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,125,9130GT17.9% 92.7 / 11.2 39G218W (GGG→TGG) PP_5374choline/carnitine/betaine transporter family protein
Reads supporting (aligned to +/- strand):  ref base G (15/17);  new base T (4/3);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

TGAGCATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTG  >  NC_002947/6125770‑6126056
                                                                                                                                               |                                                                                                                                                
tGAGCATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATAGGACATGTGGTGGACATCATTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCTCTTGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGa                                                                                                                                            <  2:439673/150‑1 (MQ=255)
  aGCATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGc                                                                                                                                          <  2:151835/150‑1 (MQ=255)
             gATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGa                                                                                                                               >  2:385323/1‑150 (MQ=255)
             gATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTCGACATCCTTGCCGTGGTGATTAGCGCATTTGGTGTATTCCAGTCGCTGGGGCTGGGCGTCGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGa                                                                                                                               <  2:171232/150‑1 (MQ=255)
                       ccGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCTGGCTCAGCCCGGTGttcatt                                                                                                                              >  1:153491/1‑137 (MQ=255)
                       ccGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCTGGCTCAGCCCGGTGttcatt                                                                                                                              <  2:153491/141‑5 (MQ=255)
                             aTTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGACTGAGCCAGGTGTTTGACCTGCCCATCAGCCTg                                                                                                               <  2:240860/150‑1 (MQ=255)
                                ggTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgggg                                                                                                            >  2:139176/1‑150 (MQ=255)
                                  tGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGt                                                                                                          <  2:328033/150‑1 (MQ=255)
                                            tATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCAcc                                                                                                <  2:150817/150‑1 (MQ=255)
                                                 gCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGc                                                                                                                                          <  2:426246/103‑1 (MQ=255)
                                                 gCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGc                                                                                                                                          >  1:426246/1‑103 (MQ=255)
                                                         gACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGGCCATCAGCCTGGGGGTGCAGATAATCTTAGATCggaaga                                                                                   >  2:97252/1‑145 (MQ=255)
                                                         gACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCa                                                                                  >  2:134769/1‑150 (MQ=255)
                                                         gACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCa                                                                                  <  1:15789/150‑1 (MQ=255)
                                                               tggtggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTAccc                                                                            >  1:404810/1‑150 (MQ=255)
                                                                 gtggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGTCTATCATCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTc                                                                          <  1:139176/150‑1 (MQ=255)
                                                                          cTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTAcc                                                                 >  2:493585/1‑150 (MQ=255)
                                                                              ccGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAgg                                                                                                                                      <  1:485590/78‑1 (MQ=255)
                                                                              ccGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAgg                                                                                                                                      >  2:485590/1‑78 (MQ=255)
                                                                                           cgcgTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCt                                                                                               >  1:282548/1‑104 (MQ=255)
                                                                                           cgcgTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCt                                                                                               <  2:282548/104‑1 (MQ=255)
                                                                                            gcgTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGc                                               <  1:385323/150‑1 (MQ=255)
                                                                                                          ccAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGt                                                                        >  1:292748/1‑111 (MQ=255)
                                                                                                          ccAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGt                                                                        <  2:292748/111‑1 (MQ=255)
                                                                                                           cAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTg                                                                                                         >  1:416162/1‑78 (MQ=255)
                                                                                                           cAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTg                                                                                                         <  2:416162/78‑1 (MQ=255)
                                                                                                           cAGTCGCTGGACCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgggg                                                                                                            >  2:58617/1‑75 (MQ=255)
                                                                                                           cAGTCGCTGGACCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgggg                                                                                                            <  1:58617/75‑1 (MQ=255)
                                                                                                                     gcctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTtacta                                                                   >  2:447837/1‑105 (MQ=255)
                                                                                                                     gcctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTtacta                                                                   <  1:447837/105‑1 (MQ=255)
                                                                                                                                     gtcaaaCACCTGGCTCAGCCCGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCgtgt                                           <  1:303473/109‑1 (MQ=255)
                                                                                                                                     gtcaaaCACCTGGCTCAGCCCGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCgtgt                                           >  2:303473/5‑113 (MQ=255)
                                                                                                                                        gAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCtgtt   >  2:76952/1‑150 (MQ=255)
                                                                                                                                         aaCACCTGGCTCAGCCCGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCtgttg  >  2:180137/1‑150 (MQ=255)
                                                                                                                                           cacCTGGCTCAGCCCGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGa                                                                                           >  2:295655/1‑59 (MQ=37)
                                                                                                                                           cacCTGGCTCAGCCCGGTGTTTGACCGGCCCATCAGCCTGGGGGTGCAGATCACCTTGa                                                                                           <  1:295655/59‑1 (MQ=25)
                                                                                                                                              cGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGaa                              <  2:418813/117‑1 (MQ=255)
                                                                                                                                              cGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGaa                              >  1:418813/1‑117 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TGAGCATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTG  >  NC_002947/6125770‑6126056

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