Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,343,839 T→C T176A (ACA→GCA)  yciW ← CMD domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,343,8390TC95.5% 68.9 / ‑4.6 22T176A (ACA→GCA) yciWCMD domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (8/13);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTGACGCGTCGGCGTAGTTTTGTATCTCCAGCCCGGGTAACCAGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGTCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATCACACAGTCCCAGGCAGTAAGACCGACGT  >  NZ_CP009273/1343755‑1343910
                                                                                    |                                                                       
gTGACGCGTCGGCGTAGTTTTGTATCTCCAGCCCGGGTAACCAGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGTCcgcg                                                                    >  1:784949/1‑90 (MQ=255)
          ggCGTAGTTTTGTATCTCCAGCCCGGGTAACCAGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGcc                                                          >  1:464509/1‑90 (MQ=255)
              tAGTTTTGTATCTCCAGCCCGGGTAACCAGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATa                                                      <  1:511717/90‑1 (MQ=255)
              tAGTTTTGTATCTCCAGCCCGGGTAACCAGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATa                                                      <  2:1054902/90‑1 (MQ=255)
                         ctcCAGCCCGGGTAACCAGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGAt                                           <  2:464509/90‑1 (MQ=255)
                          tcCAGCCCGGGTAACCAGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGAtt                                          >  1:194064/1‑90 (MQ=255)
                                        ccAGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATcac                            <  2:272634/90‑1 (MQ=255)
                                          aGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATcacac                          <  1:871209/90‑1 (MQ=255)
                                          aGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATcacac                          <  2:386554/90‑1 (MQ=255)
                                          aGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATcacac                          <  1:267507/90‑1 (MQ=255)
                                          aGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATcacac                          <  1:267506/90‑1 (MQ=255)
                                           gCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTa                                    >  1:590389/1‑79 (MQ=255)
                                           gCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTa                                    <  2:590389/79‑1 (MQ=255)
                                             gTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATCACACAGt                       >  1:197431/1‑90 (MQ=255)
                                                  ggATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCTCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATCACACAGTCCCAg                  >  1:152430/1‑90 (MQ=255)
                                                  ggATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGaaa                                                             <  1:753743/47‑1 (MQ=255)
                                                  ggATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGaaa                                                             <  1:565906/47‑1 (MQ=255)
                                                  ggATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGaaa                                                             >  2:565906/1‑47 (MQ=255)
                                                  ggATGCCCGAGATAAGCCTGAAATGTCGCAATTGCCCGCGCCTGaaa                                                             >  2:753743/1‑47 (MQ=255)
                                                             aTAAGCCTGAAATGTCGCAATTGCCCTCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATCCCACAGTCCCAGGCAGTAAGAcc       <  2:152430/90‑1 (MQ=255)
                                                                          gTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATCACACAGTCCCAGGCAGTAAGACCGACGt  >  1:935830/1‑82 (MQ=255)
                                                                          gTCGCAATTGCCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATCACACAGTCCCAGGCAGTAAGACCGACGt  <  2:935830/82‑1 (MQ=255)
                                                                                    |                                                                       
GTGACGCGTCGGCGTAGTTTTGTATCTCCAGCCCGGGTAACCAGCGTACCGGATGCCCGAGATAAGCCTGAAATGTCGCAATTGTCCGCGCCTGAAAGCCAATAAAACCAATGATTTGGCTAATGATCACACAGTCCCAGGCAGTAAGACCGACGT  >  NZ_CP009273/1343755‑1343910

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