Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,381,779 (G)5→6 pseudogene (959/969 nt) ycjV → putative ABC transporter ATP‑binding protein YcjV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,381,7741.G91.7% 67.6 / 0.5 25pseudogene (954/969 nt)ycjVputative ABC transporter ATP‑binding protein YcjV
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (15/7);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

AGCCCATCGTGATGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTA  >  NC_000913/1381650‑1381901
                                                                                                                             |                                                                                                                                 
aGCCCATCATGATGGCAATACGCCCGGAAGCTATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTggggggca                                                                                                                          <  1:135040/133‑4 (MQ=255)
  cccATCGTGATGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAAGCGGTGCGGTATTGATGATCTACACCGCGGTgggggggcacg                                                                                                                        >  2:135040/1‑128 (MQ=255)
      tCGTGATGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGtt                                                                                                                    >  2:10992/1‑133 (MQ=255)
           aTGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTgg                                                                                                               >  1:79005/1‑133 (MQ=255)
           aTGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCAACACCACTGTTAGGGGGCACGGGTTAGTgg                                                                                                               >  1:225018/1‑133 (MQ=255)
             gggAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTc                                                                                                             >  2:93448/1‑133 (MQ=255)
                aaTACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGTAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGt                                                                                                          <  1:335367/133‑1 (MQ=255)
                 aTACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGAGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTg                                                                                                         >  2:490217/1‑133 (MQ=255)
                          gAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGt                                                                                                >  2:259605/1‑133 (MQ=255)
                          gAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGt                                                                                                >  2:12856/1‑133 (MQ=255)
                                  tCATCCGGACGCGCAGGAGGAAAATAACATATCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGAt                                                                                        <  1:209997/133‑1 (MQ=255)
                                           cgcgCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCa                                                                               >  2:1061/1‑133 (MQ=255)
                                            gcgcAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGGGCGATAGATGATTCTCATGCAg                                                                              >  1:235266/1‑133 (MQ=255)
                                              gcAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGga                                                                            <  1:1061/133‑1 (MQ=255)
                                                        aaTAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCa                                                                  >  1:339069/1‑133 (MQ=255)
                                                             cATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTAtt                                                             >  1:440741/1‑133 (MQ=255)
                                                              aTTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATAATTATCATGCATGATAAAATAGCACTATTc                                                            >  1:242549/1‑133 (MQ=255)
                                                                      cAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGat                                                    <  1:10434/133‑1 (MQ=255)
                                                                                                  gCGGAATTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAg                        >  2:338925/1‑133 (MQ=255)
                                                                                                      aaTTTATGCTCTACACCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAAc                    <  1:118045/133‑1 (MQ=255)
                                                                                                                  acacCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt        >  2:331784/1‑133 (MQ=255)
                                                                                                                  acacCACGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt        <  2:440741/133‑1 (MQ=255)
                                                                                                                   caccacGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt        <  1:93989/132‑1 (MQ=255)
                                                                                                                   caccacGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt        >  2:93989/1‑132 (MQ=255)
                                                                                                                       acGGTTGGGGCGCACGAGTTAGTGGTCCGTGCTGGTGCGTTGAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCt   <  2:225018/133‑1 (MQ=255)
                                                                                                                        cGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAACTGTCATTTGTTTGATGTAGAAACGGAAATATCAATTCGCTa  >  2:174500/1‑133 (MQ=255)
                                                                                                                        cGGTTGGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTa  >  2:314474/1‑133 (MQ=255)
                                                                                                                        cGGTTGGGGGACACGGGTTAGTGGCCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTa  <  2:242549/133‑1 (MQ=255)
                                                                                                                             |                                                                                                                                 
AGCCCATCGTGATGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTA  >  NC_000913/1381650‑1381901

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