Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,256 (C)8→10 pseudogene (3831/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,594,2481.C98.1% 187.4 / ‑3.6 53pseudogene (3839/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (30/22);  total (31/22)
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.97e-01
*NC_0009131,594,2482.C98.1% 186.5 / ‑3.6 53pseudogene (3839/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (30/22);  total (31/22)
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.97e-01

AGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑‑CCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCA  >  NC_000913/1594111‑1594376
                                                                                                                                           ||                                                                                                                                  
aGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑‑cc                                                                                                                                  >  1:488699/1‑140 (MQ=255)
aGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑‑ccc                                                                                                                                 <  2:561297/140‑1 (MQ=255)
 gCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑‑ccc                                                                                                                                 >  2:887933/1‑140 (MQ=255)
  cGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑‑cccc                                                                                                                                <  1:887933/140‑1 (MQ=255)
    ccccGCATTGCCCCCCCCAGGCTTAACGTTCTTATCACCTGCATTGGGGACATTCCCCGACACCCCCCGGCTGTTATGAAAATTAAGGTGTTGCCTTCCCAAATAAATGGCGGGTCGGGACCCGAAAGCAATTT‑‑cccccc                                                                                                                              <  2:431030/137‑1 (MQ=255)
      acGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAACCAATTT‑‑cccccc‑‑cc                                                                                                                          >  2:217249/1‑140 (MQ=255)
                   cccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCcgt                                                                                                               >  1:511607/1‑140 (MQ=255)
                    ccAGCCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCcgtc                                                                                                              >  1:258487/1‑140 (MQ=255)
                     cAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCcgtcg                                                                                                             <  1:217249/140‑1 (MQ=255)
                     cAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCcgtcg                                                                                                             <  1:586307/140‑1 (MQ=255)
                               gTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATa                                                                                                   <  1:405474/140‑1 (MQ=255)
                                      tCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCcgcg                                                                                            >  1:67334/1‑140 (MQ=255)
                                         ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                                                         <  2:576405/140‑1 (MQ=255)
                                         ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                                                         <  2:776177/140‑1 (MQ=255)
                                            gCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCg                                                                                      <  2:511607/140‑1 (MQ=255)
                                             cAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgt                                                                                     >  2:540226/1‑140 (MQ=255)
                                                 ggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                                                 <  2:909258/140‑1 (MQ=255)
                                                    gACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAag                                                                              >  2:377366/1‑140 (MQ=255)
                                                      cATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAagag                                                                            <  2:641656/140‑1 (MQ=255)
                                                       aTTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGt                                                                           >  2:871875/1‑140 (MQ=255)
                                                          tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                                                 <  2:905307/131‑1 (MQ=255)
                                                          tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                                                 >  1:905307/1‑131 (MQ=255)
                                                          tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAg                                                                        >  2:514237/1‑140 (MQ=255)
                                                            ccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGa                                                                      >  1:823978/1‑140 (MQ=255)
                                                            ccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGa                                                                      <  1:377366/140‑1 (MQ=255)
                                                            ccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGa                                                                      <  1:439487/140‑1 (MQ=255)
                                                                 aTCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCg                                                                 <  2:258487/140‑1 (MQ=255)
                                                                   cacaCCGCTGTTATGGATATTAATGTGTTGACTTGCTATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGt                                                               >  2:450093/1‑140 (MQ=255)
                                                                          cTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGt                                                        >  2:334655/1‑140 (MQ=255)
                                                                           tGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCGCGATATCCGTCGGACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtt                                                       >  2:800999/1‑140 (MQ=255)
                                                                               aTGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGt                                                   <  2:823978/140‑1 (MQ=255)
                                                                                tGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTc                                                  <  2:488699/140‑1 (MQ=255)
                                                                                   atatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca                                               <  1:540226/140‑1 (MQ=255)
                                                                                    tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACAt                                              >  2:458118/1‑140 (MQ=255)
                                                                                       tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCac                                           >  2:704043/1‑140 (MQ=255)
                                                                                       tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCac                                           >  2:852775/1‑140 (MQ=255)
                                                                                           gtgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgct                                       <  1:514237/140‑1 (MQ=255)
                                                                                           gtgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgct                                       <  2:853349/140‑1 (MQ=255)
                                                                                            tgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgcta                                      >  1:203246/1‑140 (MQ=255)
                                                                                              ttGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctacc                                    >  2:889334/1‑140 (MQ=255)
                                                                                              ttGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGCCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctacc                                    >  2:621197/1‑140 (MQ=255)
                                                                                                     gCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATa                             <  1:852775/140‑1 (MQ=255)
                                                                                                               tGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCa                   <  1:704043/140‑1 (MQ=255)
                                                                                                                ggCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCAt                  >  2:343678/1‑140 (MQ=255)
                                                                                                                ggCAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGGCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCAt                  >  2:765584/1‑140 (MQ=255)
                                                                                                                  cAGAGTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCg                <  1:334655/140‑1 (MQ=255)
                                                                                                                   agagTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCCCCGCTACCGCTGATACTGTGCGTCATCGc               >  2:110334/1‑140 (MQ=255)
                                                                                                                   agagTCGGAACCGAAAGCAATTT‑‑CCCCCCCGCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTTAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGc               >  1:696497/1‑140 (MQ=255)
                                                                                                                     agTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGccc             >  2:614220/1‑140 (MQ=255)
                                                                                                                     agTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGccc             >  2:904506/1‑140 (MQ=255)
                                                                                                                     agTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGccc             >  1:762808/1‑140 (MQ=255)
                                                                                                                     agTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGccc             >  1:892430/1‑140 (MQ=255)
                                                                                                                     agTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGccc             >  2:671959/1‑140 (MQ=255)
                                                                                                                     agTCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGccc             >  2:6518/1‑140 (MQ=255)
                                                                                                                       tCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAg           >  1:537565/1‑140 (MQ=255)
                                                                                                                       tCGGAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAg           >  2:158598/1‑140 (MQ=255)
                                                                                                                          ggAACCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCt         <  1:158598/140‑1 (MQ=255)
                                                                                                                             aCCGAAAGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACg      <  2:762808/140‑1 (MQ=255)
                                                                                                                                 aaaGCAATTTCCCCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCa  >  1:285138/1‑140 (MQ=255)
                                                                                                                                           ||                                                                                                                                  
AGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑‑CCCCCC‑‑CCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCA  >  NC_000913/1594111‑1594376

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: