Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,280,104 A→C pseudogene (2297/2532 nt) yejO ← autotransporter YejO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,280,1040AC94.6% 115.4 / ‑1.2 37pseudogene (2297/2532 nt)yejOautotransporter YejO
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (19/16);  total (20/17)
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.63e-01

GACCTTATACCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCAGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGC  >  NZ_CP009273/2279975‑2280236
                                                                                                                                 |                                                                                                                                    
gACCTTATACCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTAc                                                                                                                             >  1:291437/1‑139 (MQ=255)
 aCCTTATACCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTAcc                                                                                                                            >  2:50391/1‑139 (MQ=255)
     tataCCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGttt                                                                                                                        >  2:388091/1‑139 (MQ=255)
            gCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGt                                                                                                                 >  2:125296/1‑139 (MQ=255)
            gCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGCCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGt                                                                                                                 >  2:99721/1‑139 (MQ=255)
             ccTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTc                                                                                                                >  1:258995/1‑139 (MQ=255)
             ccTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTc                                                                                                                >  2:145815/1‑139 (MQ=255)
                 aTAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCTTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCAtg                                                                                                            >  1:201122/1‑139 (MQ=255)
                         cACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCt                                                                                                    <  2:306967/139‑1 (MQ=255)
                           cTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTAc                                                                                                                             <  2:431723/112‑1 (MQ=255)
                           cTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTAc                                                                                                                             >  1:431723/1‑112 (MQ=255)
                                   aGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGtt                                                                                          >  2:396194/1‑139 (MQ=255)
                                        catcatTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGc                                                                                     >  2:232592/1‑139 (MQ=255)
                                         atcatTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCt                                                                                    >  2:154211/1‑139 (MQ=255)
                                             ttattaAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAg                                                                                <  1:232592/139‑1 (MQ=255)
                                             ttattaAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAg                                                                                <  1:125296/139‑1 (MQ=255)
                                              tattaAATTTGCCACCTTCATTAACTTTCCCCTGGTTAGAATCGGCTTATCCCTGACCCACGGCCGTTTACAGCCAGGTTCCCCGCGTCGAACCGTTTTGAAGGTCCTTGTGTTAGCTTACCGCCGTTCCCGCTTCAg                                                                                <  1:172690/138‑1 (MQ=255)
                                               attaAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCCCGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAgc                                                                               <  2:176199/138‑1 (MQ=255)
                                                  aaaTTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGcgcg                                                                            <  1:388091/138‑1 (MQ=255)
                                                   aaTTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGt                                                                           >  2:201298/1‑138 (MQ=255)
                                                        gCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATc                                                                     >  2:203956/1‑139 (MQ=255)
                                                        gCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCAGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATc                                                                     >  1:289350/1‑139 (MQ=255)
                                                         ccATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATcc                                                                    <  1:50391/139‑1 (MQ=255)
                                                                                 ttAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCa                                            <  2:258995/139‑1 (MQ=255)
                                                                                 ttAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCAGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCa                                            <  2:289350/139‑1 (MQ=255)
                                                                                    gAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCATg                                          <  2:291437/138‑1 (MQ=255)
                                                                                    gAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCATg                                          <  2:421252/138‑1 (MQ=255)
                                                                                         ggCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGttt                                    >  2:19004/1‑139 (MQ=255)
                                                                                          gCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCTTGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATcc                                                                    <  1:389256/106‑1 (MQ=255)
                                                                                          gCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATcc                                                                    >  2:389256/1‑106 (MQ=255)
                                                                                           cGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATc                                                                     >  2:144591/1‑104 (MQ=255)
                                                                                           cGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGa                                  >  2:72386/1‑139 (MQ=255)
                                                                                           cGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTGGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATc                                                                     <  1:144591/104‑1 (MQ=255)
                                                                                             tATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATa                                <  1:32968/139‑1 (MQ=255)
                                                                                                                cTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATc             >  2:437480/1‑139 (MQ=255)
                                                                                                                cTTTCAGCCAGGGTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATc             <  1:437480/139‑1 (MQ=255)
                                                                                                                           ggTTCCCGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGc  <  1:203956/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                    
GACCTTATACCAGCCTGATAAACGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCAGCGTCGTACCGTTTTGCAGGTCCATGTGATAGCTTACCGCCGTTCCCGCTTCAGCGCGTAATATCCGGGTATTTCCCACATCCGCGCGCATGCCGTTTGATAACGTGTAGTCCTGACCATCTGTGGTAAAGC  >  NZ_CP009273/2279975‑2280236

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