Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 214,743 A→C Y174D (TAC→GAC)  proS ← proline‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273214,7430AC95.0% 70.4 / ‑3.1 20Y174D (TAC→GAC) proSproline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (10/9);  total (11/9)
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.95e-01

ACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGC  >  NZ_CP009273/214679‑214828
                                                                |                                                                                     
aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                              <  2:502037/90‑1 (MQ=255)
  ggCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGt                                                            <  1:589744/90‑1 (MQ=255)
   gcgcGGAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGAt                                                                       >  1:386223/1‑77 (MQ=255)
   gcgcGGAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGAt                                                                       <  2:386223/77‑1 (MQ=255)
        gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa                                                      <  2:465032/90‑1 (MQ=255)
        gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa                                                      <  1:245916/90‑1 (MQ=255)
        gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa                                                      >  1:421106/1‑90 (MQ=255)
               aGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                              <  1:255417/75‑1 (MQ=255)
               aGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                              >  2:255417/1‑75 (MQ=255)
               aGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTaa                                               >  1:571831/1‑90 (MQ=255)
                                    ttGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTcg                          >  2:324727/1‑90 (MQ=255)
                                      gCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTcgcg                        <  1:579657/90‑1 (MQ=255)
                                                gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATg              >  1:4759/1‑90 (MQ=255)
                                                gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTAGCGGGAACGCATg              >  1:623461/1‑90 (MQ=255)
                                                       ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGaa                    >  1:1731/1‑77 (MQ=255)
                                                       ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGaa                    <  2:1731/77‑1 (MQ=255)
                                                       ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa       >  1:464507/1‑90 (MQ=255)
                                                           aTCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa                                                      <  2:604844/39‑1 (MQ=39)
                                                           aTCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa                                                      >  1:604844/1‑39 (MQ=39)
                                                           aTCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAcg   >  2:497515/1‑90 (MQ=255)
                                                            tcgtcgGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAcgc  <  2:594708/90‑1 (MQ=255)
                                                                |                                                                                     
ACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGC  >  NZ_CP009273/214679‑214828

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: