Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,597,313 (C)7→9 pseudogene (774/5422 nt) yneO ← AIDA‑I family autotransporter YneO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,597,3061.C97.1% 113.2 / ‑3.2 35pseudogene (781/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (19/15);  total (20/15)
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.98e-01
*NC_0009131,597,3062.C97.1% 116.0 / ‑3.2 35pseudogene (781/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (19/15);  total (20/15)
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.98e-01

TGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑CCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGC  >  NC_000913/1597168‑1597432
                                                                                                                                           ||                                                                                                                                
tGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa‑‑c                                                                                                                                 <  1:506987/140‑1 (MQ=255)
     atatTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cccc‑‑cc                                                                                                                          <  1:499336/140‑1 (MQ=255)
                gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAAt                                                                                                                 >  2:144999/1‑140 (MQ=255)
                gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATAGTATTTTGAGCCGTTACCAACGCCCC‑‑CCCTCGCGaca                                                                                                                 >  1:421453/1‑138 (MQ=255)
                 ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATa                                                                                                                >  2:619738/1‑140 (MQ=255)
                    ttACCGAAATATCCATTGTAGACTCTGCGTAACGTGATGGCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTc                                                                                                             <  2:294187/140‑1 (MQ=255)
                           aaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGACGTAACCAA‑‑CCCCCACCCTCGCGAATAGTCGCCTccc                                                                                                      >  1:294187/1‑140 (MQ=255)
                           aaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTATCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTccc                                                                                                      <  1:84011/140‑1 (MQ=255)
                               tCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAGTAATCCCCCCCCCAgg                                                                                                  >  1:196756/1‑139 (MQ=255)
                                    tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTc                                                                                                             >  2:60078/1‑124 (MQ=255)
                                    tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTc                                                                                                             <  1:60078/124‑1 (MQ=255)
                                         agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTccc                                                                                                      <  1:205589/126‑1 (MQ=255)
                                         agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTccc                                                                                                      >  2:205589/1‑126 (MQ=255)
                                         agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTa                                                                                        >  1:504655/1‑140 (MQ=255)
                                           tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                                                      >  2:430294/1‑140 (MQ=255)
                                           tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                                                      >  1:776620/1‑140 (MQ=255)
                                             tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTTGCATCCCCATTTCCTTGATTACCaa                                                                                    >  1:418249/1‑140 (MQ=255)
                                             tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                                                                                    >  2:125156/1‑140 (MQ=255)
                                             tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                                                                                    >  1:513965/1‑140 (MQ=255)
                                                           aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                                                      >  2:646735/1‑140 (MQ=255)
                                                              gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                                                   <  2:739311/140‑1 (MQ=255)
                                                                           gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttt                                                      >  1:171883/1‑140 (MQ=255)
                                                                           gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttt                                                      >  2:29987/1‑140 (MQ=255)
                                                                             acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                                                    >  2:568079/1‑140 (MQ=255)
                                                                               aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                                                  <  2:776620/140‑1 (MQ=255)
                                                                                 tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTGCCCAGTTCCTTGATTACCAATTTGCAATTCAATTGCGGAATCATTATTTTttattt                                                >  2:185002/1‑139 (MQ=255)
                                                                                 tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatta                                                >  2:69224/1‑140 (MQ=255)
                                                                                 tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTACTTGATTACCAAGTTGAAATTCAATCGAGGAATCATTATTTTttatta                                                >  2:280004/1‑140 (MQ=255)
                                                                                     ggTTCCGAAACCGGTGGCATTGACACCGATAGTGGAATTTTCAGCCGTGACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTGCCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCa                                            <  2:418249/140‑1 (MQ=255)
                                                                                      gTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCAc                                           >  2:521578/1‑140 (MQ=255)
                                                                                                  gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCttt                               <  1:44367/140‑1 (MQ=255)
                                                                                                        ttGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACg                         <  1:521578/140‑1 (MQ=255)
                                                                                                             aCCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAc                    <  2:171883/140‑1 (MQ=255)
                                                                                                             aCCGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAc                    <  1:87685/140‑1 (MQ=255)
                                                                                                               cGATAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCt                  <  1:843402/140‑1 (MQ=255)
                                                                                                                 aTAATCGTATTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg                <  2:504655/140‑1 (MQ=255)
                                                                                                                       gtttttcCAGCGTTAACAAC‑‑CCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACt        <  1:534613/137‑1 (MQ=255)
                                                                                                                         tttttCAGTCGTAACCGACCCCCC‑‑CCCTCCTCATTCGTCGCCTCCCCATTTCCTTGATTACCAACTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCTCCCTTTTCAACGACAACCTGGCCATTACt        <  2:421453/139‑1 (MQ=255)
                                                                                                                         aTTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACt        <  1:619738/140‑1 (MQ=255)
                                                                                                                         aTTTTCAGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACt        <  1:646735/140‑1 (MQ=255)
                                                                                                                               aGCCGTAACCAACCCCCC‑‑CCCTCGCGAATAGTCGCCGCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTCTTTATTAGCCACTCGCCACCCTTTTCAACGCCAACCTCCCCATTACTGCCCGc  >  2:195606/1‑140 (MQ=255)
                                                                                                                                           ||                                                                                                                                
TGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑CCCC‑‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGC  >  NC_000913/1597168‑1597432

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