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,7430AC97.2% 132.3 / ‑3.6 36Y174D (TAC→GAC) proSproline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (16/19);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

TAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAA  >  NZ_CP009273/214659‑214824
                                                                                    |                                                                                 
tAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGttt                                                                              <  1:174049/90‑1 (MQ=255)
      cGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCa                                                                        <  1:504910/90‑1 (MQ=255)
                  gTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgag                                                            >  1:502635/1‑90 (MQ=255)
                   tACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgaga                                                           <  2:367447/90‑1 (MQ=255)
                    aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                          <  1:677794/90‑1 (MQ=255)
                    aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                          <  1:646816/90‑1 (MQ=255)
                    aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                          <  1:44262/90‑1 (MQ=255)
                       gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                       >  1:754092/1‑90 (MQ=255)
                       gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                       >  1:244565/1‑90 (MQ=255)
                       gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                       >  1:549619/1‑90 (MQ=255)
                       gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                       >  1:540020/1‑90 (MQ=255)
                       gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                       >  1:507115/1‑90 (MQ=255)
                            gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa                                                  <  2:781641/90‑1 (MQ=255)
                            gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa                                                  <  1:125544/90‑1 (MQ=255)
                            gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa                                                  <  2:549619/90‑1 (MQ=255)
                            gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa                                                  <  1:556935/90‑1 (MQ=255)
                                  cAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTa                                            <  1:705270/90‑1 (MQ=255)
                                     ccccATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGc                                         <  1:406390/90‑1 (MQ=255)
                                         aTGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCt                                     >  1:63388/1‑90 (MQ=255)
                                         aTGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCt                                     >  2:778967/1‑90 (MQ=255)
                                                 gAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAg                             >  1:188200/1‑90 (MQ=255)
                                                           cTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAATCATCTTTCATCAGGAATTCGCgg                   <  1:1270/90‑1 (MQ=255)
                                                            tGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCggg                  <  2:754092/90‑1 (MQ=255)
                                                             gTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGa                 >  1:476668/1‑90 (MQ=255)
                                                              tAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgaga                                                           <  1:125665/47‑1 (MQ=255)
                                                              tAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgaga                                                           >  2:125665/1‑47 (MQ=255)
                                                                  ccGCATACATTGCATCGTCGGTTTCCTGCAGGGAt                                                                   <  2:341857/35‑1 (MQ=38)
                                                                  ccGCATACATTGCATCGTCGGTTTCCTGCAGGGAt                                                                   >  1:341857/1‑35 (MQ=38)
                                                                    gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATg          >  1:207477/1‑90 (MQ=255)
                                                                    gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATg          >  1:307317/1‑90 (MQ=255)
                                                                    gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATg          >  1:339745/1‑90 (MQ=255)
                                                                           ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                          >  1:753554/1‑35 (MQ=38)
                                                                           ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                          <  2:753554/35‑1 (MQ=255)
                                                                           ttGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa   >  1:451883/1‑90 (MQ=255)
                                                                            tGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGaa  <  2:502635/90‑1 (MQ=255)
                                                                            tGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGaa  <  1:778967/90‑1 (MQ=255)
                                                                                    |                                                                                 
TAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAA  >  NZ_CP009273/214659‑214824

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: