Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,723,366 T→C A60A (GCT→GCC)  U712_19195 → putative transporter ywcJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,723,3660TC97.7% 157.4 / NA 44A60A (GCT→GCC) U712_19195putative transporter ywcJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (10/33);  minor base A (1/0);  total (11/33)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAACTAAACATATTTGTTTCTAGCAAAATCCGTTATGTGCTCAGATCCATTTTGGCCAGCATGTTTATCGGATTCGGCATTACAGCCGCAAGCAAAACAGGCAGTTATTTCTCTATGGCGGATTCTCCGTTTGCTTTTCCGGCAGCCGCTGTCACATTCGGGGCCGCTATTCTGCTGATTGCTTACGGAGGCGGAGATTTATTTACTGGCAACACTTTTTATTTCACTTATACCGCGCTCCGAAAAAAAATTCGCTGGCGAGACA  >  CP006881/3723230‑3723496
                                                                                                                                        |                                                                                                                                  
aaaaacaaaacaTATTCGCTTCAAGCAAAATCCGTTATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCctt                                                                                                                                  >  5:161842/3‑139 (MQ=255)
           atatTCGCTTCAAGCAAAATCCGTTATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCg                                                                                                                       <  2:59961/139‑1 (MQ=255)
                cgcttCAAGCAAAATCCGTTATGTGCTTCGGGCCATTTTGGCCAGTATATATATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGt                                                                                                                   <  7:15550/135‑1 (MQ=255)
                  cttCAAGCAAAATCCGTTATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTcac                                                                                                                >  3:226247/2‑139 (MQ=255)
                        gCAAAATCCGTTATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCgg                                                                                                          <  1:126195/139‑1 (MQ=255)
                           aaaTCCGTTATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGc                                                                                                       <  2:166619/139‑1 (MQ=255)
                             aTCCGTTATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCg                                                                                                     <  8:11211/139‑1 (MQ=255)
                             aTCCGTTATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCg                                                                                                     <  1:232462/139‑1 (MQ=255)
                             aTCCGTGATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGAAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCg                                                                                                     <  6:250228/139‑1 (MQ=255)
                                cGTTATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCt                                                                                                   <  2:137607/138‑1 (MQ=255)
                                  ttATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCATCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTAtt                                                                                                <  3:1942/139‑1 (MQ=255)
                                  ttATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGc                                                                                                       >  8:63995/1‑132 (MQ=255)
                                  ttATGTGCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGc                                                                                                       <  7:63995/132‑1 (MQ=255)
                                     tgtgCTTCGGTCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTct                                                                                              <  4:83870/138‑1 (MQ=255)
                                         cttcggtCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCtgatga                                                                                         >  4:139987/7‑139 (MQ=255)
                                            cggtCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTg                                                                                      >  3:287509/4‑139 (MQ=255)
                                            cggtCCATTTTGGCCAGTAAATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGTCAGCTATTTCTTTATGGCTGATACTCCGTTTCCCTTTACGGCAGCCGCTGTCACTTTCGGGGCCGCTATACTGATGAtt                                                                                       >  7:60297/4‑138 (MQ=255)
                                              gtCCATTTTGGCCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCt                                                                                    <  5:257016/138‑1 (MQ=255)
                                                  aTTTTGGGCAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACg                                                                                <  1:139924/139‑1 (MQ=255)
                                                          cAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGga                                                                        <  6:283885/139‑1 (MQ=255)
                                                          cAGTATATTTATTGGTTTTGGCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGga                                                                        <  5:70842/139‑1 (MQ=255)
                                                                            tggCATTACAGCCGCAAGCAAAAGAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCaa                                                       <  3:217767/137‑1 (MQ=255)
                                                                            tggCATTACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCaa                                                       <  7:77393/137‑1 (MQ=255)
                                                                                  tACAGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGACACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACAGGCAACCCcttt                                                >  7:267889/1‑139 (MQ=255)
                                                                                  tACAGACGCAAGCAAAACAGGCAGCTATTACTTTATGGCTGATTCTCCGTTTGCATTTCAGGCAGCCGCTGTCACTTTCGGGGCCGCTAATCTGATGATTGCTTACGGAGGAGTAGATTTATTTACCGGCAACACcttt                                                >  2:96978/1‑139 (MQ=255)
                                                                                     aGCCGGAAGGAAAAGAGGCAGCTAGTTCATTCTGGCTGAGTCTCCGTTTGCCTTTCCGCCAGCCTCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACTGAGGCGGAGATTTATTTACCGGCAACACCTTTTAt                                             <  5:18486/139‑1 (MQ=255)
                                                                                     aGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTAt                                             <  8:158837/139‑1 (MQ=255)
                                                                                     aGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTAt                                             <  2:11649/139‑1 (MQ=255)
                                                                                     aGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTAt                                             <  3:4344/139‑1 (MQ=255)
                                                                                     aGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGGTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTAt                                             <  5:178959/139‑1 (MQ=255)
                                                                                     aGCCGCAAGCAAAACAGGCAGCTATTTCTTTATGGCTGATGCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTAt                                             <  5:292965/139‑1 (MQ=255)
                                                                                          caagcaaAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCAc                                        <  1:46159/139‑1 (MQ=255)
                                                                                            agcaaAACAGGCAGCTATTTCTTTATGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCAcc                                       >  1:81769/1‑137 (MQ=255)
                                                                                                                    aTGGCTGATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGAATGCATACGGAGGCTGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAAta               >  8:282270/1‑137 (MQ=255)
                                                                                                                       gctgATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTctgatg                                                                                          <  1:112791/57‑1 (MQ=255)
                                                                                                                       gctgATTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTctgatg                                                                                          >  2:112791/4‑60 (MQ=255)
                                                                                                                           aTTCTCCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCCCTCCGGAAAAAAATCACCTGGCg       >  1:199121/1‑139 (MQ=255)
                                                                                                                              ctcCGTTTGCCTTTCCGTCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcg     <  4:206930/138‑3 (MQ=255)
                                                                                                                              ctcCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTTTTTCACCTATACCGCGCTGCGGAAAAAAATCAGCTGTCgcg     <  2:4707/138‑3 (MQ=255)
                                                                                                                              ctcCGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcg     <  6:111464/138‑3 (MQ=255)
                                                                                                                                 cgattGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGaca  <  1:277389/135‑1 (MQ=255)
                                                                                                                                 cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGaca  <  6:228274/138‑1 (MQ=255)
                                                                                                                                 cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGaca  <  4:226247/138‑1 (MQ=255)
                                                                                                                                 cGTTTGCCTTTCCGGCAGCCGCTGTCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCGCGaca  <  4:2627/138‑1 (MQ=255)
                                                                                                                                 cGTTTGCCTTTCCGGCAGCCGCTGGCACTTTCGGGGCCGCTATTCTGATGATTGCTTACGGAGGCGGAGATTTATTTACCGGCAACACCTTTTATTTCACCTATACCGCGCTCCGGAAAAAAATCAGCTGGCgcga    <  4:83653/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
ACAAACTAAACATATTTGTTTCTAGCAAAATCCGTTATGTGCTCAGATCCATTTTGGCCAGCATGTTTATCGGATTCGGCATTACAGCCGCAAGCAAAACAGGCAGTTATTTCTCTATGGCGGATTCTCCGTTTGCTTTTCCGGCAGCCGCTGTCACATTCGGGGCCGCTATTCTGCTGATTGCTTACGGAGGCGGAGATTTATTTACTGGCAACACTTTTTATTTCACTTATACCGCGCTCCGAAAAAAAATTCGCTGGCGAGACA  >  CP006881/3723230‑3723496

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