Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,597,313 (C)7→8 pseudogene (675/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,597,3061.C98.0% 174.0 / ‑2.9 50pseudogene (682/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (26/23);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

CAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCG  >  NC_000913/1597157‑1597451
                                                                                                                                                      |                                                                                                                                                  
cAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa‑‑c                                                                                                                                                  >  2:294841/1‑151 (MQ=255)
     cATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cccccc                                                                                                                                             >  1:163713/1‑151 (MQ=255)
         attGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cccccccctc                                                                                                                                         <  1:533894/151‑4 (MQ=255)
          ttGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cccccccctcg                                                                                                                                        >  2:98140/1‑147 (MQ=255)
           tGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cccccccctcgc                                                                                                                                       >  2:245702/1‑146 (MQ=255)
            ggAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTcgcg                                                                                                                                       >  1:193386/1‑151 (MQ=255)
            ggAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTcgcg                                                                                                                                       >  1:592155/1‑151 (MQ=255)
               aaTATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGaa                                                                                                                                     >  1:756887/1‑150 (MQ=255)
               aaTATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAAt                                                                                                                                    >  1:179522/1‑151 (MQ=255)
               aaTATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAAt                                                                                                                                    >  2:742573/1‑151 (MQ=255)
                   ttGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACAGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATa                                                                                                                                   >  1:236167/1‑148 (MQ=255)
                      aCTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGc                                                                                                                              >  2:814599/1‑150 (MQ=255)
                         gtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTcc                                                                                                                          <  2:345927/151‑1 (MQ=255)
                              aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGt                                                                                                                     <  1:748792/151‑1 (MQ=255)
                                 aCCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTc                                                                                                                   <  1:976006/150‑1 (MQ=255)
                                  ccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                                                                                                                 <  1:705262/151‑1 (MQ=255)
                                   cGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                                                                                >  1:788256/1‑151 (MQ=255)
                                      aaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAt                                                                                                             <  2:755315/151‑1 (MQ=255)
                                        tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTgtt                                                                                                             >  2:737143/1‑147 (MQ=255)
                                            cATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                                                                                                       <  1:741234/151‑1 (MQ=255)
                                                     gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTCTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt                                                                                              >  2:844568/1‑151 (MQ=255)
                                                          cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                                                                         <  2:842665/151‑1 (MQ=255)
                                                          cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                                                                         <  2:575197/151‑1 (MQ=255)
                                                                gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa                                                                                   <  1:599210/151‑1 (MQ=255)
                                                                   atgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATc                                                                                 >  2:809348/1‑150 (MQ=255)
                                                                     gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                                                                              <  2:37438/151‑1 (MQ=255)
                                                                      aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                                                                              >  2:814497/1‑150 (MQ=255)
                                                                        agagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattatt                                                                           >  2:848747/1‑151 (MQ=255)
                                                                          agTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttt                                                                         <  2:444541/151‑1 (MQ=255)
                                                                               tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                                                                       <  2:725825/148‑1 (MQ=255)
                                                                                gATCCTGCTCATTCAGGGTTCCGATACCGGTCGCTTTGCCACCGATAATCGTATTTTCAGCCTTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatta                                                                   <  2:210061/151‑1 (MQ=255)
                                                                                gATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                                                                    >  1:216931/1‑150 (MQ=255)
                                                                                 aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                                                                    >  1:397511/1‑149 (MQ=255)
                                                                                 aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg                                                                  <  2:156881/151‑1 (MQ=255)
                                                                                       cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTc                                                            <  2:401305/151‑1 (MQ=255)
                                                                                           ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGcc                                                         <  2:521353/150‑1 (MQ=255)
                                                                                            tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc                                                       >  1:671109/1‑151 (MQ=255)
                                                                                              aGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCCTCGCTAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACCCGCCAcc                                                      >  2:780613/1‑151 (MQ=255)
                                                                                               gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCt                                                    >  1:249410/1‑151 (MQ=255)
                                                                                                    ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTc                                                <  1:26782/150‑1 (MQ=255)
                                                                                                      gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAAc                                             <  2:852284/151‑1 (MQ=255)
                                                                                                      gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAAc                                             >  2:498200/1‑151 (MQ=255)
                                                                                                           cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACaa                                        <  2:703422/151‑1 (MQ=255)
                                                                                                           cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACaa                                        <  1:294875/151‑1 (MQ=255)
                                                                                                            ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAc                                       >  2:354050/1‑151 (MQ=255)
                                                                                                             gTGGCATTGCCACCCATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAc                                       <  1:35184/150‑1 (MQ=255)
                                                                                                                gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg                                   >  1:271017/1‑151 (MQ=255)
                                                                                                                  aTTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGcc                                 <  1:915485/151‑1 (MQ=255)
                                                                                                                            aTAATCGTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTCGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGccc                       >  2:925635/1‑151 (MQ=255)
                                                                                                                                  gTATTTTCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTg                  >  2:79946/1‑150 (MQ=255)
                                                                                                                                  gTATTTTCAGCCGTAACCAA‑‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGAt                <  1:468713/151‑1 (MQ=255)
                                                                                                                                     ttttCAGCCGTAACCAACCCCCCCCTCGCGAATAGGCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCCTTATTTTTTATTAGCCACTTCACACCCTTTTCAACGAAAACCTGGCCATTACTGCCCCCCTGAt                >  2:659248/1‑149 (MQ=255)
                                                                                                                                       ttCAGCCGTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGcc            >  2:116230/1‑151 (MQ=255)
                                                                                                                                              gTAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAAt      >  2:609203/1‑150 (MQ=255)
                                                                                                                                               tAACCAACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCTCCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATg     <  2:144903/150‑1 (MQ=255)
                                                                                                                                                aaccaaccCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGc   <  1:98156/151‑1 (MQ=255)
                                                                                                                                                 accaaccCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGcg  <  2:438684/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCG  >  NC_000913/1597157‑1597451

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