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.C96.7% 103.8 / ‑3.9 30pseudogene (3740/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (13/16);  total (13/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

GTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑‑CCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAGTGGAGATGGTGAGC  >  NC_000913/1594106‑1594391
                                                                                                                                               |                                                                                                                                                 
gTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑‑cccccccc                                                                                                                                         >  1:629275/1‑151 (MQ=255)
    tAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑‑cccccccccacg                                                                                                                                     >  2:844384/1‑147 (MQ=255)
        gACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGatat                                                                                                                                  <  2:430687/151‑1 (MQ=255)
          cacGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATcc                                                                                                                                <  2:177724/151‑1 (MQ=255)
            cGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATcc                                                                                                                                >  1:780964/1‑149 (MQ=255)
                  gCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTAcc                                                                                                                        >  2:554440/1‑151 (MQ=255)
                       ccccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGAt                                                                                                                   <  2:593738/151‑1 (MQ=255)
                         ccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATAc                                                                                                                 <  2:803808/151‑1 (MQ=255)
                           aGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATAcgc                                                                                                               <  2:265286/151‑1 (MQ=255)
                              cATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATAcgccgc                                                                                                            >  1:698910/1‑151 (MQ=255)
                               aTAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCcgcg                                                                                                           >  2:210171/1‑151 (MQ=255)
                                       cATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgtt                                                                                                   <  1:646179/151‑1 (MQ=255)
                                           tCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                                                                >  2:780055/1‑150 (MQ=255)
                                             aCCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                                                                               >  2:579954/1‑150 (MQ=255)
                                               cTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAag                                                                                             >  1:229325/1‑149 (MQ=255)
                                                    gtggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAgg                                                                                      <  1:17344/151‑1 (MQ=255)
                                                                   gaaataccacCTCTGTGATGGATATTAATGTGTTTACTTGTCATCTTTATTGCAGAGTTGGAATCGAAAGCAATGT‑‑CCTCCCCCCACGATATCCGTCGTACCTTGATCCGCCTCGGCACCGTTGTTAAGAGTCTTGATCCTGGTaatcatc                                                                       <  1:758348/144‑7 (MQ=255)
                                                                   gACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAccc                                                                          >  2:442817/1‑149 (MQ=255)
                                                                        cacaCCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGt                                                                  <  2:119979/151‑1 (MQ=255)
                                                                           acCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcac                                                               <  2:717163/151‑1 (MQ=255)
                                                                            ccGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca                                                              <  2:659364/151‑1 (MQ=255)
                                                                                tGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCac                                                          >  2:456207/1‑151 (MQ=255)
                                                                                 gTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCac                                                          <  1:509294/150‑1 (MQ=255)
                                                                                    aTGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgct                                                      >  1:618287/1‑151 (MQ=255)
                                                                                             aaTGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGAt                                             <  1:608761/151‑1 (MQ=255)
                                                                                                                  aaTGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGc                         <  2:428755/150‑1 (MQ=255)
                                                                                                                  aaTGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCt                        >  1:719328/1‑151 (MQ=255)
                                                                                                                     ggCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACg                     <  2:249859/151‑1 (MQ=255)
                                                                                                                      gCAGAGTCGGAACCGAAAGCCATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTAc                      <  1:184722/149‑1 (MQ=255)
                                                                                                                           gTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAGt               <  1:844326/151‑1 (MQ=255)
                                                                                                                             cGGAACCGAAAGCAATTTCCCCCCCCCCCGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAGTgg             >  1:307323/1‑151 (MQ=255)
                                                                                                                                          aaTTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAGTGGAGATGGTGAGc  >  1:340892/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑‑CCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAGTGGAGATGGTGAGC  >  NC_000913/1594106‑1594391

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