Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 707,707 A→G 25.9% N398S (AAC→AGC)  glnS → glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb707,7070AG25.9% 138.0 / 57.7 81N398S (AAC→AGC) glnSglutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base G (9/12);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

ATGGGCAGCCGTCAGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGATTAAGGCAGAA  >  W3110S.gb/707637‑707777
                                                                      |                                                                      
aTGGGCAGCCGTCAGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTaa                                                                        >  1:1845356/1‑71 (MQ=255)
   ggCAGCCGTCAGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAAc                                                                       >  1:1685040/1‑69 (MQ=255)
     cAGCCGTCAGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAGc                                                                   <  1:4374130/71‑1 (MQ=255)
          gTCAGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAGCAGTa                                                               >  1:3697402/1‑70 (MQ=255)
           tCAGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAGCAGt                                                                <  1:1020497/68‑1 (MQ=255)
           tCAGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAGCAGt                                                                <  1:3428612/68‑1 (MQ=255)
           tCAGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAGCAGTACa                                                             <  1:2627472/71‑1 (MQ=255)
           tCAGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACa                                                             <  1:2031645/71‑1 (MQ=255)
             aGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAGCAGTACaaa                                                           >  1:321596/1‑71 (MQ=255)
             aGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACaa                                                            >  1:3604632/1‑70 (MQ=255)
                 gCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACaaa                                                           >  1:3715587/1‑67 (MQ=255)
                 gCCGTTAAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTc                                                       >  1:660665/1‑71 (MQ=255)
                  ccGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCt                                                      >  1:364920/1‑71 (MQ=255)
                  ccGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCt                                                      >  1:226515/1‑71 (MQ=255)
                   cGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAGCAGTACAAACGTCTg                                                     >  1:392973/1‑71 (MQ=255)
                     tttAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAg                                                                    <  1:1620125/54‑1 (MQ=255)
                     tttAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAg                                                                    <  1:3556115/54‑1 (MQ=255)
                     tttAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAg                                                                    <  1:3458585/54‑1 (MQ=255)
                     tttAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAg                                                                    <  1:3254658/54‑1 (MQ=255)
                      ttAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGt                                                   >  1:752368/1‑70 (MQ=255)
                      ttAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTg                                                  >  1:1173338/1‑71 (MQ=255)
                      ttAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTg                                                  >  1:3626650/1‑71 (MQ=255)
                      ttAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTg                                                  >  1:1929327/1‑71 (MQ=255)
                        aGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGAAGTACAAACGTCTGGTGc                                                 <  1:3877296/70‑1 (MQ=255)
                             tGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACg                                                         <  1:3372980/57‑1 (MQ=255)
                              gagaTTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTa                                           >  1:2301056/1‑70 (MQ=255)
                               agaTTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAg                                                                    <  1:1621600/44‑1 (MQ=255)
                               agaTTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAg                                                                    <  1:4743078/44‑1 (MQ=255)
                               agaTTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAg                                                                    <  1:609098/44‑1 (MQ=255)
                               agaTTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAg                                                                    <  1:2089946/44‑1 (MQ=255)
                               agaTTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAGCAGTACAAACGTCTgg                                                    >  1:2353162/1‑60 (MQ=255)
                               agaTTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACaaa                                                           >  1:4120075/1‑53 (MQ=255)
                               agaTTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACaaa                                                           >  1:1407126/1‑53 (MQ=255)
                                gaTTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaag                                        <  1:2291013/71‑1 (MQ=255)
                                 aTTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAg                                                                 <  1:4755183/45‑1 (MQ=255)
                                  tttGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAAc                                                                       <  1:4745059/38‑1 (MQ=255)
                                  tttGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAAc                                                                       <  1:3555068/38‑1 (MQ=255)
                                   ttGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAGCAGTACAAACGTCTGGTGCTg                                               >  1:4493862/1‑61 (MQ=255)
                                   ttGGATTGATCGCGCCGAGTTCCGCGAAGAAGCTAGCAAGCAGTACAAACGTCTGGTGCTGGGTAaagaag                                     >  1:2202265/1‑71 (MQ=255)
                                    tGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaagaag                                     >  1:3727299/1‑70 (MQ=255)
                                    tGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaagaag                                     >  1:3932910/1‑70 (MQ=255)
                                     ggATTGATCGCGCCGATTTCCGCGAAGAAGCTAGCAAGCAGTa                                                               >  1:3703881/1‑43 (MQ=255)
                                       aTTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaag                                        <  1:642083/64‑1 (MQ=255)
                                       aTTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaag                                        <  1:446420/64‑1 (MQ=255)
                                       aTTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaag                                        <  1:2275136/64‑1 (MQ=255)
                                        ttGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGc                                                                   <  1:2213521/36‑1 (MQ=255)
                                         tGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaagaag                                     >  1:3609350/1‑65 (MQ=255)
                                         tGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaagaa                                      >  1:4288168/1‑64 (MQ=255)
                                         tGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaagaa                                      >  1:678774/1‑64 (MQ=255)
                                         tGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaagaa                                      >  1:4293967/1‑64 (MQ=255)
                                         tGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaag                                        <  1:385620/62‑1 (MQ=255)
                                          gATCGCGCCGATTTCCGCGAAGAAGCTAGCAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCt                              >  1:2747374/1‑71 (MQ=255)
                                          gATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGt                                            <  1:140918/57‑1 (MQ=255)
                                           aTCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaagaa                                      >  1:5083050/1‑62 (MQ=255)
                                           aTCGCGCCGATTTC‑‑CGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaagaa                                      >  1:4410110/1‑60 (MQ=255)
                                             cgcgCCGATTTCCGCGAAGAAGCTAGCAAGCAGTACAAACGTCTGGTGCTGGGTAaagaa                                      >  1:1632199/1‑60 (MQ=255)
                                              gcgcCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAAc                                                          <  1:3324994/39‑1 (MQ=255)
                                              gcgcCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCg                           >  1:410661/1‑70 (MQ=255)
                                              gcgaCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAAc                                                          <  1:128059/39‑1 (MQ=255)
                                               cgcCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTa                         <  1:159725/71‑1 (MQ=255)
                                               cgcCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTa                         <  1:4953329/71‑1 (MQ=255)
                                                gcCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTa                         <  1:153747/70‑1 (MQ=255)
                                                  cGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaagaag                                     >  1:464704/1‑56 (MQ=255)
                                                  cGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATg                      >  1:5062101/1‑71 (MQ=255)
                                                   gATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGc                     >  1:2413279/1‑71 (MQ=255)
                                                   gATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGc                     >  1:1176186/1‑71 (MQ=255)
                                                   gATTTCCGCGAAGAAGATAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGc                     >  1:3267300/1‑71 (MQ=255)
                                                    aTTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAAt                       <  1:2518244/68‑1 (MQ=255)
                                                    aTTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAAt                       <  1:3842677/68‑1 (MQ=255)
                                                    aTTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAAt                       <  1:3325574/68‑1 (MQ=255)
                                                      ttCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTaaa                                         >  1:2234489/1‑48 (MQ=255)
                                                            gaagaaGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGa             >  1:125604/1‑70 (MQ=255)
                                                            gaagaaGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGAt            <  1:3004980/71‑1 (MQ=255)
                                                             aagaagCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGAtt           <  1:1848622/71‑1 (MQ=255)
                                                              agaagCTAACAAGCTGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAAt                       <  1:1902348/58‑1 (MQ=255)
                                                               gaagCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGATTaa         <  1:2847531/71‑1 (MQ=255)
                                                                 agCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGATTAAgg       <  1:5085969/71‑1 (MQ=255)
                                                                  gCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGc                                  <  1:2327840/43‑1 (MQ=255)
                                                                  gCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTAtgt               <  1:2200250/62‑1 (MQ=255)
                                                                   cTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaag                                        <  1:3430576/36‑1 (MQ=255)
                                                                   cTAACAAGCAGTACAAACGTCTGGTGCTGGGTAaag                                        <  1:943185/36‑1 (MQ=255)
                                                                   cTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGATTa          <  1:809893/66‑1 (MQ=255)
                                                                     aaCAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGATTAAg        >  1:1528763/1‑66 (MQ=255)
                                                                      aCAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGATTATGGCAGa   <  1:657994/70‑1 (MQ=255)
                                                                      aCAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGATTAAGGCAGaa  >  1:3554467/1‑71 (MQ=255)
                                                                      aCAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGATTAAGGCAGaa  >  1:3235483/1‑71 (MQ=255)
                                                                      aCAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGATTAAGGCAGa   <  1:714796/70‑1 (MQ=255)
                                                                      |                                                                      
ATGGGCAGCCGTCAGGTGCCGTTTAGCGGTGAGATTTGGATTGATCGCGCCGATTTCCGCGAAGAAGCTAACAAGCAGTACAAACGTCTGGTGCTGGGTAAAGAAGTGCGTCTGCGTAATGCTTATGTGATTAAGGCAGAA  >  W3110S.gb/707637‑707777

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: