Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,596,499 C→A 9.0% pseudogene (1489/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,596,4990CA 9.0% 111.0 / 2.8 51pseudogene (1489/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base A (3/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

TATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATAC  >  NC_000913/1596352‑1596642
                                                                                                                                                   |                                                                                                                                               
tATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCa                                                                                                                                                >  1:69753/1‑149 (MQ=255)
  tCACGCGCTTTAACTATACCTCCGGTTGTAATATTCAGTTCAACCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCAtt                                                                                                                                              <  2:178896/149‑1 (MQ=255)
     cgcgCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAAATATCCCCTTACCGTAGGCATTtcc                                                                                                                                           >  2:939397/1‑149 (MQ=255)
      gcgcTTTAACTATACCTACGGTTGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTAAGTAGATGATGTCTTTAAATTCCAGAGACTCTCCGTTGAAATATTAACTATCCCCTTACCTGAGTCATTtcct                                                                                                                                          >  2:6035/1‑149 (MQ=255)
       cgcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGAGTGTCCGTTGAAATATTAACTTTCCCCTTACCGGAGGCATTtcctc                                                                                                                                         >  2:1023901/1‑149 (MQ=255)
       cgcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcctc                                                                                                                                         <  2:297640/149‑1 (MQ=255)
       cgcTTTAACTATACCTCCGCTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGAAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCAATtcctc                                                                                                                                         <  1:157038/149‑1 (MQ=255)
        gcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcctcc                                                                                                                                        >  2:790039/1‑149 (MQ=255)
          tttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTGTAAATTCCAGAGACTGTCCGTTGATATATTAACTTTCCCCTTACCGGAGGCATTTCCTCCGa                                                                                                                                      >  1:827876/1‑149 (MQ=255)
           ttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGAt                                                                                                                                     <  1:964786/149‑1 (MQ=255)
            tAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTACTACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGATATATTAACTATCCCCTTACCGGAGGCAGTTCCTACGATa                                                                                                                                    >  1:1003483/1‑149 (MQ=255)
              aCTATACCTCCGGTGGTAATATTCAGTTCACCCGTTCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGAAATGCCTCCGGTaag                                                                                                                                  <  1:239202/149‑2 (MQ=255)
              aCTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATACCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTATCGGAGGAAATGCCTCCGGTaag                                                                                                                                  <  1:227052/149‑2 (MQ=255)
                    ccTCCGGTGGTAATTTTCAGTTCACCCGTTCCTAATACACCGAACTGTATGAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGTCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGTCATTTCCTCCGATAAATCCATa                                                                                                                            <  2:636822/149‑1 (MQ=255)
                       ccGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGcc                                                                                                                         <  2:223061/149‑1 (MQ=255)
                           tggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGtt                                                                                                                     <  2:651850/149‑1 (MQ=255)
                             gtAATAGTCAGATCACCCGTGCCTAATACACCGACCTGTAGCAATTGTTCGATAGTACATGATCTCATTAAATTCCAGAGACTGTACGTTGAAATCTTAACTAACCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGa                                                                                                                   <  2:272029/149‑1 (MQ=255)
                              tAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAATTCCATAGCCCGTTgta                                                                                                                  >  2:877754/1‑147 (MQ=255)
                                 tatTCAGTTCACCCGTGCCTAATACACCGACATGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCAAGAGACTGTCAGTTGAAATATTAACTATCCCCTTACCGTAGGAATTTCCTCCGATAAAACAAAAGCCAGTTGAATTa                                                                                                               >  2:999540/1‑149 (MQ=255)
                                   ttCAGTTCACCCGTGCCTAAAACCCCGACCTGAAGCAATTGTGCGTTAGGAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAATTCCAGTGCCCGTTGAATTAAg                                                                                                             >  1:972579/1‑149 (MQ=255)
                                      aGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGt                                                                                                          <  2:139587/149‑1 (MQ=255)
                                           aCCCGTGCCTAATACACCGACATGTAGCAATTGTGCGTTAGAAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTAACCCCTTACCGGAGGCATTTCCTGCGATAAATCCCTAGCCCGTTGAATTAAGGGTgtccc                                                                                                     >  2:882843/1‑145 (MQ=255)
                                               gTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGGTGAAATATTATCTATACCCTTACCGGAGGCATTTCCTCCGATTAATCCTGAGCCCGGTGAAGTTCGGGTGGCGCCAtt                                                                                                 >  1:1021348/1‑149 (MQ=255)
                                                   cTAATACACCGCCCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAATTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCTATAAATCCATAGCCCGTTTAATTAAGGGTGGCACCATTACTg                                                                                             >  1:278922/1‑149 (MQ=255)
                                                    tAATACACCGCCCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTACAGAGACTGTCCGATGAAAAATTAACTATCCCCTTACCGGAGGAATTTCCTCCGATAAATCAATAGCCCGTTGCATTAAGGGTGACACCATTACTGa                                                                                            >  2:706356/1‑149 (MQ=255)
                                                        acacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACAtt                                                                                        >  2:113101/1‑149 (MQ=255)
                                                         cacCGACCTGTAGCAATTGTTCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGAAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACAttt                                                                                       >  2:325199/1‑149 (MQ=255)
                                                          acCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTa                                                                                      <  1:214859/149‑1 (MQ=255)
                                                           ccGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTAc                                                                                     <  1:959977/149‑1 (MQ=255)
                                                                    aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCt                                                                            <  2:37531/149‑1 (MQ=255)
                                                                      cAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTAc                                                                          >  1:652847/1‑149 (MQ=255)
                                                                            tgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCCTTGAATTAAGGGTGGCACCATAACTGACATTTACGATACCGCTACCAAAAt                                                                    >  1:294127/1‑149 (MQ=255)
                                                                            tgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTACGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAAt                                                                    <  2:995877/149‑1 (MQ=255)
                                                                            ggCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGTAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAAt                                                                    <  1:939397/148‑1 (MQ=255)
                                                                               gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGAGGCACCATTACTGACATATACGATACCGCTACCAAAAttcc                                                                 >  2:224982/1‑147 (MQ=255)
                                                                                    tagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                            >  1:988940/1‑149 (MQ=255)
                                                                                    tagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACAGACa                                                            >  2:889152/1‑149 (MQ=255)
                                                                                         gatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTTACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAt                                                       <  1:790039/149‑1 (MQ=255)
                                                                                                ttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcaca                                                >  2:911705/1‑149 (MQ=255)
                                                                                                ttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcaca                                                <  2:988940/149‑1 (MQ=255)
                                                                                                     ttCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGAAATTTACGATACCGATACCAAAATTACCGACATTGATACCCACAGTATc                                           >  1:943015/1‑149 (MQ=255)
                                                                                                                       gAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTTCCGACATTGATACCCACAGTATCGAACGATGATGTCGcatc                         <  2:203982/149‑1 (MQ=255)
                                                                                                                       gAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatc                         <  2:69753/149‑1 (MQ=255)
                                                                                                                        aaaTATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAACCACAGTATCGAACGATGATGTCGcatca                        >  1:988012/1‑148 (MQ=255)
                                                                                                                             ttAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAg                   >  2:27490/1‑149 (MQ=255)
                                                                                                                                   aTCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAg             <  1:873065/149‑1 (MQ=255)
                                                                                                                                    tCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAgt            <  1:113101/149‑1 (MQ=255)
                                                                                                                                        cTTACCGGAGGCATTTCCTCCTATAAATCCATAGCCCGTTGAATTTAGGGTGGCACCATTACTGAGATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTACCGATCGATGATGTCGCATCCATCAGGTTCAGTGTGc        <  1:6035/149‑1 (MQ=255)
                                                                                                                                           atcGGAGGAAATGCCTCCGGTAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATAACCGACATTGATACCCACAGTATCGAACGATCAGGTCGCATCCATCAGGTTCAGTGTGCCGa     >  1:985743/3‑149 (MQ=255)
                                                                                                                                            ccGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGAt    <  1:139590/149‑1 (MQ=255)
                                                                                                                                              gtaGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTATGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATAc  <  1:999540/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATAC  >  NC_000913/1596352‑1596642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑