Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,723,948 C→T P254P (CCC→CCT U712_19195 → putative transporter ywcJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,723,9480CT98.3% 199.4 / NA 58P254P (CCC→CCTU712_19195putative transporter ywcJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (39/18);  minor base G (0/1);  total (39/19)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCACCCTGATACAGTGACAGTAATGGGAGCAGTCAGAAACTTAATTCCCGTTACACTCGGCAACCTGACCGCGGGAATGGTCATGATGGGCTGGATGTATTACACGCTCAATCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTT  >  CP006881/3723832‑3724081
                                                                                                                    |                                                                                                                                     
gAGCACCCTGATACAGTGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCa                                                                                                                   <  2:225915/137‑1 (MQ=255)
       attaTACAGTTAAACTGATGGGAGAAGACTTAAAGTTAATCCCGGTTACGCTCTGCAATCTGACCGCGGGAATATTTATGTTGGGCTGGAGGTACTAAACAGTTAATCCTGTTCAATAAAAAAACTTCCAGAGATCACt                                                                                                          <  6:230799/136‑1 (MQ=255)
          aTACAGTGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTccc                                                                                                       >  4:284232/1‑139 (MQ=255)
            aCAGTGACACTGATGGTAGCAGTCAGAAAATTAATCCCGGTTACGCTGGGCAAGCTGACGGCGGTAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTTATCAATAAAAAAACTTCCAGAGATCACTCCCTg                                                                                                     <  6:21423/139‑1 (MQ=255)
                tGACACTTATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGCCAATCTGCCCGCGGGAAGAGTTTGTATGGTCTGGATGTACTACACCCTGAATCCTGATCAATAAAAAACCTTCCAGAGATCACTCCCTGAAAg                                                                                                 >  6:44598/1‑139 (MQ=255)
                tGACACTGATGGGAGCAGTCAGAAACTTAATCCGCGTTACGCTCGGCAATCTGACAGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCTAg                                                                                                 >  8:153093/1‑139 (MQ=255)
                tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAg                                                                                                 >  1:90842/1‑139 (MQ=255)
                tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAg                                                                                                 >  5:284841/1‑139 (MQ=255)
                tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAg                                                                                                 >  8:17682/1‑139 (MQ=255)
                tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAg                                                                                                 >  5:21423/1‑139 (MQ=255)
                tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAg                                                                                                 >  4:159420/1‑139 (MQ=255)
                tGACACTGATGGGAGCAGTCAGAAACATAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAg                                                                                                 >  5:230799/1‑139 (MQ=255)
                  acactgaTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGt                                                                                                >  8:233630/7‑138 (MQ=255)
                    cctgaTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGCCCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTTAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCaa                                                                                             <  1:211762/135‑1 (MQ=255)
                    actgaTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCa                                                                                                            <  3:89479/120‑1 (MQ=255)
                    actgaTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCa                                                                                                            >  4:89479/5‑124 (MQ=255)
                    actgaTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCaa                                                                                             <  3:159420/135‑1 (MQ=255)
                    actgaTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCa                                                                                              >  2:245913/5‑138 (MQ=255)
                           ggAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGt                                                                                       <  1:245913/138‑1 (MQ=255)
                                 gTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCg                                                                                <  6:284841/139‑1 (MQ=255)
                                    aGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCga                                                                               >  5:295551/1‑137 (MQ=255)
                                            aaTCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCa                                                                                                                   <  6:117957/93‑1 (MQ=255)
                                            aaTCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTg                                                                      >  7:192262/1‑138 (MQ=255)
                                            aaTCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTg                                                                      >  1:115879/1‑138 (MQ=255)
                                             atccccGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGa                                                                     <  7:17682/135‑1 (MQ=255)
                                                  cGTTACGCTCGGCAATCTGACAGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGc                                                               >  5:89295/1‑139 (MQ=255)
                                                    ttACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTg                                                             >  3:249740/1‑139 (MQ=255)
                                                        gcTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACa                                                         <  7:233630/138‑1 (MQ=255)
                                                            ggCAATCTGACCGCGGGAATAGTTAGGATGGGCTGGATGTACTACACACTGAATCCTGAGCAAAAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATg                                                       <  7:52939/137‑1 (MQ=255)
                                                             gCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAg                                                    <  7:212227/139‑1 (MQ=255)
                                                             gCAAGCTGCCCGCGGTAATAGTTATGATGGGCTGGATGAACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAg                                                    <  3:99462/139‑1 (MQ=255)
                                                                        gcgGTAATAGTTAGGAGGGGCTGGATGTACTACACACTGAATCCGGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGc                                         <  4:249740/139‑1 (MQ=255)
                                                                                  ttatgatgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCCCTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGa                               >  7:20743/3‑139 (MQ=255)
                                                                                          ggCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGc                        >  8:204072/1‑138 (MQ=255)
                                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAg                                                    <  3:32747/108‑1 (MQ=255)
                                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAg                                                    >  4:32747/1‑108 (MQ=255)
                                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGc                        <  5:30472/136‑1 (MQ=255)
                                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGc                        >  6:30472/1‑136 (MQ=255)
                                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCaaa                     >  2:48449/1‑139 (MQ=255)
                                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCa                       <  7:197525/137‑1 (MQ=255)
                                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCa                       >  8:197525/1‑137 (MQ=255)
                                                                                                    actACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAAtt              >  8:48099/3‑138 (MQ=255)
                                                                                                         acactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTaa                <  4:33709/125‑1 (MQ=255)
                                                                                                         acactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTaa                >  3:33709/7‑131 (MQ=255)
                                                                                                          cactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATAtt        >  7:95835/6‑138 (MQ=255)
                                                                                                           actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCa      >  8:68661/5‑139 (MQ=255)
                                                                                                           actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCa      >  2:2373/5‑139 (MQ=255)
                                                                                                           actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCa      >  3:37568/5‑139 (MQ=255)
                                                                                                           actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATCTAATTGATATTCa      >  4:22936/5‑139 (MQ=255)
                                                                                                             tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCatt    >  7:271466/3‑137 (MQ=255)
                                                                                                             tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCatt    >  3:14535/3‑137 (MQ=255)
                                                                                                             tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCatt    >  5:336897/3‑137 (MQ=255)
                                                                                                             tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCatt    >  5:8000/3‑137 (MQ=255)
                                                                                                             tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCACTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATCTCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCatt    >  8:316617/3‑137 (MQ=255)
                                                                                                               aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCAtttt  >  2:4102/1‑139 (MQ=255)
                                                                                                               aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCAtttt  >  2:19962/1‑139 (MQ=255)
                                                                                                               aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCAtttt  >  7:137462/1‑139 (MQ=255)
                                                                                                               aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCGGGACCTGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAAAGTAATTCATATTCAtttt  >  7:202971/1‑139 (MQ=255)
                                                                                                                    |                                                                                                                                     
GAGCACCCTGATACAGTGACAGTAATGGGAGCAGTCAGAAACTTAATTCCCGTTACACTCGGCAACCTGACCGCGGGAATGGTCATGATGGGCTGGATGTATTACACGCTCAATCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTT  >  CP006881/3723832‑3724081

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