Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,256 (C)8→7 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,2490C.98.3% 268.7 / ‑2.6 61pseudogene (3739/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base . (30/29);  minor base A (0/1);  total (30/31)
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

CACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTA  >  NC_000913/1594150‑1594340
                                                                                                   |                                                                                           
cACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTc                                                                                             >  2:1395084/1‑100 (MQ=255)
cACCTGCACTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTCATGTGTTGACTTGCCACATTAATGGCAGAGTCGGAACCGAAAGCAATTTc                                                                                             >  1:32944/1‑100 (MQ=255)
      naGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccccc                                                                                       >  1:180026/2‑100 (MQ=255)
        gtggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTACTGGCAGAGTCGGAACCGAAAGCAATTTcccccccac                                                                                     >  2:1274412/1‑98 (MQ=255)
        gtggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccac                                                                                     >  2:1349930/1‑98 (MQ=255)
         tggtgACATTTCCCGACATCACACCGCTGTTATTGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccacg                                                                                    <  2:1443796/100‑4 (MQ=255)
          ggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccacga                                                                                   >  2:1684442/1‑96 (MQ=255)
               cATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATcc                                                                             >  2:1534574/1‑100 (MQ=255)
                 tttCCCGATATCATACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCcg                                                                            <  2:586886/99‑1 (MQ=255)
                 tttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCcgt                                                                           >  1:22045/1‑100 (MQ=255)
                  ttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCcgtc                                                                          <  2:233325/100‑1 (MQ=255)
                   tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCcgtcg                                                                         >  1:1563325/1‑100 (MQ=255)
                    cccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCcgtcgt                                                                        >  2:812344/1‑100 (MQ=255)
                        aCATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTAccc                                                                    <  2:498748/100‑1 (MQ=255)
                        aCATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTAccc                                                                    <  1:501302/100‑1 (MQ=255)
                         cATCACACCGCTGTTCTGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTAccc                                                                    <  1:1062915/99‑1 (MQ=255)
                          aTCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTg                                                                  <  2:252782/100‑1 (MQ=255)
                             acacCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATa                                                               <  1:921183/100‑1 (MQ=255)
                             acacCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATa                                                               >  1:1503980/1‑100 (MQ=255)
                              cacCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATAc                                                              >  2:593147/1‑100 (MQ=255)
                                    tGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCcgcg                                                        <  2:1041076/100‑1 (MQ=255)
                                      ttATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCgg                                                       <  2:1452214/100‑1 (MQ=255)
                                        aTGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                     >  2:995798/1‑99 (MQ=255)
                                            atatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgtt                                                >  2:1118055/1‑100 (MQ=255)
                                            atatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgtt                                                <  1:633031/100‑1 (MQ=255)
                                             tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg                                               >  2:1491734/1‑100 (MQ=255)
                                             tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg                                               >  1:1599123/1‑100 (MQ=255)
                                             tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg                                               >  2:1351141/1‑100 (MQ=255)
                                               ttAATGTGTTGACTTGTCGTATTAATGGCAGAGTCGGAATCGAAAGCAATTT‑CCCCCCCACGCTATCCGGCGTACCCTGATACGCCGCGGAATCgttgtt                                             <  2:630740/100‑1 (MQ=255)
                                               ttAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGTAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                             <  2:463826/100‑1 (MQ=255)
                                                tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                            >  1:1023530/1‑100 (MQ=255)
                                                 aaTGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                            >  2:1786073/1‑99 (MQ=255)
                                                    gtgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAagag                                        >  2:928025/1‑100 (MQ=255)
                                                     tgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAagag                                        <  2:1721427/99‑1 (MQ=255)
                                                      gtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTc                                      >  1:1044822/1‑100 (MQ=255)
                                                      gtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAACCGTTGTTAAGAGTc                                      >  1:832114/1‑100 (MQ=255)
                                                         gACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAgg                                   >  1:571808/1‑100 (MQ=255)
                                                         gACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAgg                                   >  1:1049012/1‑100 (MQ=255)
                                                                cATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGGCGCGGAATCGTTGTTAAGAGTCAGGATCCCgg                            >  1:1378262/1‑100 (MQ=255)
                                                                cATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGCATCGTCGCTACGAGCCAGAATCCCgg                            >  2:1541910/1‑100 (MQ=255)
                                                                  tatTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCTCGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTa                          <  2:222715/100‑1 (MQ=255)
                                                                   atTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAc                         <  2:1641487/100‑1 (MQ=255)
                                                                       aTGGCAGAGTCGGAACCTAAAGCAATTT‑CCCACCCACGATATCCGTCGTACCCTGATACGCCGTGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAg                     <  2:1280080/100‑1 (MQ=255)
                                                                       aTGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAg                     >  1:1352705/1‑100 (MQ=255)
                                                                        tGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGCAATCGTTGTTTAGAGTCAGGAACCCGTAACCCAg                     >  2:1013032/1‑99 (MQ=255)
                                                                         ggCAGAGTCGGAACCGAAAGCAATTT‑ACCCCCCACGATATCCGTCGTACCCTGATACGCAGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTCCCCAGtt                   >  1:592538/1‑100 (MQ=255)
                                                                            agagTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTg                <  1:182451/100‑1 (MQ=255)
                                                                            agagTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTCATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttttn                <  1:1349920/100‑2 (MQ=255)
                                                                            agagTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCCTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTg                <  1:1274401/100‑1 (MQ=255)
                                                                              agTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTc              <  1:741208/100‑1 (MQ=255)
                                                                              agTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTc              <  1:1089677/100‑1 (MQ=255)
                                                                               gTCGGAACCGAAAGCAATTT‑CCCCCCCTCGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTca             <  2:571812/100‑1 (MQ=255)
                                                                               gTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTca             >  1:485155/1‑100 (MQ=255)
                                                                               gTCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTca             <  1:743970/100‑1 (MQ=255)
                                                                                tCGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcac            <  2:1352325/100‑1 (MQ=255)
                                                                                 cGGAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca           >  2:1301469/1‑100 (MQ=255)
                                                                                 cGGAACCGAAAGCAA‑TTACCCCCCCACGAGATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca           <  1:492166/100‑1 (MQ=255)
                                                                                   gAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATc         >  2:1420780/1‑100 (MQ=255)
                                                                                   gAACCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATc         <  2:705272/100‑1 (MQ=255)
                                                                                     aCCGAAAGCAATTT‑CCTCCCCATGTTATCCGTCGTACCCTGATACGCCGCGGAATCGTGGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACAt          <  2:1523808/97‑1 (MQ=255)
                                                                                     aCCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCac       <  1:1491723/100‑1 (MQ=255)
                                                                                     aCCGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCac       <  2:1599134/100‑1 (MQ=255)
                                                                                      ccaaaaTCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCacc      <  2:1092414/97‑1 (MQ=255)
                                                                                      ccGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCacc      <  1:675777/100‑1 (MQ=255)
                                                                                       cGAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccg     >  1:1501047/1‑100 (MQ=255)
                                                                                        gAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgc    >  2:592302/1‑100 (MQ=255)
                                                                                        gAAAGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgc    >  1:458159/1‑100 (MQ=255)
                                                                                          aaGCAATTT‑CCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgcta  >  1:780892/1‑100 (MQ=255)
                                                                                                   |                                                                                           
CACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTA  >  NC_000913/1594150‑1594340

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑