Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 716,388 C→T pseudogene (210/651 nt) ybfG ← pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913716,3880CT97.1% 108.4 / NA 36pseudogene (210/651 nt)ybfGpseudogene
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/16);  minor base A (1/0);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTG  >  NC_000913/716254‑716521
                                                                                                                                      |                                                                                                                                     
ccaTTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGtaattaac                                                                                                                                   <  1:746592/139‑1 (MQ=255)
     aaGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTa                                                                                                                              <  2:174502/139‑1 (MQ=255)
        agCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatca                                                                                                                           >  1:611603/1‑139 (MQ=255)
                       tGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTCGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                            <  1:554058/139‑1 (MQ=255)
                        gCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                                           <  1:840922/139‑1 (MQ=255)
                               cGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATa                                                                                                    <  2:353043/139‑1 (MQ=255)
                                       aTGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTACGGAGCCATAGACAAgcg                                                                                            >  2:107545/1‑139 (MQ=255)
                                              aaGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCAt                                                                                     >  2:152894/1‑139 (MQ=255)
                                              aaGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCAt                                                                                     <  1:152894/139‑1 (MQ=255)
                                                     cACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGtaattaac                                                                                                                                   <  2:652531/86‑1 (MQ=255)
                                                     cACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGtaattaac                                                                                                                                   >  1:652531/1‑86 (MQ=255)
                                                              gAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTc                                                                     <  2:330475/139‑1 (MQ=255)
                                                                         gTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCa                                                                                          >  2:242334/1‑107 (MQ=255)
                                                                         gTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCa                                                                                          <  1:242334/107‑1 (MQ=255)
                                                                                       aCTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                            >  1:500438/1‑75 (MQ=255)
                                                                                       aCTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAAGTAACCATTATCAATCAGCAGACCTCCg                                                                                                            <  2:500438/75‑1 (MQ=255)
                                                                                                cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCACACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAAt                                   <  1:254274/139‑1 (MQ=255)
                                                                                                cGTGGCAACCGGGGGTGACAAGATCCTGCTATGCATAAATAACCATTATCAATAAGCAGACCTCCGGAGCCATAGACAAGCGTACCCATCCGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCAAAAt                                   >  2:277817/1‑139 (MQ=255)
                                                                                                   ggCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                            >  1:655505/1‑63 (MQ=255)
                                                                                                   ggCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                            <  2:655505/63‑1 (MQ=255)
                                                                                                      aaCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTttat                             >  2:709149/1‑139 (MQ=255)
                                                                                                       aCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                            <  2:711762/59‑1 (MQ=255)
                                                                                                       aCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                            >  1:711762/1‑59 (MQ=255)
                                                                                                             ggTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTaa                          >  2:124328/1‑135 (MQ=255)
                                                                                                             ggTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTaa                          <  1:124328/135‑1 (MQ=255)
                                                                                                                      aTCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                                           >  1:635536/1‑45 (MQ=255)
                                                                                                                      aTCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                                           <  2:635536/45‑1 (MQ=255)
                                                                                                                      aTCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGttt                                            <  2:152687/108‑1 (MQ=255)
                                                                                                                      aTCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGttt                                            >  1:152687/1‑108 (MQ=255)
                                                                                                                      aTCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATCTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAGATCTTTATTAACCGGTAAGAGTTc             >  1:468303/1‑139 (MQ=255)
                                                                                                                                gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcg                                                                                            <  1:551392/50‑1 (MQ=255)
                                                                                                                                gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcg                                                                                            >  2:551392/1‑50 (MQ=255)
                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATa                                                                                                    <  2:266163/41‑1 (MQ=255)
                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATa                                                                                                    >  1:266163/1‑41 (MQ=255)
                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCt                                                                 >  2:424547/1‑76 (MQ=255)
                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCt                                                                 <  1:424547/76‑1 (MQ=255)
                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTtg  >  1:103970/1‑139 (MQ=255)
                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTtg  >  2:471397/1‑139 (MQ=255)
                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTtg  >  2:432509/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                     
CCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTG  >  NC_000913/716254‑716521

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