Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,723,952 A→C K256Q (AAA→CAA)  U712_19195 → putative transporter ywcJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,723,9520AC96.2% 182.2 / NA 53K256Q (AAA→CAA) U712_19195putative transporter ywcJ
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (35/16);  total (35/18)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCACCCTGATACAGTGACAGTAATGGGAGCAGTCAGAAACTTAATTCCCGTTACACTCGGCAACCTGACCGCGGGAATGGTCATGATGGGCTGGATGTATTACACGCTCAATCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAA  >  CP006881/3723833‑3724086
                                                                                                                       |                                                                                                                                      
aGCACCCTGATACAGTGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCa                                                                                                                        >  5:46520/1‑136 (MQ=255)
         aTACAGTGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTccc                                                                                                            >  6:262133/1‑139 (MQ=255)
               tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAg                                                                                                      >  7:58899/1‑139 (MQ=255)
               tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAg                                                                                                      >  5:60294/1‑139 (MQ=255)
               tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAAGAAAAAAACTTCCAGAGATCACTCCCTGCAAg                                                                                                      >  6:228981/1‑139 (MQ=255)
                 acactgaTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGt                                                                                                     >  5:283832/7‑138 (MQ=255)
                 acactgaTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGt                                                                                                     >  5:128886/7‑138 (MQ=255)
                   actgaTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCa                                                                                                   >  5:146866/5‑138 (MQ=255)
                       aTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCg                                                                                              <  6:81479/139‑1 (MQ=255)
                          gtaGCAGTAATAAACTTAATCCCCGTTACGCTCGGCAATATGACCGCGGTAATAGTTATGATGGGCTGGATGTACTACAAACTGAATCATGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGt                                                                                            <  8:219277/136‑1 (MQ=255)
                            aGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTAt                                                                                         >  6:288317/1‑139 (MQ=255)
                                 tCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTcc                                                                                      >  1:226922/1‑137 (MQ=255)
                                    gAAACTAAATCCCCGTAACGCGCGGCAATCTGCCCGCGGCAATAGTTATGATGGGCGGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCgagag                                                                                 <  3:184931/139‑1 (MQ=255)
                                     aaaCTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCga                                                                                    >  6:152988/1‑135 (MQ=255)
                                     aaaCTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCga                                                                                    >  5:74011/1‑135 (MQ=255)
                                            atccccGTTACGCTCGGCAATCTGACCGCGGTAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGa                                                                          <  6:276204/135‑1 (MQ=255)
                                               cccGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGAc                                                                      <  5:288317/139‑1 (MQ=255)
                                                 cGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGc                                                                    <  6:61855/139‑1 (MQ=255)
                                                     acgcgcGGAAAGCTGCCCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTgg                                                                 <  7:127944/133‑1 (MQ=255)
                                                       gcTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACa                                                              <  8:213303/138‑1 (MQ=255)
                                                            gCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAg                                                         <  2:8214/139‑1 (MQ=255)
                                                                       gcgGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGc                                              <  5:262133/139‑1 (MQ=255)
                                                                       gcgGGAATAGATAGGAGGGGCGGGAAGTACTACACACTGAATCCTGATAAATAAAAACACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGc                                              <  2:125359/139‑1 (MQ=255)
                                                                                 ttatgatgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATg                                                            <  4:88573/113‑1 (MQ=255)
                                                                                 ttatgatgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATg                                                            >  3:88573/3‑115 (MQ=255)
                                                                                   atgatgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAg                                                         >  7:16973/1‑116 (MQ=255)
                                                                                   atgatgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAg                                                         <  8:16973/116‑1 (MQ=255)
                                                                                      atgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCg                               >  4:106619/1‑139 (MQ=255)
                                                                                      atgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCg                               >  1:101943/1‑139 (MQ=255)
                                                                                        gggCGGCATGAACTACCCACTGAATCCTGATAAATAAAAACACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGc                             <  3:18509/139‑1 (MQ=255)
                                                                                         ggCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCga                                                                                    >  6:5712/1‑83 (MQ=255)
                                                                                         ggCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCga                                                                                    <  5:5712/83‑1 (MQ=255)
                                                                                         ggCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATAAAGCCGGCCTCGGATTCCGGc                             >  3:131012/1‑138 (MQ=255)
                                                                                          gCTGGATGTAATACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCaa                           <  6:155888/139‑1 (MQ=255)
                                                                                           cTGGATGTACTACACACTGAATCCTGATCAATAAAAAACCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCaaa                          >  8:5250/1‑139 (MQ=255)
                                                                                           cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATTTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCaaa                          >  1:44570/1‑139 (MQ=255)
                                                                                                   actACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAAtt                   >  1:137268/3‑138 (MQ=255)
                                                                                                     tACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGa                 >  5:62542/1‑138 (MQ=255)
                                                                                                        acactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATAtt             >  4:57446/7‑139 (MQ=255)
                                                                                                          actgaaTCCTTATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGc                             <  6:161474/117‑1 (MQ=255)
                                                                                                          actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGc                             >  5:161474/5‑121 (MQ=255)
                                                                                                          actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCa           >  2:132287/5‑139 (MQ=255)
                                                                                                          actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCa           >  4:3710/5‑139 (MQ=255)
                                                                                                            tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATg                                                            <  6:247490/86‑1 (MQ=255)
                                                                                                            tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATg                                                            >  5:247490/3‑88 (MQ=255)
                                                                                                            tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCatt         >  4:172588/3‑137 (MQ=255)
                                                                                                            tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCatt         >  1:286867/3‑137 (MQ=255)
                                                                                                            tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCatt         >  6:17971/3‑137 (MQ=255)
                                                                                                              aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCAtttt       >  7:196252/1‑139 (MQ=255)
                                                                                                              aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCAtttt       >  5:143240/1‑139 (MQ=255)
                                                                                                              aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCAtttt       >  8:159667/1‑139 (MQ=255)
                                                                                                              aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCAtttt       >  1:34827/1‑139 (MQ=255)
                                                                                                                    gATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaa  >  8:190186/1‑138 (MQ=255)
                                                                                                                       |                                                                                                                                      
AGCACCCTGATACAGTGACAGTAATGGGAGCAGTCAGAAACTTAATTCCCGTTACACTCGGCAACCTGACCGCGGGAATGGTCATGATGGGCTGGATGTATTACACGCTCAATCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAA  >  CP006881/3723833‑3724086

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: ‑