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,9520AC98.2% 198.1 / NA 54K256Q (AAA→CAA) U712_19195putative transporter ywcJ
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (44/9);  total (44/10)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACAGTAATGGGAGCAGTCAGAAACTTAATTCCCGTTACACTCGGCAACCTGACCGCGGGAATGGTCATGATGGGCTGGATGTATTACACGCTCAATCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAA  >  CP006881/3723848‑3724086
                                                                                                        |                                                                                                                                      
tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTcc                                                                                                             >  8:325434/1‑132 (MQ=255)
tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTcc                                                                                                             <  7:325434/132‑1 (MQ=255)
tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAg                                                                                                      >  6:245221/1‑139 (MQ=255)
tGACACTGATGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCACTGCAAg                                                                                                      >  5:217126/1‑139 (MQ=255)
    actgttGGGCGCAGTCATAAAGTTAATCCCCGTTACGCTCGGCAATCTCACCGCGGGAATAGTTATCAGGGGCTGGATGTACAACACACTGAATCCTGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCaa                                                                                                  <  4:28646/134‑1 (MQ=255)
    actgaTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCa                                                                                                   >  2:326880/5‑138 (MQ=255)
        aTGGGAGCAGTCAGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGATTAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCg                                                                                              >  3:209976/1‑139 (MQ=255)
                    aGAAACTTAATCCCCGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATACAAACACTTCCAGAGATCACTCCCTGCCAGGCAAAGCGGTTATTCCga                                                                                    >  5:216478/1‑137 (MQ=255)
                             agacccGTTACGCTCGGAATACTTCCCGCGGTAATATTTAGTATGGGCTTCATTAAATCCACACTGAATCCTCAGCAAAAAAAAAACTTCCACAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGaa                                                                         <  6:217126/136‑1 (MQ=255)
                                  cGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGa                                                                          <  6:89882/133‑1 (MQ=255)
                                  cGTTACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGa                                                                          >  5:89882/1‑133 (MQ=255)
                                    ttACGCTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTg                                                                  >  1:149901/1‑139 (MQ=255)
                                    ttACGCTCGGCAATCTGACAGAGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCATTATCAATAAAAAAACTTCCAGAGATAACTCCCTGCAAGTCAAAGCGGTTATTCCTACAGATCTTGAAGACGCTg                                                                  >  5:264352/1‑139 (MQ=255)
                                        gcTCGCCAATCTGACCCGGGGAATCGTTATTATGGGCTGGATGAACTACACACTGAATCCTGATCAAAAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACa                                                              <  1:330610/138‑1 (MQ=255)
                                         cTCGGCAATCTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACAt                                                             >  5:296547/1‑139 (MQ=255)
                                                atcTGACCGCGGGAATAGTTATGATGGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGt                                                      >  1:290338/3‑139 (MQ=255)
                                                                    atgatgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAg                                                         <  7:103376/116‑1 (MQ=255)
                                                                    atgatgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAg                                                         >  8:103376/1‑116 (MQ=255)
                                                                    atgatgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGatt                                  >  8:239183/1‑137 (MQ=255)
                                                                       atgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCt                                        <  5:305640/130‑1 (MQ=255)
                                                                       atgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCt                                        >  6:305640/1‑130 (MQ=255)
                                                                       atgGGCTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCg                               >  6:188965/1‑139 (MQ=255)
                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGCTCACTCCCTGCAAGTCACAGCGGTTATTCAGCGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCaaa                          >  8:143016/1‑139 (MQ=255)
                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGt                                                                                            <  4:96276/73‑1 (MQ=255)
                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGt                                                                                            >  3:96276/1‑73 (MQ=255)
                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCaaa                          >  5:305339/1‑139 (MQ=255)
                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCaaa                          >  4:244464/1‑139 (MQ=255)
                                                                            cTGGATGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCaaa                          >  2:86475/1‑139 (MQ=255)
                                                                                 tGTACTACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTaa                     >  2:134784/1‑139 (MQ=255)
                                                                                     ctACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGa                 >  3:19404/2‑139 (MQ=255)
                                                                                      tACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGa                 >  7:150292/1‑138 (MQ=255)
                                                                                      tACACACTGAATCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGa                 >  6:325907/1‑138 (MQ=255)
                                                                                         acactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGc                             <  2:13085/117‑1 (MQ=255)
                                                                                         acactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGc                             >  1:13085/7‑123 (MQ=255)
                                                                                         acactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGatat              >  5:125617/7‑138 (MQ=255)
                                                                                         acactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGatat              >  5:99576/7‑138 (MQ=255)
                                                                                         acactgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATAtt             >  6:277170/7‑139 (MQ=255)
                                                                                           actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCa           >  8:46539/5‑139 (MQ=255)
                                                                                           actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCa           >  7:145783/5‑139 (MQ=255)
                                                                                           actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCa           >  8:290634/5‑139 (MQ=255)
                                                                                           actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCa           >  3:331323/5‑139 (MQ=255)
                                                                                           actgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCa           >  3:262904/5‑139 (MQ=255)
                                                                                             tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCatt         >  3:75098/3‑137 (MQ=255)
                                                                                             tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCatt         >  2:197386/3‑137 (MQ=255)
                                                                                             tgaaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAAAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCatt         >  4:206753/3‑137 (MQ=255)
                                                                                               aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGt                                                                                            >  7:280122/1‑54 (MQ=255)
                                                                                               aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGt                                                                                            <  8:280122/54‑1 (MQ=255)
                                                                                               aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCAtttt       >  3:107229/1‑139 (MQ=255)
                                                                                               aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCAtttt       >  1:298401/1‑139 (MQ=255)
                                                                                               aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATGCCGGCAAATGTAATTGATATTCAtttt       >  6:183993/1‑139 (MQ=255)
                                                                                                    tgATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaa   >  2:13565/2‑138 (MQ=255)
                                                                                                     gATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaa  >  6:303449/1‑138 (MQ=255)
                                                                                                     gATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaa  >  6:268443/1‑138 (MQ=255)
                                                                                                     gATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaa  >  5:301768/1‑138 (MQ=255)
                                                                                                        |                                                                                                                                      
TGACAGTAATGGGAGCAGTCAGAAACTTAATTCCCGTTACACTCGGCAACCTGACCGCGGGAATGGTCATGATGGGCTGGATGTATTACACGCTCAATCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAA  >  CP006881/3723848‑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: