Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,877,816 T→C pseudogene (408/1482 nt) yeaV → BCCT family transporter YeaV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,877,8160TC96.5% 90.5 / ‑3.5 29pseudogene (408/1482 nt)yeaVBCCT family transporter YeaV
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (15/13);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

GGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTA  >  NZ_CP009273/1877729‑1877901
                                                                                       |                                                                                     
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                     >  1:335899/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                     >  1:170877/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                     >  2:601247/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                     >  2:299408/1‑90 (MQ=255)
ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                                     >  1:896145/1‑90 (MQ=255)
     cATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTccgcccgc                                                                                >  2:894031/1‑90 (MQ=255)
       ttAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCaa                                                                              <  1:767005/90‑1 (MQ=255)
       ttAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCaa                                                                              <  1:894031/90‑1 (MQ=255)
          aTCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGc                                                                           >  1:924999/1‑90 (MQ=255)
             aTGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCa                                                                        <  1:475779/90‑1 (MQ=255)
             aTGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCa                                                                        <  1:503561/90‑1 (MQ=255)
             aTGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCa                                                                        <  1:97618/90‑1 (MQ=255)
               ggCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATg                                                                      >  1:373501/1‑90 (MQ=255)
               ggCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATg                                                                      >  1:720728/1‑90 (MQ=255)
                                  cGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCAt                                                                       <  1:531325/70‑1 (MQ=255)
                                  cGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCAt                                                                       >  2:531325/1‑70 (MQ=255)
                                         aCAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTcc                                            <  2:474303/90‑1 (MQ=255)
                                         aCAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTcc                                            <  1:601247/90‑1 (MQ=255)
                                         aCAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTcc                                            <  2:779053/90‑1 (MQ=255)
                                              ggTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATc                                       >  2:766621/1‑90 (MQ=255)
                                                   gAGCCTTTCCGGCATTATTGCTACTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCAc                                  >  1:825440/1‑90 (MQ=255)
                                                        tttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCg                             >  2:633355/1‑90 (MQ=255)
                                                         ttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCgg                            >  1:576770/1‑90 (MQ=255)
                                                              gCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCAc                       <  1:995861/90‑1 (MQ=255)
                                                               cATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACt                      <  1:405791/90‑1 (MQ=255)
                                                                   attGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGAcc                  <  2:836681/90‑1 (MQ=255)
                                                                      gctgctATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCAtt               <  1:952368/90‑1 (MQ=255)
                                                                           tATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCt          <  2:372026/90‑1 (MQ=255)
                                                                                   gCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTa  >  1:455142/1‑90 (MQ=255)
                                                                                   gCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTa  >  2:605004/1‑90 (MQ=255)
                                                                                       |                                                                                     
GGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCTTGTTGTTA  >  NZ_CP009273/1877729‑1877901

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: