Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,363,696 T→C F16S (TTC→TCC)  pspD → phage shock protein PspD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,363,6960TC94.6% 106.2 / ‑0.4 37F16S (TTC→TCC) pspDphage shock protein PspD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (16/19);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

TCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCT  >  NZ_CP009273/1363610‑1363781
                                                                                      |                                                                                     
tCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCaa                                                                                    >  2:579441/1‑90 (MQ=255)
      tACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGc                                                                              >  1:479067/1‑90 (MQ=255)
      tACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGc                                                                              >  1:98559/1‑90 (MQ=255)
       aCGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCa                                                                             <  1:1056816/90‑1 (MQ=255)
            gCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCaa                                                                        <  1:187372/90‑1 (MQ=255)
             ccGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAg                                                                       <  2:479067/90‑1 (MQ=255)
               gTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCt                                                                     >  2:624565/1‑90 (MQ=255)
               gTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCt                                                                     <  2:1279194/90‑1 (MQ=255)
                 ttCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTgg                                                                   <  2:1409687/90‑1 (MQ=255)
                    cGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc                                                                <  2:754372/90‑1 (MQ=255)
                    cGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc                                                                <  1:1346780/90‑1 (MQ=255)
                    cGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc                                                                <  1:457053/90‑1 (MQ=255)
                         aCTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActtctt                                                           >  1:121748/1‑90 (MQ=255)
                         aCTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActtctt                                                           >  1:542423/1‑90 (MQ=255)
                          cTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTa                                                          >  2:1196623/1‑90 (MQ=255)
                           tgtgAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTAc                                                         <  1:714729/90‑1 (MQ=255)
                           tgtgAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTAc                                                         <  1:579441/90‑1 (MQ=255)
                              gAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGc                                                      <  1:173297/90‑1 (MQ=255)
                                ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAc                                                    <  1:511519/90‑1 (MQ=255)
                                ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAc                                                    <  1:1284839/90‑1 (MQ=255)
                                ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCAc                                                    <  1:1431264/90‑1 (MQ=255)
                                   aaGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTgc                                                 <  1:1407143/90‑1 (MQ=255)
                                      ttATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCTAAATTAGCAGGCAAg                                                                       >  1:524180/1‑65 (MQ=255)
                                      ttATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCTAAATTAGCAGGCAAg                                                                       <  2:524180/65‑1 (MQ=255)
                                      ttATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAg                                                                       <  2:821860/65‑1 (MQ=255)
                                      ttATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAg                                                                       >  1:821860/1‑65 (MQ=255)
                                          gAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATgg                                          >  1:360459/1‑90 (MQ=255)
                                               cTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc                                                                >  1:444662/1‑63 (MQ=255)
                                               cTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTAc                                                                <  2:444662/63‑1 (MQ=255)
                                               cTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCcgg                                     >  1:640300/1‑90 (MQ=255)
                                                    tGGCAACAGGGCGGGCAAAAGGTAAAGCCTGGTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCgggg                                >  2:1187425/1‑89 (MQ=255)
                                                                gggCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTa                                             >  2:738616/1‑65 (MQ=255)
                                                                gggCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTa                                             <  1:738616/65‑1 (MQ=255)
                                                                           aaaGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActt                                                              >  2:1326228/1‑37 (MQ=39)
                                                                           aaaGCCTGGTTCCAAATTAGCAGGCAAGCTGGTActt                                                              <  1:1326228/37‑1 (MQ=39)
                                                                                  ggTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGGGGTAAAAACAGTTGCt  >  2:472303/1‑90 (MQ=255)
                                                                                  ggTTCCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGGGTGGCGGGGGGGGGGGTAAAAAAAGtttgt  >  2:1018529/1‑87 (MQ=255)
                                                                                      |                                                                                     
TCACGTTACGTAGCCGTTTCCGTCAACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCT  >  NZ_CP009273/1363610‑1363781

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