Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,452,339 A→G intergenic (‑3/+17) yfcS ← / ← yfcU putative periplasmic pilin chaperone/pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,452,3390AG97.5% 282.2 / ‑2.3 81intergenic (‑3/+17)yfcS/yfcUputative periplasmic pilin chaperone/pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (38/41);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATG  >  NC_000913/2452242‑2452434
                                                                                                 |                                                                                               
ggCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATc                                                                                               <  1:608000/100‑1 (MQ=255)
ggCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATc                                                                                               <  1:285632/100‑1 (MQ=255)
 gCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCa                                                                                              >  1:1121007/1‑100 (MQ=255)
 gCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCa                                                                                              >  1:1441743/1‑100 (MQ=255)
 gCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCa                                                                                              >  1:1315754/1‑100 (MQ=255)
  cGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCAt                                                                                             >  2:1227119/1‑100 (MQ=255)
  cGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCAt                                                                                             >  1:170654/1‑100 (MQ=255)
     tACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTcc                                                                                          <  1:1313109/100‑1 (MQ=255)
      aCAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCt                                                                                         <  1:625364/100‑1 (MQ=255)
      aCAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCt                                                                                         >  2:238971/1‑100 (MQ=255)
      aCAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCt                                                                                         <  1:792261/100‑1 (MQ=255)
         gACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTaa                                                                                      >  1:658223/1‑100 (MQ=255)
              cAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACg                                                                                 <  2:608901/100‑1 (MQ=255)
              cAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACg                                                                                 >  2:1125514/1‑100 (MQ=255)
               aGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGt                                                                                >  2:879220/1‑100 (MQ=255)
                gggCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGtt                                                                               >  1:708083/1‑100 (MQ=255)
                 ggCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTa                                                                              >  1:1516489/1‑100 (MQ=255)
                   cAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTAtt                                                                            >  1:1439138/1‑100 (MQ=255)
                      aCAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGt                                                                         <  1:124220/100‑1 (MQ=255)
                       cAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGtt                                                                        <  2:1389056/100‑1 (MQ=255)
                       cAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGtt                                                                        >  2:879709/1‑100 (MQ=255)
                       cAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGtt                                                                        >  2:1356049/1‑100 (MQ=255)
                         aaCTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAg                                                                      <  1:521099/100‑1 (MQ=255)
                            tGGTAGCGGTAAGCAGTAACACCAGCGTCATAGCGCCCAATTTGGTAGAGCCATGTAAATCAGACATGGGTCATTCCTTAAGCTCGTTATTGGTCAGCt                                                                    <  2:261339/99‑1 (MQ=255)
                             ggTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTga                                                                  >  1:29669/1‑100 (MQ=255)
                             ggTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTga                                                                  >  2:899664/1‑100 (MQ=255)
                                 gcgcTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTgaaaca                                                              >  2:285721/1‑100 (MQ=255)
                                  cgcTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTgaaacag                                                             <  2:1441733/100‑1 (MQ=255)
                                   gcTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTgaaacaga                                                            >  1:1018555/1‑100 (MQ=255)
                                     tAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTgaaacagaaa                                                          <  2:223457/100‑1 (MQ=255)
                                       aGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTgaaacagaaaca                                                        >  1:631675/1‑100 (MQ=255)
                                        gaaGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGaaacagaaacag                                                       >  1:1057110/3‑100 (MQ=255)
                                         cAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAgg                                                      >  1:695761/1‑100 (MQ=255)
                                           gTAACGCCATCGTCATAGCACCCAATTTGGCAGAGCAAATTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAgggg                                                    <  2:1429475/100‑1 (MQ=255)
                                           gTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAgggg                                                    <  1:79201/100‑1 (MQ=255)
                                           gTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAgggg                                                    >  2:1423791/1‑100 (MQ=255)
                                            tAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAggggg                                                   <  2:1439128/100‑1 (MQ=255)
                                            tAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAAACAGACATGGGTCATTCCTTACTCACGTTATTGGTTAGCTGAAACAGAAACAggggg                                                   <  1:392813/100‑1 (MQ=255)
                                             aaCGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGa                                                  >  1:608106/1‑100 (MQ=255)
                                                 ccAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCgg                                              <  2:836439/100‑1 (MQ=255)
                                                  ttccgatctAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggt                                             <  2:289867/92‑1 (MQ=255)
                                                    gCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtgg                                           >  1:1559376/1‑100 (MQ=255)
                                                     cGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtggg                                          <  2:232197/100‑1 (MQ=255)
                                                     cGTCATAACACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtggg                                          <  1:821247/100‑1 (MQ=255)
                                                      gTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtgggt                                         >  1:875825/1‑100 (MQ=255)
                                                       tCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtgggtg                                        <  2:1567308/100‑1 (MQ=255)
                                                       tCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtgggtg                                        <  1:682311/100‑1 (MQ=255)
                                                         aTAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTc                                      <  1:1054469/100‑1 (MQ=255)
                                                           aGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtagatcggaa                                    >  1:289868/1‑91 (MQ=255)
                                                           aGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCAc                                    <  1:1227126/100‑1 (MQ=255)
                                                           aGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCAc                                    <  1:101078/100‑1 (MQ=255)
                                                            gCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCAcc                                   <  2:170654/100‑1 (MQ=255)
                                                            gCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCAcc                                   >  1:1185385/1‑100 (MQ=255)
                                                             cACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCt                                  >  2:1419641/1‑100 (MQ=255)
                                                                 cAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTg                              <  2:1161132/100‑1 (MQ=255)
                                                                 cAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTg                              <  2:29669/100‑1 (MQ=255)
                                                                 cAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTg                              <  1:1423801/100‑1 (MQ=255)
                                                                    tttGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGt                           <  1:251912/100‑1 (MQ=255)
                                                                     ttGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTc                          <  2:564614/100‑1 (MQ=255)
                                                                      tGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCt                         <  2:695758/100‑1 (MQ=255)
                                                                        gCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGc                       >  2:1097019/1‑100 (MQ=255)
                                                                         nagagCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCt                      >  2:71561/2‑100 (MQ=255)
                                                                         cAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCt                      >  2:71813/1‑100 (MQ=255)
                                                                         cAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCt                      <  1:1356058/100‑1 (MQ=255)
                                                                         cAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCt                      >  2:486004/1‑100 (MQ=255)
                                                                          agagCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTc                     >  1:938832/1‑100 (MQ=255)
                                                                           gagCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTcc                    <  1:1242109/100‑1 (MQ=255)
                                                                           gagCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTcc                    >  2:1403453/1‑100 (MQ=255)
                                                                            agCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAACCAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTcc                    <  1:704566/99‑1 (MQ=255)
                                                                            agCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCt                   >  2:327450/1‑100 (MQ=255)
                                                                             gCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTcc                    <  1:626517/98‑1 (MQ=255)
                                                                               agaGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGaa                >  1:1062789/3‑100 (MQ=255)
                                                                               aaaGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGaa                >  2:1140633/1‑100 (MQ=255)
                                                                                aaGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAAt               <  2:986718/100‑1 (MQ=255)
                                                                                aaGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAAt               >  2:824662/1‑100 (MQ=255)
                                                                                aaGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAAt               <  2:81301/100‑1 (MQ=255)
                                                                                aaGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAAt               <  2:1315745/100‑1 (MQ=255)
                                                                                 aGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATg              <  1:1125521/100‑1 (MQ=255)
                                                                                  gTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGa             <  2:209129/100‑1 (MQ=255)
                                                                                    aaaTCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACa           <  2:1185378/100‑1 (MQ=255)
                                                                                             aTGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATg  <  1:674958/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATG  >  NC_000913/2452242‑2452434

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