Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,630 G→C 12.7% pseudogene (1239/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,6300GC12.7% 127.5 / 6.5 48pseudogene (1239/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (27/10);  new base C (3/3);  total (35/13)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTG  >  NC_000913/1466481‑1466778
                                                                                                                                                     |                                                                                                                                                    
aTTGATATCTCCGGCGACAGCGCATCGGTTGATAACTCTGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACg                                                                                                                                                      <  2:1320913/150‑1 (MQ=17)
aTTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACg                                                                                                                                                      <  2:1093342/150‑1 (MQ=255)
          ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTTAACAACGAAGGCGAGAGCGCAATCACCAACGGTGGCACCGGCACGCAGATAAACg                                                                                                                                            >  2:948807/1‑150 (MQ=255)
                    cgcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTGAACAACGAAGGCGAGAGCGCCATCACCAACGGGGGCACCGGCACTCAGATCAACGGTCACGACGc                                                                                                                                  >  2:1118379/1‑150 (MQ=14)
                     gcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGcc                                                                                                                                 >  1:1283216/1‑150 (MQ=255)
                     gcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGcc                                                                                                                                 >  1:1283228/1‑150 (MQ=255)
                                            cATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                          >  1:448906/1‑150 (MQ=255)
                                              tGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGGGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                                                                                                        >  2:352835/1‑150 (MQ=255)
                                                  cGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCGAACAACAGCGGCAAAACCACCGt                                                                                                    >  1:194263/1‑150 (MQ=255)
                                                   gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTCAAGACGCCACCGCGAAAAACAAAGGAAAAACCACAGtt                                                                                                   >  1:61762/1‑150 (MQ=14)
                                                    tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTATCGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTg                                                                                                  >  1:956548/1‑150 (MQ=255)
                                                                aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                      <  2:1283228/150‑1 (MQ=255)
                                                                aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                      <  2:1283216/150‑1 (MQ=255)
                                                                     aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCGGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAAGGCAAGGATTCCAcc                                                                                 >  1:166286/1‑150 (MQ=255)
                                                                      tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                                                <  1:275998/150‑1 (MQ=255)
                                                                                 gTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAtt                                                                     >  2:1406310/1‑150 (MQ=255)
                                                                                  ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAGCGGTACGGAAATTg                                                                    >  2:427090/1‑150 (MQ=255)
                                                                                     aCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                                                 >  1:200616/1‑150 (MQ=255)
                                                                                       ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTGAAGACGCCACCGCAAAAAACAACGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGCAAATAGCTGGc                                                               >  2:909619/1‑150 (MQ=17)
                                                                                             cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGAGGCACCGGAACGCAGATTAACGGTCATGACGCCACGGCAAAAAACAACGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACGGAAATACCTGGCAATatc                                                         >  2:527134/1‑148 (MQ=17)
                                                                                             cAGGCGGTTGTTAACAACGAAGGCGAGAGCACCATCACCAAAGGCGGCACCGGCACTCAGATAAACGGTAAAGACGCCACCGCTAAAAACAGAGGAAAAACCACAGTTGAAGGAAATGATTCCACCGGTACCAAAATAGCTGGCAATatc                                                         >  2:527143/1‑148 (MQ=12)
                                                                                                   gttgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAg                                                   >  2:1297656/1‑150 (MQ=255)
                                                                                                   gttgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCACTAACGGGAAg                                                   >  2:1297663/1‑150 (MQ=255)
                                                                                                   gttgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAACGGGAAg                                                   >  2:1297679/1‑150 (MQ=255)
                                                                                                       ttAACAACGAAGGCGAGAGCGCCATCACCAACAGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGa                                               >  2:301111/1‑150 (MQ=255)
                                                                                                          ataacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGACATGGCGGGCAATAACGGGACGGTGATTc                                            >  2:2538/3‑150 (MQ=11)
                                                                                                            aacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAg                                          >  2:610480/1‑150 (MQ=255)
                                                                                                            aacGAAGGCGAGAGCGCAATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATCACGCCGCGGCAAACAACTACTTCACAACCACCGTTGAAGGCTAGGATTCCACCGGTACGGAAATTGCTGGCGATCTCGGGCGGGTGATTCAg                                          >  2:837151/1‑150 (MQ=14)
                                                                                                               gAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAc                                       <  1:1379651/150‑1 (MQ=255)
                                                                                                               gAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAc                                       >  1:965733/1‑150 (MQ=255)
                                                                                                                     gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAGGACGCCACGGCAATCAACATAGGCAAAACCACCGTTGTCGGCAAGGATTCCACCGGTACGGAAATTGCTTGCAATGACGGGCAGGTGATTCTGGACGGCGAt                                 >  2:741562/1‑150 (MQ=17)
                                                                                                                        cgctacATCAGCAGCGGTTGCACCGGACCGGAGATTAACTGCGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTTCTGGCAATAACGGGAATGTGATTCAGGACGGCGATCTg                              <  1:80485/145‑1 (MQ=17)
                                                                                                                           gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                           >  2:1019784/1‑150 (MQ=255)
                                                                                                                           cccATCACCAACGGTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                           <  2:632814/149‑1 (MQ=255)
                                                                                                                           accATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                           <  1:802248/149‑1 (MQ=255)
                                                                                                                           accATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                           <  2:500041/149‑1 (MQ=255)
                                                                                                                             cATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                         >  1:558573/1‑150 (MQ=255)
                                                                                                                              aTCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTc                        >  1:581916/1‑150 (MQ=255)
                                                                                                                              aTCACCAACGGAGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTc                        >  1:11102/1‑150 (MQ=255)
                                                                                                                                  ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa                                  <  1:1009601/136‑1 (MQ=255)
                                                                                                                                  ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa                                  <  1:1009489/136‑1 (MQ=255)
                                                                                                                                  ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa                                  >  2:1009489/1‑136 (MQ=255)
                                                                                                                                  ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa                                  >  2:1009601/1‑136 (MQ=255)
                                                                                                                                       ggTGGCACCGGCACGCAGATTACGGGTGATGCCGCCACGGCAAACAACAACGGCAAAACCCCCGTTGACGGCAAGGATGCCCCCGGTCCGGAAAGTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcggc               <  1:2538/150‑1 (MQ=17)
                                                                                                                                       ggTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGAGGGCAATACCGGCAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcggc               <  1:2540/150‑1 (MQ=14)
                                                                                                                                           gCACCGGCACTCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGATGGCAATAACGGGAAGGTGATTCAGGACGGAGATCTGGATGTCAGAGGCGGCGGTc           >  2:1323121/1‑150 (MQ=255)
                                                                                                                                           gCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc           >  2:1323136/1‑150 (MQ=255)
                                                                                                                                               cGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg       >  1:659520/1‑150 (MQ=255)
                                                                                                                                                 gCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTa     <  1:301111/150‑1 (MQ=255)
                                                                                                                                                    cGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg  >  1:690300/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTG  >  NC_000913/1466481‑1466778

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