Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,132 (C)8→7 pseudogene (3955/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,1250C.89.2% 121.9 / 11.8 37pseudogene (3962/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base . (17/16);  total (19/18)
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.93e-01

CATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGA  >  NC_000913/1593991‑1594234
                                                                                                                                      |                                                                                                             
catcaGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccc                                                                                                          <  2:12589/140‑1 (MQ=255)
     ggCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccccaggc                                                                                                     >  1:295220/1‑136 (MQ=255)
          ttGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATaa                                                                                                >  1:629899/1‑140 (MQ=255)
           tGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACg                                                                                              >  2:244127/1‑140 (MQ=255)
               tCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCa                                                                                          >  2:52423/1‑140 (MQ=255)
                ccGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCat                                                                                         >  2:348578/1‑140 (MQ=255)
                ccGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCat                                                                                         <  1:52423/140‑1 (MQ=255)
                   gTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATc                                                                                      <  1:348578/140‑1 (MQ=255)
                   gTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATc                                                                                      >  2:360753/1‑140 (MQ=255)
                            ccGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAgt                                                                             >  1:362171/1‑140 (MQ=255)
                             cGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCat                                                                                         <  2:460505/127‑1 (MQ=255)
                             cGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCat                                                                                         >  1:460505/1‑127 (MQ=255)
                             cGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAgtg                                                                            >  1:438061/1‑140 (MQ=255)
                                        aGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                                                 <  2:362171/140‑1 (MQ=255)
                                         gTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTccc                                                                <  1:529136/140‑1 (MQ=255)
                                          tACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCg                                                               >  1:446176/1‑140 (MQ=255)
                                               tCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATc                                                          >  2:324824/1‑140 (MQ=255)
                                                ccccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATca                                                         <  1:360753/140‑1 (MQ=255)
                                                  ccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcaca                                                       <  1:466706/140‑1 (MQ=255)
                                                   cGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                      <  1:588142/140‑1 (MQ=255)
                                                   cGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                      <  2:319583/140‑1 (MQ=255)
                                                    ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACAcc                                                     >  2:111100/1‑140 (MQ=255)
                                                     gTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                                    <  2:438061/140‑1 (MQ=255)
                                                     gTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                                    <  1:244127/140‑1 (MQ=255)
                                                      tttCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                                    <  2:629899/140‑1 (MQ=255)
                                                      tttCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGc                                                   >  2:606208/1‑140 (MQ=255)
                                                          cccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                               <  1:606208/140‑1 (MQ=255)
                                                          cccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                               <  1:324824/140‑1 (MQ=255)
                                                          cccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                               <  2:446176/140‑1 (MQ=255)
                                                                      ttCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGa                                          >  1:385942/1‑134 (MQ=255)
                                                                      ttCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGa                                          <  2:385942/134‑1 (MQ=255)
                                                                        cATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtgt                                 >  1:206652/1‑140 (MQ=255)
                                                                              aaCCGTGCCACCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGAc                           >  1:93816/1‑140 (MQ=255)
                                                                                ccGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACtt                         >  1:146801/1‑140 (MQ=255)
                                                                                   tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGata                                        <  2:30355/122‑1 (MQ=255)
                                                                                   tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGata                                        >  1:30355/1‑122 (MQ=255)
                                                                                   tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc                      <  2:146801/140‑1 (MQ=255)
                                                                                   tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGCTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTccc                      >  1:276280/1‑140 (MQ=255)
                                                                                                       ccgccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGa  >  2:223542/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                             
CATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGA  >  NC_000913/1593991‑1594234

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