Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,256 (C)8→9 pseudogene (3732/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,594,2481.C89.8% 198.6 / 13.4 59pseudogene (3740/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base C (27/26);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑CCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAG  >  NC_000913/1594111‑1594377
                                                                                                                                             |                                                                                                                                   
aGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑cc                                                                                                                                  >  2:1133228/1‑139 (MQ=255)
    acacGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑cccccc                                                                                                                              >  2:491539/1‑139 (MQ=255)
               cccccccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGatat                                                                                                                    <  2:417699/139‑1 (MQ=255)
                  ccccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCATACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCcg                                                                                                                 <  1:207550/139‑1 (MQ=255)
                  ccccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCcg                                                                                                                 >  1:403545/1‑139 (MQ=255)
                            aaCGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTg                                                                                                       >  1:171891/1‑139 (MQ=255)
                                         ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                                                          <  1:939288/139‑1 (MQ=255)
                                         ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑CCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                                                          >  2:592369/1‑138 (MQ=255)
                                         ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑CCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                                                          <  1:592369/138‑1 (MQ=255)
                                           tGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATc                                                                                        <  1:491539/139‑1 (MQ=255)
                                           tGCAGTGGTGACATTTCCCGACATCACACCGCTGTCATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATc                                                                                        <  1:13691/139‑1 (MQ=255)
                                             cAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgt                                                                                      <  2:403545/139‑1 (MQ=255)
                                                    gACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAag                                                                               <  1:670478/139‑1 (MQ=255)
                                                          tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                                                          <  2:757193/122‑1 (MQ=255)
                                                          tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                                                          >  1:757193/1‑122 (MQ=255)
                                                          tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAg                                                                         <  2:956418/139‑1 (MQ=255)
                                                               aCATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATcc                                                                    <  1:311762/139‑1 (MQ=255)
                                                                   cacaCCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑cccc                                                                                                                                >  2:953600/1‑75 (MQ=255)
                                                                   cacaCCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑cccc                                                                                                                                <  1:953600/75‑1 (MQ=255)
                                                                     cacCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAc                                                              >  2:830646/1‑139 (MQ=255)
                                                                          cTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGt                                                         >  1:613974/1‑139 (MQ=255)
                                                                           tGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtt                                                        >  1:702762/1‑139 (MQ=255)
                                                                             tttttggtttttaTGTGTTTACTTTCCATTTTAATGGCCAGATTCGAACCCGAAACAAATTCCC‑‑CCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtt                                                        <  2:1258536/127‑1 (MQ=255)
                                                                             ttATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑cccccc                                                                                                                              >  1:1284754/1‑67 (MQ=255)
                                                                             ttATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑cccccc                                                                                                                              <  2:1284754/67‑1 (MQ=255)
                                                                                    tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑CCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCa                                                           <  2:773589/126‑1 (MQ=255)
                                                                                    tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑CCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCa                                                           >  1:773589/1‑126 (MQ=255)
                                                                                     atTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATc                                              >  1:716417/1‑139 (MQ=255)
                                                                                       tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATAcg                                                                                                  <  2:1427736/85‑1 (MQ=255)
                                                                                       tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATAcg                                                                                                  >  1:1427736/1‑85 (MQ=255)
                                                                                          tgtgTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGat                                                                                                                      <  2:519388/62‑1 (MQ=255)
                                                                                          tgtgTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGat                                                                                                                      >  1:519388/1‑62 (MQ=255)
                                                                                          tgtgTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgc                                         >  1:805321/1‑139 (MQ=255)
                                                                                           gtgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgct                                        <  2:613974/139‑1 (MQ=255)
                                                                                             gtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATTCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                                                  >  1:1078268/1‑95 (MQ=255)
                                                                                             gtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctac                                      <  1:369951/139‑1 (MQ=255)
                                                                                             gtTGACTTGCCATATTAAGGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑CCCCCCCGCACGATATTCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                                                  <  2:1078268/95‑1 (MQ=255)
                                                                                              ttGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                                                                 <  2:336807/95‑1 (MQ=255)
                                                                                              ttGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                                                                 >  1:336807/1‑95 (MQ=255)
                                                                                                     gCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                                                          >  1:201477/1‑79 (MQ=255)
                                                                                                     gCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                                                          <  2:201477/79‑1 (MQ=255)
                                                                                                      ccATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                                                  >  1:1416809/1‑86 (MQ=255)
                                                                                                      ccATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                                                  <  2:1416809/86‑1 (MQ=255)
                                                                                                      ccATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctac                                      >  2:822838/1‑130 (MQ=255)
                                                                                                      ccATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctac                                      <  1:822838/130‑1 (MQ=255)
                                                                                                        atatTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACtg                           >  1:1415031/1‑139 (MQ=255)
                                                                                                          atTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGTTACCCAGTTTTGTCACATCACCGCTACCGCTGATACtgtg                         <  1:1174182/139‑1 (MQ=255)
                                                                                                              aTGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                                                  >  2:884038/1‑78 (MQ=255)
                                                                                                              aTGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                                                  <  1:884038/78‑1 (MQ=255)
                                                                                                               tGGCAGAGTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCa                    >  2:1139243/1‑139 (MQ=255)
                                                                                                                      agTCGGAACCGCAAGCAATT‑‑TCCCCCCCCCAGGCAATCCGCCGTACCCTGATACGCCGAGGAAGCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCCTCATCGccc              >  2:133756/1‑139 (MQ=255)
                                                                                                                      agTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGcc               >  2:1083807/1‑139 (MQ=255)
                                                                                                                      agTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAaga                                                                              >  1:407139/1‑75 (MQ=255)
                                                                                                                      agTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAaga                                                                              <  2:407139/75‑1 (MQ=255)
                                                                                                                      agTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAcc                                                             >  1:203125/1‑92 (MQ=255)
                                                                                                                      agTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAcc                                                             <  2:203125/92‑1 (MQ=255)
                                                                                                                      agTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTc                     <  2:702471/132‑1 (MQ=255)
                                                                                                                      agTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTc                     >  1:702471/1‑132 (MQ=255)
                                                                                                                      agTCGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGccc              >  2:1289977/1‑139 (MQ=255)
                                                                                                                         cGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctac                                      <  2:526401/112‑1 (MQ=255)
                                                                                                                         cGGAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctac                                      >  1:526401/1‑112 (MQ=255)
                                                                                                                           gAACCGAAAGCAATT‑‑TCCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCt          <  1:1083807/139‑1 (MQ=255)
                                                                                                                           gAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctac                                      <  1:1125143/110‑1 (MQ=255)
                                                                                                                           gAACCGAAAGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctac                                      >  2:1125143/1‑110 (MQ=255)
                                                                                                                                  aaGCAATT‑‑TCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAg  >  2:153115/1‑139 (MQ=255)
                                                                                                                                             |                                                                                                                                   
AGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATT‑‑T‑‑CCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAG  >  NC_000913/1594111‑1594377

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