Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,071,321 A→G 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.0% 222.4 / ‑2.3 66intergenic (+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 (1/1);  new base G (30/34);  total (31/35)
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.55e-01

GATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAAC  >  NC_000913/2071225‑2071417
                                                                                                |                                                                                                
gATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGata                                                                                               >  1:981673/1‑100 (MQ=255)
gATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGata                                                                                               <  2:298845/100‑1 (MQ=255)
  tGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGatatn                                                                                             <  2:572903/100‑2 (MQ=255)
  tGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATc                                                                                             >  1:1215244/1‑100 (MQ=255)
  tGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTGCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATc                                                                                             <  1:405080/100‑1 (MQ=255)
   gCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCt                                                                                            >  1:469160/1‑100 (MQ=255)
          ctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTat                                                                                     <  2:62012/100‑1 (MQ=255)
           tGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATACAACGATCGATATCTATTTTatc                                                                                    <  2:755988/100‑1 (MQ=255)
           tGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatc                                                                                    <  2:641411/100‑1 (MQ=255)
            gACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatcg                                                                                   <  1:832986/100‑1 (MQ=255)
            gACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatcg                                                                                   >  2:971022/1‑100 (MQ=255)
              ctctGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatcgat                                                                                 <  2:460699/100‑1 (MQ=255)
                ctGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatcgatcg                                                                               >  1:1083005/1‑100 (MQ=255)
                       tGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTatat                                                                        <  1:314859/100‑1 (MQ=255)
                       tGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTatat                                                                        <  2:1280182/100‑1 (MQ=255)
                       tGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTatat                                                                        >  1:1096384/1‑100 (MQ=255)
                       tGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTAATCGATCGTTtata                                                                         <  2:686299/99‑1 (MQ=255)
                        gCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatc                                                                       >  1:1156839/1‑100 (MQ=255)
                        gCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatc                                                                       >  1:1155960/1‑100 (MQ=255)
                         cTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcg                                                                      >  1:1191941/1‑100 (MQ=255)
                         cTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcg                                                                      >  2:422446/1‑100 (MQ=255)
                           nacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                    >  2:733983/2‑100 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                    >  1:469473/1‑100 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                    >  2:337120/1‑100 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                    >  2:750033/1‑100 (MQ=255)
                            acaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgatc                                                                   <  1:137266/100‑1 (MQ=255)
                                  gCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAg                                                             <  2:976592/100‑1 (MQ=255)
                                  gCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAg                                                             <  2:1291469/100‑1 (MQ=255)
                                   cACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGc                                                            >  1:442393/1‑100 (MQ=255)
                                    aCTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCt                                                           <  1:237479/100‑1 (MQ=255)
                                    aCTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCt                                                           <  1:62086/100‑1 (MQ=255)
                                       gTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtcat                                                        >  2:1271193/1‑100 (MQ=255)
                                       gTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaat                                                        <  2:1191930/100‑1 (MQ=255)
                                        ttttCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaata                                                       >  2:641742/1‑100 (MQ=255)
                                        ttttCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaata                                                       <  2:1223780/100‑1 (MQ=255)
                                                aCTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCtt                                               <  2:1068510/100‑1 (MQ=255)
                                                   gTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGt                                            >  2:62397/1‑100 (MQ=255)
                                                   gTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGt                                            >  2:484870/1‑100 (MQ=255)
                                                      tatctTGTCCGTGGAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAg                                         >  1:884443/3‑100 (MQ=255)
                                                          ttGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAAc                                     <  1:1031458/100‑1 (MQ=255)
                                                           tGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACa                                    >  1:981459/1‑100 (MQ=255)
                                                            gTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACAt                                   >  1:1223438/1‑100 (MQ=255)
                                                                  gCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAg                             <  2:469469/100‑1 (MQ=255)
                                                                   cAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGa                            >  1:500099/1‑100 (MQ=255)
                                                                   cAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGCCAGTAACATGCACAGa                            >  2:102922/1‑100 (MQ=255)
                                                                    aaTAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGAt                           <  1:792211/100‑1 (MQ=255)
                                                                     aTAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGAt                           >  1:486927/1‑99 (MQ=255)
                                                                        gCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGt                       >  1:176/1‑100 (MQ=255)
                                                                         cTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTa                      <  1:484874/100‑1 (MQ=255)
                                                                         cTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTa                      <  1:1208307/100‑1 (MQ=255)
                                                                          tCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTAc                     <  2:1056077/100‑1 (MQ=255)
                                                                           cAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACa                    <  2:1179689/100‑1 (MQ=255)
                                                                           cAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACa                    <  2:340247/100‑1 (MQ=255)
                                                                              taataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGaa                 <  1:669426/100‑1 (MQ=255)
                                                                              taataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGaa                 <  2:1156828/100‑1 (MQ=255)
                                                                               aataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGaaa                >  2:1243217/1‑100 (MQ=255)
                                                                               aataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGaaa                <  2:1082995/100‑1 (MQ=255)
                                                                                 taGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGa              >  1:269696/1‑100 (MQ=255)
                                                                                 taGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGa              <  1:746469/100‑1 (MQ=255)
                                                                                  aGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGan             <  2:868324/100‑2 (MQ=255)
                                                                                  aGAATAAAAAGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAAAAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGAc             <  2:462828/100‑1 (MQ=255)
                                                                                   gAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACa            >  2:343136/1‑100 (MQ=255)
                                                                                     aTAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACAtt          <  2:877602/100‑1 (MQ=255)
                                                                                         aaCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAgg      <  2:1005273/100‑1 (MQ=255)
                                                                                          aCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAggg     >  2:160867/1‑100 (MQ=255)
                                                                                           cgatcgatATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGa    >  1:1152246/1‑100 (MQ=255)
                                                                                             aTCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGaac  >  2:570936/1‑100 (MQ=255)
                                                                                                |                                                                                                
GATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAAC  >  NC_000913/2071225‑2071417

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