Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,141,748 T→C 34.1% V12A (GTG→GCG) ‡ PP_4526 → MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,141,7480TC34.1% 28.7 / 18.6 29V12A (GTG→GCG) ‡PP_4526MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (11/8);  new base C (5/5);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GGCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTGGGCGCGGTGGAC  >  NC_002947/5141610‑5141893
                                                                                                                                          |                                                                                                                                                 
ggCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCACGGCGCtggcga                                                                                                                                        <  2:326619/150‑5 (MQ=255)
          aaaaCAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGCCCATggc                                                                                                                              >  2:331490/1‑149 (MQ=255)
          aaaaCAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCACGGCGCtggc                                                                                                                                          >  2:510429/1‑136 (MQ=255)
          aaaaCAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCACGGCGCtggc                                                                                                                                          <  1:510429/138‑3 (MQ=255)
                    aacTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGCCCATGGCCATcatcagc                                                                                                                    <  1:331490/150‑7 (MQ=255)
                          atcatcATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCACGGCGCtggcga                                                                                                                                        <  1:141664/124‑5 (MQ=255)
                          atcatcATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCACGGCGCtggcga                                                                                                                                        >  2:141664/1‑120 (MQ=255)
                                         ttttATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCccgcaactg                                                                                               >  1:268710/1‑150 (MQ=255)
                                               tGAAAAATAAACCCCCCTGACCTGCAATCGATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTTATCCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCccgcaactgccgcaa                                                                                         <  1:215794/150‑1 (MQ=255)
                                                gAAAAATAAACCGCCCTAACCTGCAATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCccgcaactgccgcaac                                                                                        >  1:387873/1‑150 (MQ=255)
                                                       aaaCCGCCTTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGc                                                                                 <  2:497657/150‑1 (MQ=255)
                                                             cccTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCtgatg                                                                                                                                         <  2:675326/88‑1 (MQ=255)
                                                             cccTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCtgatg                                                                                                                                         >  1:675326/1‑88 (MQ=255)
                                                                             aTCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACa                                                           >  2:751994/1‑150 (MQ=255)
                                                                             aTCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACa                                                           <  1:751994/150‑1 (MQ=255)
                                                                                            aGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCgctggcctgggg                                            >  1:26319/1‑148 (MQ=255)
                                                                                                  ccTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAgg                                                                            <  2:137753/112‑1 (MQ=255)
                                                                                                  ccTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAgg                                                                            >  1:137753/1‑112 (MQ=255)
                                                                                                               cACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTgg                         >  1:193335/1‑150 (MQ=255)
                                                                                                                            gCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTggg            >  2:205262/1‑150 (MQ=255)
                                                                                                                              tcatcAGCGCCACGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTGGgcg          >  1:789314/5‑150 (MQ=255)
                                                                                                                              tCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTgg                         >  1:535713/1‑135 (MQ=255)
                                                                                                                              tCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTgg                         <  2:535713/135‑1 (MQ=255)
                                                                                                                              tCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTGGgcg          >  1:238790/1‑150 (MQ=255)
                                                                                                                                 tcAGCGCCACGGCGCTGATGATGGCCATGTGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTc                                            <  1:56006/112‑1 (MQ=255)
                                                                                                                                 tcAGCGCCACGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTc                                            >  2:56006/2‑113 (MQ=255)
                                                                                                                                 tcAGCGCCACGGCGCTGATGATGGCCATGGACATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTAt                    <  1:679706/136‑1 (MQ=255)
                                                                                                                                 tcAGCGCCACGGCGCTGATGATGGCCATGGACATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTAt                    >  2:679706/2‑137 (MQ=255)
                                                                                                                                      gcCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTGGGCGCGGTGGAc  <  2:268710/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
GGCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTGGGCGCGGTGGAC  >  NC_002947/5141610‑5141893

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