Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,071,321 A→G 100% intergenic (+110/‑218) yeeP → / → flu pseudogene, CP4‑44 prophage; 50S ribosome‑binding GTPase family;Phage or Prophage Related; putative histone/CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,071,3210AG97.3% 276.3 / ‑2.6 85intergenic (+110/‑218)yeeP/flupseudogene, CP4‑44 prophage; 50S ribosome‑binding GTPase family;Phage or Prophage Related; putative histone/CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (39/43);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGG  >  NC_000913/2071225‑2071414
                                                                                                |                                                                                             
gATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGata                                                                                            >  2:590053/1‑100 (MQ=255)
  tGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATc                                                                                          <  2:1493442/100‑1 (MQ=255)
    ccctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTa                                                                                        <  1:1983091/100‑1 (MQ=255)
    ccctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTa                                                                                        <  1:1636959/100‑1 (MQ=255)
    ccctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTa                                                                                        >  1:1380413/1‑100 (MQ=255)
      ctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTAtt                                                                                      >  2:769112/1‑100 (MQ=255)
         cctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTa                                                                                   <  2:2180006/100‑1 (MQ=255)
              ctctGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatcgat                                                                              <  1:2190936/100‑1 (MQ=255)
               tctGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatcgatc                                                                             <  2:2036505/100‑1 (MQ=255)
                 tGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGt                                                                           >  2:1025467/1‑100 (MQ=255)
                 tGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGt                                                                           >  1:1772123/1‑100 (MQ=255)
                  gAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATATAATAAAACGATCGATATCTATTTTATCGATCGtt                                                                          >  1:1133309/1‑100 (MQ=255)
                  gAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGtt                                                                          <  1:2107590/100‑1 (MQ=255)
                    gtttGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAAAATAAAACGCTCGATATCTATTTTATCGATCGTTta                                                                        <  1:1135260/97‑1 (MQ=255)
                    gatctcTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTta                                                                        <  1:210917/95‑1 (MQ=255)
                    gTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTta                                                                        <  2:1036477/100‑1 (MQ=255)
                    gTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTta                                                                        <  2:852422/100‑1 (MQ=255)
                      cTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGGTCGTTtata                                                                      <  2:2296564/100‑1 (MQ=255)
                      cTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTtata                                                                      <  1:1036909/100‑1 (MQ=255)
                       tGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTTCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTtata                                                                      <  2:1251036/99‑1 (MQ=255)
                        gCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatcgatcgatataaac                                                                    <  2:1503524/100‑9 (MQ=255)
                         cTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTtaagatc                                                                   >  2:210914/1‑95 (MQ=255)
                         cTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcg                                                                   >  2:537685/1‑100 (MQ=255)
                          tCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGGTATCTATTTTATCGATCGTTTATatcga                                                                  <  1:1600728/100‑1 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                 >  2:2387372/1‑100 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                 >  1:2132500/1‑100 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                 >  2:7184/1‑100 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                 >  2:1411531/1‑100 (MQ=255)
                            acaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatagat                                                                 >  2:2086930/1‑99 (MQ=255)
                             caAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgatcg                                                               >  1:1504370/1‑100 (MQ=255)
                               aaaGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATcgatcgat                                                             <  2:1694663/100‑1 (MQ=255)
                                   cACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGc                                                         <  1:204816/100‑1 (MQ=255)
                                    aCTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCt                                                        <  2:1196502/100‑1 (MQ=255)
                                       gTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAACtaat                                                     >  1:860307/1‑100 (MQ=255)
                                        ttttCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCCATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaata                                                    <  1:1686303/100‑1 (MQ=255)
                                        ttttCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaata                                                    <  1:1025478/100‑1 (MQ=255)
                                        ttttCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaata                                                    <  1:1514143/100‑1 (MQ=255)
                                             gTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAAc                                               >  2:1743296/1‑100 (MQ=255)
                                                aCTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCtt                                            <  1:1411546/100‑1 (MQ=255)
                                                  tGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTg                                          <  2:1995901/100‑1 (MQ=255)
                                                  tGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTg                                          <  1:118185/100‑1 (MQ=255)
                                                  tGTCTCTCTTGTCCATGCAATAGCTCAATAATAGAATAAAACGATCGATGTCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTg                                          <  2:365369/100‑1 (MQ=255)
                                                   gTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGt                                         >  1:2365512/1‑100 (MQ=255)
                                                   gTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGt                                         >  1:530535/1‑100 (MQ=255)
                                                    tctcttTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGt                                         <  2:113860/99‑1 (MQ=255)
                                                     ctctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTc                                        <  2:2191468/99‑1 (MQ=255)
                                                      tctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAg                                      <  2:51685/100‑1 (MQ=255)
                                                      tctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAg                                      >  1:1333418/1‑100 (MQ=255)
                                                        tctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTa                                    <  2:726980/100‑1 (MQ=255)
                                                         ctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTaa                                   >  1:1251056/1‑100 (MQ=255)
                                                           tGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACa                                 >  2:1764275/1‑100 (MQ=255)
                                                           tGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACa                                 <  2:1504353/100‑1 (MQ=255)
                                                           tGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCAATTATATAGATCGATAAGCTAATAATAACCTTTGTCAGTAACa                                 >  2:734323/1‑100 (MQ=255)
                                                              ccGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCATTAACATGc                              <  1:1688486/100‑1 (MQ=255)
                                                              ccGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGc                              >  2:343206/1‑100 (MQ=255)
                                                              ccGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGc                              <  2:2073443/100‑1 (MQ=255)
                                                               cGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGc                              >  1:2156896/1‑99 (MQ=255)
                                                                gttcAATATATCAATAATATAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGcac                            <  2:2156871/97‑1 (MQ=255)
                                                                gTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACAt                                <  2:665448/96‑1 (MQ=255)
                                                                 tGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGcaca                           <  2:1706758/100‑1 (MQ=255)
                                                                 tGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGcaca                           >  1:1503397/1‑100 (MQ=255)
                                                                    aaTACCTCAATAATAGAATAAAACGATCAGTCTGTAGATCATCGCTCGTCTATAACCATCAACAACCTGATAACAAACTATGGCCTTAACATGCACAGaa                        >  2:1363985/1‑99 (MQ=255)
                                                                     aTAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATa                       >  1:2148660/1‑100 (MQ=255)
                                                                     aTAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATa                       >  1:222785/1‑100 (MQ=255)
                                                                        gCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACTTTTGTTAGTAACATGCACAGATACGt                    >  1:2286002/1‑100 (MQ=255)
                                                                        gCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGt                    >  1:2402349/1‑100 (MQ=255)
                                                                        gCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATAcat                    >  1:210870/1‑98 (MQ=255)
                                                                          tCAATAATAGAATAAAACGATTGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTAc                  <  1:1013505/100‑1 (MQ=255)
                                                                          tCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCCCAGATACGTa                   >  2:1297799/1‑99 (MQ=255)
                                                                          tCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTAc                  <  1:1000627/100‑1 (MQ=255)
                                                                           cAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACa                 >  1:563743/1‑100 (MQ=255)
                                                                           cAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACa                 >  1:434019/1‑100 (MQ=255)
                                                                            aataatAGAATAAAACGATCGATATCTATTTTAGCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAg                <  2:1382711/100‑1 (MQ=255)
                                                                             ataataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGa               >  2:1109857/1‑100 (MQ=255)
                                                                              taataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGaa              >  2:1414413/1‑100 (MQ=255)
                                                                                ataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAg            >  2:581195/1‑100 (MQ=255)
                                                                                ataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAg            >  1:1798855/1‑100 (MQ=255)
                                                                                ataGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACATACAGAAAg            <  2:210867/100‑1 (MQ=255)
                                                                                     aTAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACAtt       >  1:1993943/1‑100 (MQ=255)
                                                                                     aTAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACAt        >  1:856903/1‑99 (MQ=255)
                                                                                        aaaCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAg    <  1:2399257/100‑1 (MQ=255)
                                                                                         aaCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAgg   <  1:1533053/100‑1 (MQ=255)
                                                                                          aCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAggg  >  1:2063526/1‑100 (MQ=255)
                                                                                          aCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAggg  <  2:1244388/100‑1 (MQ=255)
                                                                                          aCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAggg  <  2:942992/100‑1 (MQ=255)
                                                                                                |                                                                                             
GATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGG  >  NC_000913/2071225‑2071414

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