Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,832,108 C→T 22.8% T124T (ACG→ACA astA ← arginine succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,832,1080CT22.8% 116.5 / 30.6 57T124T (ACG→ACAastAarginine succinyltransferase
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base T (6/7);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AACATAAAGCGCGATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATA  >  W3110S.gb/1832040‑1832170
                                                                    |                                                              
aacataaaGCGCGATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGc                                                                >  1:1268490/1‑69 (MQ=255)
aacataaaGCGCGATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGc                                                                >  1:2020592/1‑69 (MQ=255)
aacataaaGCGCGATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGc                                                                >  1:3607805/1‑69 (MQ=255)
aacataaaGCGCGATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGc                                                                >  1:3276880/1‑69 (MQ=255)
  cataaaGCGCGATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCTGGTCGAGAAACAGCg                                                               <  1:1552669/68‑1 (MQ=255)
   ataaaGCGCGATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTg                                                             <  1:1493127/69‑1 (MQ=255)
    taaaGCGCGATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTg                                                             >  1:1253968/1‑68 (MQ=255)
      aaGCGCGATTTCGACAGCAAATAGCCGTTGCCGTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGcac                                                          >  1:257386/1‑69 (MQ=255)
      aaGCGCGATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGcac                                                          >  1:203256/1‑69 (MQ=255)
      aaGCGCGATTTCGACAGAAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGcac                                                          >  1:300551/1‑69 (MQ=255)
         cgcgATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGc                                                       >  1:2606336/1‑69 (MQ=255)
           cgATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTc                                                     <  1:2771880/69‑1 (MQ=255)
              tttCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCgct                                                  <  1:1709177/69‑1 (MQ=255)
              tttCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCgct                                                  <  1:2680550/69‑1 (MQ=255)
              tttCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCgct                                                  <  1:695342/69‑1 (MQ=255)
              tttCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCgct                                                  <  1:23750/69‑1 (MQ=255)
              tttCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCgct                                                  <  1:2123913/69‑1 (MQ=255)
              tttCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCgct                                                  <  1:437877/69‑1 (MQ=255)
              tttCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCgct                                                  <  1:604471/69‑1 (MQ=255)
                      gCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCgg                                          <  1:69978/69‑1 (MQ=255)
                        aaaTAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGgt                                         >  1:1027948/1‑68 (MQ=255)
                        aaaTAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGgt                                         >  1:2267079/1‑68 (MQ=255)
                          aTAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGtgtg                                      <  1:1642711/69‑1 (MQ=255)
                           tAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGa                                     >  1:267551/1‑69 (MQ=255)
                            aGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTG‑‑CAGCTCGCTGCTGCCGGTGTGATc                                   >  1:1348077/1‑68 (MQ=255)
                             gCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATc                                   <  1:1838167/69‑1 (MQ=255)
                               cGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGt                                 >  1:3337943/1‑69 (MQ=255)
                               cGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGt                                 >  1:2129443/1‑69 (MQ=255)
                               cGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGt                                 >  1:1734891/1‑69 (MQ=255)
                               cGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGt                                 <  1:1803761/69‑1 (MQ=255)
                               cGTTGCCCTCTTTGCGCCAGGCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGt                                 <  1:2949082/69‑1 (MQ=255)
                                 ttGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGtt                                <  1:2076934/68‑1 (MQ=255)
                                 ttGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTa                               <  1:1813388/69‑1 (MQ=255)
                                 ttGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTa                               <  1:502509/69‑1 (MQ=255)
                                   gCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTa                               >  1:2324912/1‑67 (MQ=255)
                                   gCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTa                               >  1:1915411/1‑67 (MQ=255)
                                       tctTTGCGCCAGTCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGt                                 >  1:3722912/1‑61 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTgaga                         >  1:22571/1‑67 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTgaga                         >  1:2227460/1‑67 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTgaga                         >  1:2174575/1‑67 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTgaga                         >  1:1548267/1‑67 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTgaga                         >  1:383783/1‑67 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCg                                                    >  1:2122994/1‑40 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCg                                                    >  1:1561664/1‑40 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCg                                                    >  1:1762965/1‑40 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGtgtg                                      >  1:1355018/1‑54 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGtgtg                                      >  1:3686235/1‑54 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGtgtg                                      >  1:1932867/1‑54 (MQ=255)
                                         tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGtgtg                                      >  1:661467/1‑54 (MQ=255)
                                                  gTCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCgg              <  1:3317668/69‑1 (MQ=255)
                                                   tCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGt                                 <  1:418896/49‑1 (MQ=255)
                                                   tCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGt                                 <  1:3216094/49‑1 (MQ=255)
                                                   tCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGt                                 <  1:3543445/49‑1 (MQ=255)
                                                   tCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGt                                 <  1:347009/49‑1 (MQ=255)
                                                   tCCGGGTCGAGAAACAGTGTGCACAGCTCGCTGCTGCCGGTGTGATCGt                                 <  1:1405942/49‑1 (MQ=255)
                                                   tCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGc             <  1:442895/69‑1 (MQ=255)
                                                   tCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGGTACTGAGAAACAGCGTCGGc             <  1:1041973/69‑1 (MQ=255)
                                                     cGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGgtgt                                       >  1:3231212/1‑41 (MQ=255)
                                                             gaAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTta    <  1:385687/68‑1 (MQ=255)
                                                             gaAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTta    <  1:831856/68‑1 (MQ=255)
                                                              aaaCAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTtata  <  1:3317996/69‑1 (MQ=255)
                                                                    |                                                              
AACATAAAGCGCGATTTCGACAGCAAATAGCCGTTGCCCTCTTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATA  >  W3110S.gb/1832040‑1832170

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